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Preface 


The  purpose  of  this  research  was  to  determine  the  minimum  performance  re¬ 
quirements  for  a  ground-based,  infrared  LAJJAR  designed  to  detect  deep-space  satel¬ 
lites,  and  to  select  a  candidate  sensor  design  using  available  technology,  i  he  research 
shows  that  a  design  capable  of  high  accuracy  and  resolution  is  feasible  using  current 
technology,  with  greatly  improved  sensor  availability  compared  to  the  Ground-Based 
Electro-Optical  Deep-Space  Surveillance  (GEODSS)  system  currently  in  use. 
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search.  I  also  thank  Major  ihomas  Kelso,  my  co-advisor,  for  assistance  with  IATjtX, 
and  Michael  Wichura  for  inventing  the  PjCTeX  macros.  Finally,  I  wish  to  thank  my 
family  for  their  perseverance  throughout  this  1.500- hour  ‘academic  exercise’. 
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i  Abstract 

This  research  determines  the  minimum  performance  requirements  of  a  ground- 
based  infrared  LADAR  designed  to  detect  deep-space  satellites,  and  presents  a  can¬ 
didate  sensor  design  based  on  current  technology.  The  research  examines  LADAR 
techniques  and  detection  methods  to  determine  the  optimum  LADAR  configuration, 
and  then  assesses  the  effects  of  atmospheric  transmission,  background  radiance,  and 
turbulence  across  the  infrared  region  to  find  the  optimum  laser  wavelengths.  Diffrac¬ 
tion  theory  is  then  used  in  a  parametric  analysis  of  the  transmitted  laser  beam  and  re¬ 
ceived  signal,  using  a  Cassegrainian  telescope  design  and  heterodyne  detection.  The 
effects  of  beam  truncation  and  obscuration,  heterodyne  misalignment,  off-boresight 
detection,  and  image-pixel  geometry  are  also  included  in  the  analysis. 

The  derived  equations  are  then  used  to  assess  the  feasibility  of  several  candi¬ 
date  designs  under  a  wide  range  of  detection  conditions  including  daylight  operation 
through  cirrus.  The  results  show  that  successful  detection  is  theoretically  possi¬ 
ble  under  most  conditions  by  transmitting  a  high-power  frequency-modulated  pulse 

_L *,  ’  '  \  ■ 

train  from  an  isotopic  l_3C02  laser  radiating  at  11.17 ^fn,  and  utilizing  post-detection 
integration  and  pulse  compression  techniques.  1  >  *  . 
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GROUND-BASED  DEEP-SPACE  LADAR 
FOR  SATELLITE  DETECTION: 

A  PARAMETRIC  STUDY 


I.  Introduction 


Background 

The  mission  responsibilities  of  US  Space  Command  lie  in  the  areas  of  space 
operations  and  aerospace  defense  (81:36).  In  support  of  these  responsibilities,  US 
Space  Command  controls  a  comprehensive  network  of  space  surveillance  sensors 
known  as  the  Space  Surveillance  Network  (SSN).  Although  some  of  the  28  SSN 
sensors  are  capable  of  detecting  objects  in  deep-space  (from  5,000  km  to  about 
40,000  km  altitude),  the  only  dedicated  deep-space  capability  lies  with  the  Ground- 
Based  Electro-Optical  Deep-Space  Surveillance  System  (GEODSS)  (81:42:  62:18). 
Since  satellite  detection  by  GEODSS  is  achieved  through  the  detection  of  sunlight 
reflected  off  the  satellite,  GEODSS  operation  relies  on  clear  weather,  night-time 
conditions.  This  severe  limitation  means  that  each  of  the  four  currently  existing 
GEODSS  sensors  is  available  for  an  average  of  only  14%  of  the  time  for  surveillance 
tasks  —  a  shortfall  which  has  contributed  to  the  widening  gap  between  the  estimated 
population  of  deep-space  objects,  and  the  SSN-monitored  population  (25:1).  The 
current  escalation  in  numbers  of  deep-space  satellite  systems,  coupled  with  rising 
Soviet  activity  in  geosynchronous  orbit,  and  continued  Soviet  development  of  ASAT 
capability,  now  mandates  increased  deep-space  surveillance  (74:252;  58:19). 

o 

The  majority  of  plans  for  deep-space  capability  improvements  (30:Sec  V,2; 
58:38)  have  so  far  been  aimed  at  increasing  GEODSS  track  rates  (number  of  dif¬ 
ferent  satellite  detections  per  hour),  while  no  reports  were  found  which  investigate 
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improved  sensor  availability.  Attempts  to  improve  GEODSS  availability  would  have 
little  utility,  however,  due  to  the  inherent  limitations  of  passive  sensing  at  visible 
wavelengths.  Clearly,  an  alternative  detection  method  must  be  used  to  gain  useful 
improvements  in  sensor  availability  under  diverse  weather  and  lighting  conditions. 
Therefore,  active  detection  techniques  must  be  introduced.  These  techniques  could 
be  implemented  for  deep-space  surveillance  using  conventional  microwave  radars, 
but  the  large  size,  high  cost,  and  poor  search  capability  of  radars  preclude  them 
from  serious  consideration.  An  attractive  alternative  to  microwave  radar  is  infrared 
LADAR  (LAser  ra.DAR),  which  uses  an  infrared  (IR)  laser  as  the  active  source. 
LADAR  advantages  over  microwave  radar  include  smaller  size,  superior  accuracy, 
and  improved  resolution  (29:3;  12:1526;  73:37-2).  Benefits  of  using  IR  over  visible 
wavelengths  include  operation  in  daylight  and  successful  detection  under  a  diverse 
range  of  meteorological  conditions. 

Although  LADAR  technology  has  improved  steadily  throughout  the  past  two 
decades,  the  application  of  LADAR  techniques  to  satellite  detection  has  been  lim¬ 
ited.  ‘Enhanced’  satellites,  equipped  with  retro-directive  corner-cube  reflectors  for 
enhanced  laser  returns,  have  been  successfully  detected  with  LADAR  on  several 
occasions  (38:1742;  54:77;  20:398).  However,  the  only  LADAR  routinely  used  for 
‘unenhanced’  satellite  detection  is  the  Firepond  LADAR  operated  by  MIT  Lincoln 
Laboratory  at  Westford,  Massachusetts  (26:329;  77:148;  76:1).  The  major  Firepond 
limitation  for  deep-space  application  is  that  ‘unenhanced’  detection  has  been  lim¬ 
ited  to  medium  altitudes  (below  10,000  km)  (76:1).  No  evidence  has  been  found  of 
the  successful  demonstration  of  an  operational  LADAR  for  detecting  ‘unenhanced’ 
deep-space  satellites  under  adverse  weather,  daylight  conditions. 

Specific  Problem 

The  viability  of  an  alternative  ground-based  deep-space  senior  with  high  avail¬ 
ability  requires  investigation.  The  purpose  of  this  research,  therefore,  is  to  analyt- 
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ically  determine  the  minimum  theoretical  performance  requirements  for  a  ground- 
based,  / R  LADAR  for  deep-space  satellite  detection,  and  to  select  a  candidate  sen¬ 
sor  design  by  matching  minimum  performance  parameters  with  available  technology. 
Where  available  technology  cannot  meet  the  required  performance  specifications,  the 
corresponding  technological  shortfalls  will  be  highlighted. 

Scope  and  Limitations 

The  general  research  approach  is  limited  to  parametric  analysis  of  a  LADAR 
system  which  is  capable  of  search,  acquisition,  and  tracking  of  a  deep-space  satellite 
to  provide  range,  direction,  and  velocity  information.  Specifically,  the  research  ana¬ 
lyzes  the  characteristics  of  the  laser,  optics,  detector,  target,  background  radiance, 
and  atmospheric  transmission.  The  effects  on  performance  of  weather,  turbulence, 
and  daylight  are  also  included,  along  with  a  review  of  search  and  detection  tech¬ 
niques,  and  scanning  strategies. 

The  analysis  does  not  consider  LADAR  mechanical  construction,  or  detailed 
design  of  signal-processing  circuitry  beyond  sub-system  level.  Furthermore,  detection 
of  satellites  below  5,000  km  altitude  and  above  40,000  km  altitude  is  not  considered. 

Methodology 

Investigation  of  the  research  problem  is  divided  into  four  major  subject  ar¬ 
eas:  (1)  LADAR  techniques,  (2)  the  infrared  environment,  (3)  parametric  analysis, 
and  (4)  candidate  system  design.  Following  a  review  of  background  literature  in 
Chapter  II,  research  results  corresponding  to  the  four  subject  areas  are  presented  in 
Chapters  III  to  VI,  inclusive.  Research  conclusions  and  recommendations  are  then 
presented  in  Chapter  VII.  The  sub-objectives  of  each  chapter  are  outlined  in  the 
following  paragraphs. 


3 


LADAR  Techniques.  This  research  phase  reviews  available  LADAR  techniques 
for  long-range  high-resolution  LADAR.  The  four  major  research  sub-objectives  are 
as  follows: 

1.  Target  Characteristics.  Determine  the  approximate  size  of  deep-space  satel¬ 
lites.  and  their  typical  IR  reflectivity  characteristics. 

2.  Detection  Methods.  Review  the  application  of  coherent  and  incoherent  LADAR 
detection  schemes,  including  the  use  of  high-resolution  focal-plane  detector 
arrays. 

3.  LADAR  Technology.  Review  the  performance  characteristics  of  currently  avail¬ 
able  lasers,  optics,  and  detectors  suitable  for  LADAR  applications. 

4.  Target  Acquisition.  Determine  a  suitable  search  scheme  for  deep-space  target 
acquisition. 

The  Infrared  Environment.  The  second  research  phase  examines  the  major 
environmental  characteristics  which  affect  LADAR  reception  and  IR  propagation  in 
the  atmosphere.  The  three  major  research  sub-objectives  are  as  follows: 

1.  Background  Radiance.  Determine  the  level  of  background  solar  and  thermal 
sky-radiance  across  the  IR  region. 

2.  Atmospheric  Transmission.  Determine  the  transmittance  characteristics  of  IR 
laser  beams  through  typical  atmospheric  aerosols  for  various  weather  conditions 
and  slant  paths. 

3.  Other  Propagation  Effects.  Assess  the  influence  of  atmospheric  turbulence  and 
high-power  beam  effects  on  laser  beam  propagation  through  the  atmosphere. 

Parametric  Analysis.  The  third  phase  in  the  solution  process  involves  the  de- 

o 

velopment  of  all  design  parameter  equations  necessary  to  predict  theoretical  LADAR 
performance.  Sub-objectives  include  the  optimization  of  laser,  optics,  and  detector 
characteristics  to  achieve  maximum  system  efficiency. 
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Candidate  Design.  The  final  phase  of  research  involves  translation  of  the  de¬ 
rived  LADAR  parametric  equations  into  a  candidate  design  based  on  the  perfor¬ 
mance  capabilities  of  current  technology.  The  primary  sub-objective  is  to  determine 
whether  the  necessary  laser,  detector,  optics,  and  processing  technology  is  currently 
available  to  match  the  specified  performance  requirements. 
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II.  Literature  Review 


Introduction 

This  chapter  reviews  background  literature  pertinent  to  this  research.  In  par¬ 
ticular.  the  review  covers  the  topics  of  deep-space  surveillance  objectives,  GEODSS, 
LADAR  history,  Firepond.  and  a  brief  summary  of  current  LADAR  technology. 

Deep-Space  Surveillance  Objectives 

The  objectives  of  the  US  Space  Command  deep-space  sensor  network  are  de¬ 
fined  by  Gaposchkin  as  follows: 

1.  Support  acquisition  of  new  launches;  especially  foreign. 

2.  Monitor,  identify,  and  quantify  orbital  maneuvers. 

3.  Maintain  a  catalogue  of  all  targets  fin  deep-space. 

4.  Support  technique  development  and  target  identification.  (25:1) 

In  fulfillment  of  these  objectives,  the  SSX  tracks  a  variety  of  orbiting  space  ob¬ 
jects,  such  as  active  and  inactive  payloads,  spent  booster-rocket  stages,  and  miscel¬ 
laneous  debris  resulting  from  anomalous  satellite  events,  such  a„s  explosions  (18:306). 
SSN  tracking  data  also  helps  in  planning  space  launches  and  orbital  maneuvers  so 
that  possible  collisions  are  avoided  (18:306). 

Although  SSN  surveillance  capability  appears  comprehensive,  the  resource¬ 
consuming  objective  of  maintaining  a  catalogue  of  deep-space  objects  is  difficult  to 
achieve.  Catalogue  maintenance  places  growing  demands  on  each  sensor,  commensu¬ 
rate  with  the  rapid  growth  in  deep-space  population,  which  is  projected  to  be  2,600 
by  the  year  2000  (47:4).  Furthermore,  even  the  catalogued  objects  require  frequent 
updates  of  their  catalogued  positions  due  to  orbital  maneuvers  and  perturbations. 
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If  a  catalogued  position  becomes  outdated,  the  object  is  'lost.'  A  common  definition 
used  for  a  ’lost"  satellite  is  a  catalogued  item  with  an  orbital  element  set  greater 
than  30  to  60  days  old,  depending  on  the  type  of  object  (25:3).  Table  1  provides  a 
breakdown  of  142  "lost-  satellites,  compiled  by  Taff  for  a  selected  day  in  1986. 

Table  1.  Summary  of  Known,  Lost,  Deep-Space  Satellites  (Reproduced  from  (  7 S : o ) ) 


Object 

Soviet 

\'on-Soviet 

Payloads 

i  Synchronous 

13 

32 

i  . 

i  Semi-Synchronous 

6 

1 

Rocket  Bodies 

Synchronous 

$ 

,S 

1 

Semi -Synchronous 

2 

11 

Debris 

Synchronous 

18 

12 

Non-Synchronous 

29 

2 

TOTALS 

_ 

t>6 

CEODSS 

The  primary  SSN  sensor  for  deep-space  surveillance  is  GEODSS  (Sl:42).  Four 
GEODSS  sites  are  currently  operational,  with  locations  at  Socorro,  New  Mexico: 
Maui,  Hawaii;  Taegu  Air  Base,  South  Korea;  and  Diego  Garcia  in  the  Indian  Ocean. 
A  fifth  site  is  planned  for  location  somewhere  in  the  eastern  Atlantic  region  (47:3). 
Each  GEODSS  site  has  two  primary  telescopes  and  one  auxiliary  telescope  for  sensing 
in  the  0.3-0.86  j.im  band  (47:18).  General  telescope  characteristics  are  listed  in 
Table  2  (63:129;  65:o7-58). 

Each  telescope  mount  contains  a  sidereal  drive,  which  maintains  the  telescope 
field-of-view  fixed  on  the  star  background  so  that  orbiting  satellites  appear  as  streaks 
of  light  on  the  sensed  image  (65:57).  Using  this  method,  objects  with  a  visual 


Table  2.  GEODSS  Telescope  Characteristics 


Parameter 

PRIMARY 

AUXILIARY 

Diameter 

1  m 

40  cm 

Field-of- View 

2. 1  degrees 

6  degrees 

Focal  Length 

218  cm 

76  cm 

Max  Slew  Rate 

2  deg/sec 

15  deg /sec 

Pointing  Accuracy 

1.5  arc-second 

1.5  arc-second 

i  Scan  Rate 
| 

2.400  sq.  deg/hr 

15,000  sq.  deg/hr 

i  Main  Element 

Rite  hey-Chretien 

Schmidt  camera  i 

magnitude  of  15.5  are  routinely  tracked  (33:536),  while  detection  of  fainter  objects 
at  magnitude  16.5  is  possible  under  certain  conditions  (65:57). 

Pn-formanet  Limitations.  The  primary  gauge  of  GEODSS  performance  is 
track  rate  (56:2).  Although  extreme  nightly  variations  occur,  peak  track  rates  were 
reported  in  1987  to  be  between  100  and  150  tracks  per  telescope  per  night  (58:14). 
These  rates  could  be  considerably  improved  through  software  modifications  and  hard¬ 
ware  replacement,  but  are  ultimately  limited  by  ambient  light  and  weather  conditions 
(56:1-2).  Bright-skv  conditions,  such  as  twilight,  or  the  7-dav  period  centered  on  the 
full  moon,  seriously  affect.  GEODSS  capability  (78:7:  58:2).  Additionally,  tracking  is 
not  attempted  when  less  than  1  %  of  the  night  sky  is  clear.  This  weather  condition 
often  exists  for  30%  to  60%  of  available  tracking  time  (47:5). 

.1  railahilitij.  Koehler  reports  seasonal  data  for  each  GEODSS  location  concern¬ 
ing  hours  of  darkness  and  the  probability  of  cloud-free  viewing  (47:60).  When  data  is 
averaged  over  all  seasons  and  sites,  an  average  expected  value  for  hours  of  darkness 
per  night ,  and  probability  of  clear  skies  can  be  derived.  These  figures  are  9.375  and 
0.6.  respectively.  If  a  twilight  period  of  two  hours  per  night  is  considered  unaccept¬ 
able  for  operation,  then  the  average  available  viewing  time,  limited  by  twilight  and 
weather,  is  only  4.5  hours  each  night.  If  91  'bright-sky'  nights  are  then  assumed  to 
occur  each  year  due  to  full  moon  conditions,  average  GEODSS  availability  during 
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hours  of  darkness  is  estimated  to  be  35%.  This  translates  to  an  average  Jj-hour 
availability  of  14%.  These  figures  do  not  account  for  the  additional  limitations  due 
to  solar  illumination  of  the  orbit,  satellite  viewing  geometres,  and  svstem  failures. 
However,  the  calculations  do  support  several  reports  that,  due  to  environmental  re¬ 
strictions  on  operation.  GEODSS  is  not  suited  to  providing  time-critical  surveillance 
support  (58:20;  47:5;  30:Sec  V.2). 

Improvements.  Clearly,  sensor  availability  is  the  limiting  factor  in  improving 
GEODSS  track  rate.  Therefore,  substantial  improvement  in  GEODSS  tolerance  r < » 
adverse  weather  and  high  ambient  light  is  required  to  meet  tasking  objectives.  In¬ 
deed.  one  analysis  of  deep-space  sensor  performance  recommends  that  GEODSS  sen¬ 
sors  be  modified  for  24-nour  capability  by  shifting  operation  to  the  long- wa  veieng’  h 
IR  portion  of  the  spectrum  (30:Sec  V.3).  This  method  could  use  a  passive  1R  sen¬ 
sor  to  detect  object  thermal  radiation,  and  is  attractive  due  to  the  lower  levels  of 
background  sky  radiance  at  IR  wavelengths  when  compared  to  visible  wavelength- 
(80:49).  However,  if  the  object  is  space  debris,  or  an  inactive  payload,  the  tem¬ 
perature  difference  between  object  and  background  may  be  insufficient  for  reliable 
passive  detection. 

A  more  reliable  detection  scheme  could  be  implemented  by  using  an  active  IF. 
laser  source  to  irradiate  the  target.  The  detection  decision  would  be  then  based  on 
the  relative  level  of  received  laser  power  reflected  from  the  target.  Any  additional  IR 
signal  received  due  to  the  object's  thermal  radiation  is,  therefore,  a  bonus'  in  the  ac 
five  scheme  (23).  Processing  of  the  received  target  signal  can  also  provide  additional 
target  information  such  as  range,  velocity,  and  rotation  rate  (26:329;  2S:',)76!. 

Advantages  of  using  IR  wavelengths  over  visible  wavelengths  include  higher 
atmospheric  transmission,  lower  scattering  and  lower  scintillation  (intensitv  !luc- 
tuations)  (12:1526).  When  compared  to  microwave  sources,  IR  laser  sources  pro¬ 
vide  superior  angular,  range,  and  Doppler-frequency  (velocity)  resolution  ( 12:1526: 
73:37-2). 
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LADAR  History 

The  use  of  lasers  for  satellite  detection  was  first  demonstrated  in  1964.  when  a 
Q-switched  ruby  laser  (visible  wavelength  —  694.3  nm)  was  used  to  determine  range 
to  the  Beacon  Explorer  B  satellite  within  a  precision  of  a  few  meters  (20:398).  The 
satellite  was  fitted  with  retro-reflectors  for  enhanced  signal  reflection,  and  detec¬ 
tion  was  achieved  at  night.  Daytime  detection  using  a  ruby  laser  was  subsequently 
demonstrated  in  1967  (38:1743).  Since  then,  at  least  14  reflector-equipped  satellites 
have  been  placed  in  orbit,  providing  reference  points  for  a  variety  of  geodetic  stud¬ 
ies  (20:398:  70:391).  Satellite  detection  at  IR  wavelengths  was  further  assisted  by 
the  launch  of  the  Laser  Geodynamic  Satellite  (LAGEOS)  in  1976.  LAGEOS  was 
specifically  designed  to  support  scientific  applications,  and  consists  of  a  60-cm  di¬ 
ameter  sphere,  studded  with  426  reflectors  (20:398).  Four  of  the  reflectors  are  made 
of  germanium  for  long  wavelength  IR  use,  while  the  remainder  are  for  operation  in 
the  visible  region  (77:152).  For  LADAR  applications,  the  launch  of  LAGEOS  into  a 
stable,  well-defined  orbit  at  5,900  km  altitude  provides  a  reference  point  in  inertial 
space  for  LADAR  calibration  purposes  (20:398;  77:152). 

Although  tracking  of  enhanced  satellites  is  routinely  performed,  evidence  of 
unenhanced  detection  is  scarce.  One  example  system  is  the  Maui  Optical  Tracking 
and  Identification  Facility  (MOTIF),  which  uses  a  ruby  laser  “for  space  object  il¬ 
lumination  or  to  provide  range  information  on  selected  targets"  (1:22).  The  only 
example  found  in  the  literature  of  an  operational  satellite-tracking  IR  LADAR  is  the 
Firepond  LADAR  operated  by  MIT  Lincoln  Laboratory  (31;  32;  10). 

Firepond  LADAR 

The  Firepond  LADAR  project  commenced  development  in  the  late  1960s  with 
primary  emphasis  on  satellite  tracking  and  identification  (26:327).  Drawing  from 
CO2  laser  and  IR  detector  research  at  Lincoln  Laboratory,  a  LADAR  was  developed 
by  the  late  1970s  into  an  “automatic  angle,  range,  and  Doppler  frequency  tracking 
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coherent  laser  radar  operating  at  a  wavelength  of  10.6  fim  ...capable  of  tracking 
and  recording  a  variety  of  data  on  both  enhanced  and  unenhanced  low  earth  orbiting 
targets"  [emphasis  added]  (77:148).  The  characteristics  of  the  Firepond  LADAR  are 
taken  from  Gschwendtner  et  al.  and  summarized  in  Table  3. 

Table  3.  Firepond  LADAR  Characteristics  (Extracted  from  (26:329)) 


Parameter 

Specification 

Wavelength 

10.6  pm 

Aperture 

1.2  meters 

Beamwidth  (A/D) 

s=10  /.trad 

Rayleigh  Range  (D2/A) 

100  km 

Beam  Diameter  (at  1000  km) 

10  meters 

Doppler  Coefficient 

2  kHz/cm/sec 

Peak  Power 

1  to  10  kW 

PRF  Continuously  Variable 

20  to  200  pps 

Pulse  Width  Continuously  Variable 

10  to  1  ms 

Major  Firepond  accomplishments  are  reported  by  Gschwendtner  et  al.  as  fol¬ 
lows: 

1.  Development  of  a  high-frequency-stability,  long-range,  coherent  IR  radar. 

2.  Demonstration  of  monopulse  angle  tracking  of  enhanced  and  unenhanced  satel¬ 
lite  targets. 

3.  Determination  of  target  center-of-gravity  motion  and  rotational  dynamics  using 
Doppler  measurements  capability. 

4.  Generation  of  one-dimensional  pseudo-images  from  narrowband 
time-Doppler-intensity  data.  (26:329) 

Current  LADAR  Technology 

In  principle,  an  increase  in  LADAR  detection  range  from  low  altitude  to  deep- 
space  appears  to  affect  only  the  free-space  Moss’  term  in  the  LADAR  range  equation. 
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However,  implementation  is  not  so  simple  due  to  limitations  in  laser  transmitter 
power,  and  the  consequent  need  for  a  very  small  divergence  laser  beam  to  pro¬ 
vide  high  irradiance  at  long  range  (4:53).  This  narrow-beam  requirement  increases 
the  need  for  an  accurate  a  priori  orbital  prediction  for  the  object,  so  that  acqui¬ 
sition/search  time  is  reduced.  The  narrow-beam  requirement  also  means  that  very 
precise  pointing  and  tracking  of  the  laser/telescope  combination  is  essential,  so  that 
the  object  remains  in  the  telescope  field-of-view. 

Technological  research  and  development  has  produced  many  LADAR  improve¬ 
ments  during  the  last  decade.  Advances  in  laser  technology  have  been  considerable, 
with  current  IR  lasers  exhibiting  high  power,  stable  operation,  and  narrow  beam 
divergence  (31:  35:  60;  66).  The  optics  and  detector  components  needed  to  handle 
high  levels  of  laser  power  and  provide  long-range  resolution  have  also  undergone  ex¬ 
tensive  development  (7;  11;  3S;  64).  Finally,  the  acquisition,  tracking,  and  pointing 
sub-systems  have  been  adapted  and  refined  from  microwave  radar  systems  to  accom¬ 
modate  the  special  needs  of  electro-optical  systems  (SO).  The  detailed  capabilities 
of  current  LADAR  technology  is  the  subject  of  further  analysis  in  Chapter  III,  with 
emphasis  on  specific  deep-space  requirements. 

Summary 

The  current  deep-space  surveillance  network  is  not  achieving  the  track  rates 
necessary  to  accomplish  all  operational  objectives.  In  particular,  although  GEODSS 
detection  performance  is  sufficient  for  deep-space  surveillance,  availability  is  severely 
hampered  by  adverse  weather  and  bright-sky  conditions.  An  attractive  alternative 
to  the  GEODSS  method  of  passive  sensing  in  the  visible  region  uses  an  IR  laser 
to  irradiate  the  target.  Although  active  IR  methods  have  been  implemented  in  the 
Firepond  LADAR,  capability  is  limited  to  detection  of  unenhanced  satellites  at  low 
altitude.  Extension  of  this  low-altitude  capability  to  deep-space  appears  promising 
with  the  benefit  of  recent  developments  in  IR  laser,  detector,  and  optics  technology. 


III.  LADAR  Techniques 


Introduction 

This  chapter  describes  various  LADAR  techniques  suitable  for  a  long-range 
LADAR,  and  reviews  the  currently  available  technology  for  major  LADAR  elements. 
In  particular,  the  typical  characteristics  of  deep-space  targets  are  discussed,  followed 
by  an  outline  of  detection  methods  suitable  for  these  targets.  Current  technology  is 
then  reviewed  for  laser  transmitters,  detectors,  and  optics.  Finally,  a  target  acquisi¬ 
tion  scheme  is  proposed  to  optimize  the  target  search  phase  of  LADAR  detection. 

Target  Characteristics 

The  distribution  of  deep-space  objects  which  are  of  primary  interest  to  this 
research  can  be  categorized  according  to  three  orbital  classes:  geosynchronous,  semi- 
synchronous,  and  highly-elliptical.  Figure  1  provides  a  pictorial  representation  of 
each  orbital  class.  Successful  detection  of  targets  from  all  three  classes  depends  on, 
among  other  factors,  the  relative  viewing  geometry  between  the  ground-based  sensor 
and  the  target  object.  For  example,  only  a  fraction  of  the  geosynchronous  population 
are  'visible’  from  a  fixed  sensor  site  all  of  the  time,  while  the  remaining  fraction  may 
never  be  visible  from  the  same  site  due  to  obscuration  by  the  earth.  Conversely, 
objects  with  much  shorter  orbital  periods  (lower  altitudes)  may  be  visible  several 
times  daily  for  short  durations  from  a  single  site,  depending  on  the  exact  orbital 
characteristics. 

Clearly,  a  global  network  of  sensors  is  required  to  provide  adequate  coverage  of 
the  majority  of  deep-space  objects.  The  target  range  from  any  one  of  these  sensors 
is  then  dependent  on  the  relative  position  of  the  object  along  its  orbit.  Although 
targets  in  the  highly-elliptical  orbit  class  can  occupy  altitudes  as  high  as  40,000  km  or 
more,  successful  orbit  predictions  can  be  performed  from  observations  taken  closer  to 


13 


the  target’s  perigee,  which  is  only  a  few  hundred  kilometers  altitude.  Assuming  that 
observations  for  this  class  are  not  essential  near  apogee,  the  maximum  target  range 
experienced  by  a  deep-space  sensor  would  then  occur  when  observing  geosynchronous 
objects. 

For  an  equatorial  sensor  which  is  situated  coincident  with  the  orbital  sub- 
point  of  a  geostationary  satellite  (zero  degree  inclination),  the  target  range  is  nearly 
36,000  km.  In  general,  however,  the  satellite-sensor  geometry  is  more  complicated, 
as  shown  in  Figure  2.  Based  on  GEODSS  site  characteristics,  the  sensor  may  be 
located  at  latitudes  of  up  to  36  degrees  (30:111-22),  while  the  satellite  orbit  may  be 
inclined  by  several  degrees  to  the  equatorial  plane.  Since  the  target  range  is  also 
a  function  of  the  satellite  sub-point,  a  wide  range  of  satellite-sensor  geometries  is 
possible. 
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Figure  2.  Typical  Sensor-Target  Viewing  Geometry 

The  limits  of  sensor  coverage  can  be  defined  in  terms  of  maximum  look  angle . 
which  is  the  average  maximum  sensor  viewing  angle  for  successful  target  detection, 
measured  from  the  sensor  local  zenith.  Although  actual  values  of  maximum  look 
angle  for  successful  detection  of  a  specific  target  will  vary  due  to  atmospheric  effects, 
background  radiance,  and  target  characteristics,  an  average  maximum  sensor  look 
angle  is  appropriate  for  design  calculations.  For  example,  GEODSS  maximum  look 
angles  are  generally  less  than  60  degrees  (30:111-22).  Th;s  value  will  be  adopted  as 
a  minimum  design  objective  to  ensure  adequate  sensor  coverage.  This  equates  to 
target  ranges  for  geosynchronous  objects  of  between  36,000  and  40,000  km. 

Apart  from  target  range,  the  other  major  target  characteristic  that  affects 
detection  is  the  target’s  LADAR  cross-section ,  which  is  a  function  of  target  size 
and  orientation,  surface  shape  and  roughness,  and  surface  reflectivity.  For  LADAR 
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detection  purposes,  target  surfaces  may  be  classified  as  either  specular  o r  diffuse.  A 
specular  target  is  composed  of  surfaces  which  are  ‘flat’  at  IR  wavelengths,  and  reflect 
incident  radiation  preferentially  in  a  specific  direction.  A  diffuse,  or  rough,  target  is 


composed  of  many  independent  scatterers  and  reflects  incident  radiation  equally  in 
all  directions. 

In  general,  objects  such  as  satellites  may  be  modeled  as  hybrid  targets,  ex¬ 
hibiting  1R  reflectivity  characteristics  which  include  both  diffuse  and  specular  com¬ 
ponents.  In  terms  of  LADAR  cross-section,  the  addition  of  a  relatively  small  diffuse 
reflectivity  component  to  a  purely  specular  target  surface  causes  the  overall  reflected 
signal  to  resemble  that  from  a  purely  diffuse  target  (6:43).  Therefore,  without  loss  of 
generality,  this  analysis  assumes  that  typical  deep-space  targets  are  characterized  by 
surfaces  with  a  purely  diffuse  IR  reflectivity.  Using  a  diffuse  spherical  target  model, 
with  diameter  much  greater  than  the  incident  wavelength,  the  LADAR  cross-section, 
eric,  may  consequently  be  written  as  follows  (48:20 IB): 


07m 


(1) 


where  pXT  is  the  IR  reflectivity  of  the  target  surface,  and  A  is  the  projected  physical 
target  area. 

Typical  target  cross-sections  can  be  determined  from  typical  values  of  A  and 
ptr.  The  size  and  shape  of  deep-space  objects  varies  considerably,  with  target  char¬ 
acteristic  dimensions  ranging  from  1  to  30  meters  (84:37).  For  detection  analysis 
purposes,  a  spherical  target  with  a  diameter  of  3  meters  is  therefore  assumed  as 
typical. 

The  IR  reflectivity  of  deep-space  objects  also  varies,  depending  on  the  com¬ 
position  of  target  surface  material.  For  satellites,  surface  materials  include  solar 
cell  components,  optical  solar  reflectors,  multi-layer  insulation,  surface  paints,  and 
heat  radiators  (2).  The  majority  of  such  materials  are  designed  for  satellite  thermal- 
management  purposes,  and  are  characterized  by  high  reflectivity  in  the  optical  region 
and  high  emissivity  in  the  IR  region,  so  that  heat  absorption  and  retention  is  mini¬ 
mized.  Unfortunately,  the  goals  of  satellite  thermal  management  design  oppose  the 
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desire  for  LADAR  targets  with  high  IR  reflectivity.  Accordingly,  typical  IR  reflec¬ 
tivity  values  for  thermal  control  surfaces  are  only  about  20%  (3:9). 

Substituting  .4  =  ~(1.5)2  and  pir  =  0.2  into  Equation  (1)  yields  a  value  of 
<T;CJ  %  4  rrr  .  which  is  suitable  for  further  analysis  of  design  feasibility.  With  typical 
examples  of  target  orbits,  maximum  range,  and  LADAR  cross-section  thus  defined, 
methods  of  long-range,  high-resolution,  LADAR  detection  will  now  be  discussed. 

LADAR  Detection  Methods 

LADA.R  systems  may  be  categorized  according  to  receiver  design  as  either 
direct  detection,  or  coherent  detection  systems.  In  the  direct  detection  scheme,  de¬ 
tection  decisions  are  based  on  the  results  of  a  direct  measurement  of  received  signal 
energy,  losing  a  pulsed  transmitter  and  measuring  the  pulse  round-trip  time  gives 
the  target  range,  while  the  addition  of  shaft-position  encoders  on  the  receiver  tele¬ 
scope  mount  allows  calculation  of  the  target's  angular  position.  The  major  benefits 
of  direct  detection  are  simplicity  of  hardware  design  and  relaxed  requirements  for 
stringpnt  laser  frequency  stability  Hmvevcr,  high  levels  of  background  noise  severely 
degrade  the  performance  of  direct  detection  systems. 

In  coherent  detection  schemes,  the  received  signal  is  mixed  with  a  local  os¬ 
cillator  (LO)  signal  on  the  photodetector  surface.  Photocurrent  is  produced  in  the 
detector  at  a  beat  frequency  equal  to  the  difference  between  signal  and  local  oscillator 
frequencies.  Normally,  the  local  oscillator  is  offset  from  the  transmitter  frequency  by 
an  intermediate  frequency  (IF),  so  that  the  photocurrent  beat  frequency  equals  the 
IF  frequency  plus  any  Doppler  shift  imposed  by  target  radial  velocity.  This  coherent 
detection  scheme  is  known  as  heterodyne  detection  and  is  illustrated  for  a  pulsed 
system  in  the  LADAR  block  diagram  shown  in  Figure  3.  If  a  single  laser  is  used 
to  provide  both  transmitter  and  local  oscillator  signals,  the  system  is  referred  to  as 
an  offset  homodyne  scheme  (49:153);  if  no  offset  is  used,  the  system  is  simply  called 
. omodyne  detection  (24:318). 


LEGEND 

ft  =  Transmitter  Frequency  Az/El  =  Telescope  Azimuth  and  Elevation 

/./  =  Intermediate  Frequency  BS  =  Beam-Splitter 

Id  =  Doppler-Shift  Frequency 


Figure  3.  Block  Diagram  of  Pulse-Doppler  LADAR  using  Heterodyne  Detection 


8 


Figure  3  a  Is:  chows  the  use  of  range  gates  for  signal  timing.  Range  gates 
are  used  to  control  receiver-on  time  so  that  detector  output  is  only  enabled  during 
the  period  in  which  a  transmitted  pulse  is  expected  to  return.  For  example,  the 
pulse  round-trip  time  to  a  target  at  36,000  km  range  is  about  240  milliseconds. 
Allowing  for  atmospheric  propagation  delays,  optics  delays,  variation  in  the  target 
range,  and  a  typical  pulse  duration  of  1  /.is,  the  receiver  might  be  triggered  'on'  after 
a  delay  of  239  ms  from  pulse  transmission,  and  then  triggered  'off'  after  a  further 
2  ms  delay.  This  method  prevents  erroneous  detections  during  the  239  ms  between 
pulse  transmission  and  the  earliest  expected  pulse  return,  and  allows  more  efficient 
searching  of  pre-defined  range  cells.  The  improvement  in  detection  performance 
through  the  use  of  range  gates  will  be  further  discussed  in  Chapter  V. 

One  major  element  normally  found  in  the  post-detection  circuitry  shown  in 
Figure  3  is  the  matched  filter.  The  matched  filter  is  also  triggered  by  range  gate 
signals  and  is  the  optimum  filter  when  the  signal  is  corrupted  by  additive  white 
Gaussian  noise  (73:3-8).  The  ideal  matched  filter  is  usually  approximated  as  a  filter 
with  a  frequency  response  matched  to  the  spectrum  of  the  transmitted  pulse  (88:72). 
A  bank  of  such  filters,  with  bandpasses  which  are  contiguous  in  the  frequency  domain, 
provides  ease  of  Doppler-shift  measurement  and  closer  matching  of  the  filter  to  the 
received  Doppler-shifted  signal. 

The  main  benefits  of  coherent  detection  schemes  are  additional  target  informa¬ 
tion  through  measurement  of  Doppler  shift,  and  the  capability  of  obtaining  quantum- 
limited  detection  with  high  background  noise  (73:37-4).  Heterodyne  systems,  there¬ 
fore,  may  be  “several  orders  of  magnitude  more  sensitive  than  a  direct  detection 
system  which  is  typically  background  noise  or  amplifier  noise  limited'’  (27:93).  For 
a  typical  10.6  /m i  LADAR,  heterodyne  reception  should  allow  detection  of  signals 
which  are  100  to  1000  times  weaker  than  the  equivalent  direct  detection  receiver 
(17:290).  Further  sensitivity  improvements  are  possible  by  transmitting  a  train  of 
pulses  and  then  integrating  the  received  signal  over  several  pulses. 
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Pulse  Compi'ession.  Improvements  in  range  resolution  can  be  gained  by  using 
short  pulse  widths;  however,  many  high-power  lasers  are  characterized  by  long  pulses 
(10-1000  i is )  of  high  energy  (50-300  J).  A  common  technique  used  to  achieve  high 
range-resolution  with  long  pulses  is  pulse  compression  (72:121).  The  simplest  pulse 
compression  technique  uses  frequency  modulation  (FM)  to  generate  a  transmitter 
pulse  of  constant  amplitude  with  linearly  increasing  frequency  throughout  the  pulse, 
also  known  as  chirp.  The  matched  filter  of  Figure  3  then  acts  as  a  pulse  compression 
filter  by  adding  a  signal  time  delay  which  is  inversely  proportional  to  frequency.  The 
added  time  delay  allows  the  higher  frequencies  in  the  pulse's  trailing  edge  to  'catch 
up’  to  the  leading  edge  so  that  the  wide,  low-amplitude  received  pulse  is  compressed 
into  a  short,  high-amplitude  pulse  (73:422). 

A  typical  implementation  of  the  pulse  compression  filter  uses  a  surface  acoustic- 
wave  device  made  of  a  piezoelectric  material  such  as  lithium  niobate  (36:439).  If  the 
FM  bandwidth  in  the  transmitted  pulse  is  Bp\t ,  with  pulse  width  r,  then  the  ratio 
of  uncompressed  pulse  width  to  compressed  pulse  width  —  or  compression  ratio  — 
is  r  B/r,\f.  The  peak  power  of  the  received  pulse  is  also  increased  by  the  compression 
ratio  r  Bp\i  after  passing  through  the  filter  (72:423). 

The  pulse  compression  method  may  also  be  used  with  the  post-detection  inte¬ 
gration  scheme.  This  allows  high  resolution  to  be  combined  with  the  high  sensiti',i(-y 
possible  through  using  coherent  detection.  Although  several  advantages  are  gained 
by  using  coherent  detection  schemes,  some  disadvantages  exist.  The  primary  draw¬ 
backs  are  critical  optical  alignment  tolerances,  hardware  complexity,  and  stringent 
stability  requirements  for  the  transmitter  and  local  oscillator.  Nevertheless,  for  the 
extremely  long  ranges  under  consideration,  high  receiver  sensitivity  is  of  prime  im¬ 
portance  in  LADAR  design.  For  this  reason,  and  because  the  target’s  radial  velocity 
component  is  available  from  the  measured  Doppler  shift,  the  LADAR  system  ot  Fig¬ 
ure  3  will  be  adopted  as  the  basis  for  further  system  analysis.  The  performance 
requirements  of  each  LADAR  system  element  will  now  be  examined  in  detail. 
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LADAR  Transmitters 


The  characteristics  required  of  a  laser  which  is  suitable  for  long-range  coherent 
IR  LADAR  applications  include  single- frequency  operation,  high  pulse  energy,  high 
frequency  stability  (both  inter-pulse  and  intra-pulse),  and  high  pulse  repetition  rate. 
Several  transmitter  configurations  are  currently  in  use.  with  the  majority  utilizing 
Transversely-Excited  Atmospheric-pressure  carbon  dioxide  (TEA  C0->)  lasers  as  the 
radiation  source,  at  a  wavelength  of  10.6  /zm  (88:83;  39:2).  Two  other  important 
sources  exist  for  LADAR  transmitters:  Neodymium-doped  Yttrium  Aluminum  Gar¬ 
net  (Nd:YAG)  at  L .06  //m  (27:96).  and  Deuterium  Fluoride  (DF)  at  about  3.8 
(9:8). 

Due  to  moderate  atmospheric  absorption  at  10.6  /zm.  two  variants  of  the  ordi¬ 
nary  CO2  laser  have  been  recently  considered.  The  first  uses  isotopicallv-substituted 
l3Clb02  in  lieu  of  12C1602  in  the  lasing  gas  mixture,  which  produces  radiation  at 
about  11.1  /zm  (57:  31:  27).  Since  carbon-13  abundance  in  the  atmosphere  repre¬ 
sents  only  1.1%  of  the  total  C02  concentration  (87:5-61),  self-absorption  is  greatly 
reduced  at  11.1  /zm.  Performance  from  13C1602  lasers  is  reported  to  equal  that  of 
l2Cls02  lasers  (57:106). 

The  second  C02  variant  uses  a  frequency-doubling  crystal  such  as  Tl3AsSe3 
or  AgGaSe2  to  convert  C02  9-11  /<m  to  '.".-o  "  radiation  (56;  27). 

Conversion  efficiencies  of  40%,  or  higher,  for  frequency-doubled  C02  are  advocated 
(27:99;  56:6). 

Although  several  other  laser  sources  are  currently  under  consideration  for 
LADAR  use,  the  C02,  Nd:YAG,  and  DF  lasers  represent  the  most  technologically- 
mature  examples  of  high-power  pulsed  IR  LADAR  transmitters  currently  available. 
The  characteristics  of  each  source  will  now  be  described  individually. 

CO 2  Transmitters.  C02  LADAR  technology  is  the  oldest  and  most  mature  of 
the  sources  under  consideration  (88:83).  High  electrical  efficiencies  of  around  20%  are 
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achievable,  which  significantly  reduces  power  generation  and  handling  requirements 
when  compared  to  other  systems.  The  output  pulse  from  a  TEA  CO;  laser  usually 
consists  of  a  short,  high-amplitude  spike  of  30  to  200  ns  duration,  followed  by  a 
long,  low-amplitude  'tail'  of  1  to  5  fis  duration  (52:74).  Since  ideal  LADAR  pulse 
characteristics  include  high  energy  and  short  duration,  the  presence  of  the  long  tail 
in  the  TEA  CO:  pulse  is  undesirable.  Recent  experiments  have  been  performed 
successfully  with  intra-cavitv  SF^-He  saturable-absorber  cells  to  greatly  reduce  the 
tail  duration  without  loss  of  peak  pulse  power  (52:77). 

For  high-power  applications,  several  different  laser  arrangements  are  possible 
such  as  hybrid  TEA,  injection-locked  TEA.  and  Master  Oscillator-Power  Amplifier 
(MOPA)  configurations.  The  hybrid  TEA  set-up  consists  of  a  TEA  gain  section  and 
a  low  pressure  continuous- wave  (CW)  section  in  a  common  optical  cavity  (17:290: 
60:149).  Although  good  inter-pulse  frequency  stability  is  achieved  with  this  system, 
output  pulse  energy  and  pulse  repetition  rates  are  limited  due  to  saturation  of  the 
CW  section,  and  the  onset  of  heating  and  acoustic  effects  (60:149-150). 

The  injection-locked  technique  uses  an  external  low-power  CW  oscillator  as  an 
injection  source  into  a  TEA  cavity.  In  this  case,  laser  output  frequency  is  controlled 
by  the  TEA  cavity.  Mode  selection  in  the  TEA  section  only  requires  a  small  injected 
signal  when  the  difference  between  TEA  tuned  frequency  auu  injection  frequency  is 
small  (61:11).  Therefore,  careful  frequency  control  between  the  TEA  section  and  the 
injected  signal  is  required. 

Additional  problems  exist  with  injection-locked  systems  in  the  form  of  intra¬ 
pulse  frequency  instability,  or  ‘chirp’.  Frequency  chirp  results  primarily  from  a  phe¬ 
nomenon  known  as  laser-induced  medium  perturbation  (LIMP)  (32:56;  66:62).  The 
origin  of  LIMP  lies  in  "‘differential  heating  between  lasing  and  non-lasing  levels 
(32:56),  which  gives  rise  t.o  density  variations  in  the  lasing  medium:  the  resultant 
change  in  medium  refractive  index  causes  intra-pulse  chirp  on  the  order  of  a  few 
hundred  kilohertz  (32:56). 
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The  MO  PA  configuration  consists  of  a  low-pressure  master  oscillator  operat¬ 
ing  at  low  pulse  energies,  which  is  then  directed  iruo  a  multi-pass  FkA  amphtier 
configured  as  an  unstable  resonator  (60:151:  32:57).  MOPA  configuration-  have  -e\  - 
eral  advantages  over  the  previous  systems  and  are  characterized  by  low  level-  of 
frequency  chirp,  high  gain,  and  improved  control  of  pulse  width  and  shape  1 32:57: 
66:63).  For  these  reasons,  MOPA  systems  are  the  preferred  configuration  for  moder¬ 
ate  to  high-power  systems,  exhibiting  pulse  repetition  rates  of  several  hundred  pulses 
per  second,  with  typical  pulse  energies  of  2  to  3  joules,  and  pulse  widths  of  1  to  5  ys 
(66:63:  32:59:  60:151). 

Higher  pulse  energies  with  longer  pulse  widths  can  be  obtained  from  the  ampli¬ 
fier  by  using  different  laser  excitation  techniques.  Using  an  electron-beam  sustained 
discharge  in  a  13C02  amplifier,  for  example,  Moody  et  al.  report  a  pulse  energy  of 
200  J  and  pulse  width  of  30  /.is  (57:102). 

Nd:YAG  Transmitters.  Pulsed  N’d:YAG  Users  of  moderate  to  high-power  have 
been  used  extensively  in  a  variety  of  applications  for  many  years.  I.  nfortunately, 
the  “poor  transverse-mode  quality  and  broad  spectral  width  typical  of  NdAAG 
lasers"  (27:97)  renders  them  poor  candidates  for  high-resolution  LADAR  systems. 
Additionally,  atmospheric  scintillation  is  most  severe  at  shorter  wavelengths  (>S:83L 
reducing  the  probability  of  reliable  LADAR  operation  at  1.06  /mi. 

DF  Transmitters.  DF  lasers  have  undergone  considerable  development  in  re¬ 
cent  years,  and  have  been  implemented  successfully  in  ranging  applications,  with 
peak  powers  above  2  MW  and  pulse  lengths  less  than  10C  ns  (9:8).  However,  DP 
lasers  normally  operate  in  multi-line  mode,  with  power  divided  over  about  32  dis¬ 
crete  lines  (42:70).  Single-line  operation  is  possible,  but  at  conversion  efficiencies 
of  around  only  20%  of  multi-line  power  (67:13).  Since  reported  repetition  rates  are 
also  low  at  around  60  Hz  (9:11),  DF  lasers  require  further  scaling  improvements  to 
warrant  serious  consideration  for  high-power  LADAR  application. 
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Detectors 


As  LADAR  transmitter  technology  has  evolved  and  improved,  great  advances 
in  IR  detector  technology  have  also  been  made.  The  single  detector  elements  used  in 
the  late  1950s  were  replaced  by  linear  detector  arrays  in  the  1960s  and  early  1970s. 
with  two-dimensional  focal  plane  arrays  (FPAs)  becoming  state-of-the-art  in  the  late 
1970s  (16:60).  Since  then.  FPA  development  has  concentrated  on  improved  manu¬ 
facturing  techniques  and  array  uniformity,  and  fabrication  of  higher  array  densities. 
For  LADAR  applications,  higher  spatial  resolution  requires  a  higher  density  of  FPA 
elements,  while  reliable  detection  requires  a  highly  uniform  detector  responsivity 
across  all  FPA  elements. 

If  linear  arrays  are  used,  some  form  of  mechanical  scanning  is  usually  required 
to  cover  the  desired  field-of-view.  However,  the  use  of  two-dimensional  FPAs  allows 
the  array  to  'stare'  at  a  scene  for  a  fixed  time  before  moving,  or  'stepping',  to  the 
next  scene.  The  'step-staring1  method,  which  is  adopted  for  this  analysis,  greatly 
simplifies  LADAR  design  since  precise  mechanical  scanning  of  the  field-of-view  is  no 
longer  required,  while  sensitivity  is  increased  from  longer  signal  integration  times 
(16:61). 

An  additional  detector  requirement  for  use  with  high-speed  LADAR  targets  is 
high  bandwidth.  For  the  wavelengths  and  radial  velocity  components  of  the  targets 
under  consideration,  the  detector  must  be  able  to  respond  to  a  maximum  IF  band¬ 
width  of  about  2  GHz.  Fortunately.  FPAs  are  routinely  fabricated  with  baudwidths 
which  meet  this  requirement  (64:53). 

The  majority  of  detector  advances  have  so  far  been  concentrated  on  the  wave¬ 
length  regions  corresponding  to  the  two  major  IR  atmospheric  transmission  windows 
at  3-5  n m  and  8-12  fim.  However,  a  significant  decrease  in  available  detector  per¬ 
formance  exists  when  moving  from  short  to  long-wavelength  detector  operation.  For 
example,  state-of-the-art  FPA  technology  for  visible  wavelengths  consists  of  an  array 
of  2048  x  2048  elements;  at  3-5  /im,  current  array  limits  are  512  x  512  elements. 
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while  the  8-12  /im  region  is  currently  limited  to  12S  x  128  elements  (22:1-5).  Detec¬ 
tor  element  sizes  for  these  arrays  vary  from  10  /z m  x  10  //rn  to  50  /z m  x  50  /zm. 
depending  on  the  material  used  (22:13;  21:74). 

Various  FPA  materials,  integration  techniques,  signal  multiplexing,  and  data 
processing  methods  are  in  use  or  are  being  developed  for  IR  detectors  The  primary 
configurations  currently  in  use  for  the  3-5  /zm  band  are  made  from  Hgi_xCdrTe, 
In  Sb,  Pt  Si,  or  PbS,  while  for  the  8-12  /z m  band,  only  Hgi_xCdxTe  has  been  de¬ 
veloped  to  a  sufficient  level  (11:114).  The  value  of  x  in  the  Hg]_rCdrTe  material 
is  varied  to  achieve  a  desirable  material  bandgap,  detector  cut-off  wavelength,  and 
peak  responsivity.  For  example,  a  value  of  x  —  0.205  for  a  detector  cooled  to  77  °K 
provides  a  material  with  quantum  efficiency  exceeding  S5%  (with  anti-reflective  coat¬ 
ing),  cut-off  wavelength  of  11.7  /z m,  and  peak  responsivity  at  10.6  /zm  (64:52). 

Two-dimensional  FPAs  are  made  using  both  Charge-Coupled-Device  (CCD) 
and  Field-Effect-Transistor  (FET)  switches  to  address  the  detectors,  with  both  types 
allowing  full  two-dimensional  multiplexing  (11:114;  22:15).  However,  two  practical 
limitations  exist  for  staring  FPAs  which  arise  from  finite  CCD  storage  capacity  and 
data  readout  time.  These  limitations  are  reflected  in  a  specified  'maximum  FPA 
stare  time'  and  ‘maximum  data  output  rate’  for  FPAs.  Typical  values  for  maximum 
stare  times  are  200-2000  /z s  for  the  3-5  /zm  band,  and  5-50  /zs  for  the  8-12  /zm 
band,  while  typical  maximum  output  rates  are  reported  to  be  about  10  MHz  for  the 
3-5  /zm  band  and  2  MHz  for  the  8-12  /zm  band  (11:115). 

Optics 

The  purpose  of  LADAR  optics  is  to  efficiently  deposit  the  maximum  amount 
of  transmitter  energy  onto  the  target,  and  then  to  efficiently  focus  the  maximum 
°  amount  of  reflected  signal  from  the  target  onto  the  detector.  The  tasks  of  transmit¬ 
ting  and  receiving  energy  in  this  manner  need  not  be  performed  by  the  same  set  of 
optics,  but  is  often  performed  with  one  transmit/receive  telescope  due  to  reduced 
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cost  and  complexity.  On  the  other  hand,  coherent  LADAR  performance  is  signifi¬ 
cantly  improved  in  terms  of  signal-to-noise  ratio  by  using  a  receiver  aperture  which 
is  up  to  three  times  larger  than  the  transmitter  aperture  (83:117). 

The  aperture  used  in  IR  LADAR  can  be  considered  to  be  the  IR  equivalent 
of  an  antenna  at  microwave  frequencies.  In  the  IR  case,  several  reflectors  and/or 
lenses  are  used  in  a  telescope  configuration,  with  the  detector  or  FPA  situated  at 
the  telescope's  effective  focal  length.  Additional  beamsplitters,  mirrors,  and  lenses 
are  used  to  direct  laser  output  radiation  out  through  the  telescope  in  much  the  same 
manner  as  microwave  feed-horns.  If  a  single  telescope  is  used  for  both  the  transmit 
and  receive  functions,  a  system  of  mirrors  is  used  as  a  ‘transmit/receive'  switch  to 
alternately  direct  laser  energy  out  through  the  telescope  during  the  transmit  cycle, 
and  direct  received  radiation  onto  the  detector  during  the  receive  cycle. 

Telescope  configurations  vary  widely;  however,  analyses  have  shown  that  the 
best  performance  is  obtained  from  the  Cassegrain  configuration  (34;  85).  In  partic¬ 
ular,  the  aplanatic  version  of  the  Cassegrain  design  —  the  Ritchey  Chretien  con¬ 
figuration  —  is  favored  for  its  compactness,  superior  alignment,  and  low  aberration 
characteristics  (34:161).  An  example  of  such  a  design  is  illustrated  in  Figure  4. 
showing  a  typical  Cassegrainian  arrangement.  Typical  Ritchey-Chretien  dimensions 
are  shown  in  Table  2  in  Chapter  II  for  the  GEODSS  main  telescope.  Diameters 
of  secondary  mirrors  for  such  telescopes  range  from  20%  to  35%  of  primary  mirror 
diameters  (44:979),  which  results  in  partial  obscuration  of  the  signal.  This  effect  will 
be  examined  in  detail  in  Chapter  V. 

Primary  mirror  diameter  is  the  limiting  factor  in  manufacturing  large  tele¬ 
scopes  for  high-resolution  applications.  The  fabrication  of  mirrors  larger  than  2  to 
3  meters  is  best  performed  by  first  making  several  smaller  mirror  segments,  and  then 
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assembling  them  into  a  single  element  afterwards  (38:58).  A  prime  example  of  this 
technology  is  the  Hubble  Space  Telescope  mirror,  which  has  a  diameter  of  2.4  meters 
and  was  made  as  a  welded  assembly  of  ultra-low  expansion  fused-silica  (38:58). 


Figure  4.  Cassegrainian  Telescope  Configuration 

One  other  limitation  encountered  when  designing  LADAR  optics  is  possible 
damage  of  optical  surfaces  when  using  high-power  transmitters.  For  this  reason, 
care  must  be  taken  to  avoid  maintaining  high-power  laser  output  a s  a  small  diameter 
beam  when  passing  through  lenses  or  reflecting  off  mirrors.  For  example,  experiments 
have  shown  that  surface  reflectivity  of  copper  and  molybdenum  mirrors  is  altered 
when  CO2  laser  irradiance  exceeds  1.2  GW/cm2  (79:415). 

Acquisition,  Tracking,  and  Pointing 

Acquisition,  tracking,  and  pointing  (ATP)  functions  of  a  LADAR  system  con¬ 
sist  of  three  related  tasks.  Acquisition  involves  searching  of  a  specific  volume  in 
target  space  with  the  aim  of  identifying  a  suitable  target.  Pointing  involves  the 
accurate  positioning  of  the  telescope  boresight,  or  instrumental  line  of  sight  (LOS), 
in  the  same  direction  as  the  target  LOS.  Therefore,  pointing  error  can  be  defined  as 
the  average  difference  between  instrumental  LOS  and  target  LOS  (80:40).  Using  the 
same  terminology,  tracking  involves  minimizing  any  variation  in  pointing  error  over 
time. 

The  LADAR’s  ability  to  perform  effective  pointing  and  tracking  is  largely  de¬ 
termined  by  the  quality  of  electro-mechanical  control  systems  used  to  position  the 
telescope  and  control  the  various  beam-steering  mirrors.  Several  design  concepts 


exist  to  perform  these  functions,  including  inertial  stabilization,  mirror  stabiliza¬ 
tion,  gear-driven  gimbals,  momentum  wheels,  shaft-angle  encoders,  active  structural 
elements,  and  thermal  control  sub-systems.  Detailed  discussion  of  these  topics  is 
beyond  the  scope  of  this  research  and  is  comprehensively  treated  by  Ulich.  Masten 
and  Hilkert,  and  Stockum  et  al.  (SO;  51:  75).  Additionally,  a  detailed  analysis  of  dis¬ 
turbances  such  as  wind,  seismic,  gyroscope  noise,  electronics  noise,  and  atmospheric 
turbulence  is  performed  by  Merritt  and  Meline  (55). 

While  tracking  and  pointing  performance  is  primarily  an  electro-mechanical 
design  problem,  selection  of  an  acquisition  scheme  is  a  fundamental  step  which  affects 
the  choice  of  detector,  optics,  and  transmitted  waveform  characteristics.  In  the 
case  under  study,  the  choice  of  a  staring  FPA  and  pulsed  waveform  to  obtain  high 
resolution  and  sensitivity  places  limits  on  the  number  of  suitable  acquisition  methods. 
Also,  the  long  round-trip  time  of  at  least  0.24  seconds  to  geosynchronous  orbit  for 
transmitted  pulses  implies  a  relatively  long  search  time  for  each  stage  of  the  step- 
stare  sequence,  assuming  only  a  single  telescope  is  used. 

However,  search  time  could  be  reduced  if  two  telescopes  are  fixed  to  the  same 
mount,  with  independent  LOS  control  systems.  The  first  'master'  telescope  would 
perform  as  transmitting  antenna,  and  be  'pointed-ahead'  along  the  search  path  by 
an  angular  distance  of  #iead-  The  other  ‘slave’  telescope,  configured  as  receiving 
antenna,  would  always  lag  the  transmitter  telescope  position  by  9^^  throughout  the 
step-stare  acquisition  phase.  In  this  manner,  the  ‘stare’  process  is  being  performed 
by  the  receiver  telescope  while  several  pulse  trains  are  ‘in-flight’. 

The  determination  of  9^^  can  be  written  as  a  multiple  of  the  receiver  telescope 
field-of-view,  which  is  equivalent  to  a  multiple  of  the  number  of  steps  in  the  step- 
stare  sequence.  That  is,  if  the  time  of  flight  for  one  pulse  is  expressed  as  2 R/c,  where 
R  is  the"  target  range,  and  c  is  the  speed  of  light,  and  the  total  waveform  duration  is 
denoted  by  T,  then 

2  R 

=  — =  «r  (2) 

c  T 
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where  receiver  processing  and  telescope  stepping  delays  have  been  ignored,  and  aT 
is  the  receiver  angular  field-of-view. 

To  alleviate  interference  problems  associated  with  backscatter  of  in-flight  wave¬ 
forms,  variations  in  pulse  transmission  rate  can  be  implemented  in  a  similar  way  to 
PRF  switching  or  jitter  in  microwave  applications.  With  this  refinement,  each  trans¬ 
mitted  waveform  is  given  a  unique  pulse  spacing  which  is  synchronized  with  receiver 
range-gate  circuitry.  Therefore,  the  receiver  only  locks-on  to  the  expected  synchro¬ 
nized  return  waveform. 

This  ‘master-slave’  acquisition  method  has  the  potential  to  dramatically  re¬ 
duce  search  times  when  compared  with  a  single-telescope  design.  The  actual  search 
time  is  limited  primarily  by  the  stare  period,  which  depends  on  the  duration  of  the 
transmitted  waveform.  An  additional  benefit  gained  from  using  this  method  is  the 
potential  to  gain  performance  improvements  through  selection  of  a  receiver  aperture 
which  is  larger  than  the  transmitter  aperture. 

Summary 

This  chapter  has  examined  the  various  LADAR  techniques  available  for  design 
of  a  ground-based,  deep-space  LADAR.  Deep-space  target  characteristics  have  been 
examined,  with  selection  of  the  following  target  specifications  for  further  analysis  of 
design  feasibility: 


Maximum  Target  Range  =  40,000  km 
Maximum  Look  Angle  >  60  degrees 
LADAR  Cross  Section  of  Target  =  4  m2 

Various  detection  methods  have  been  reviewed,  with  the  heterodyne  detection 
method  selected  as  the  optimum  LADAR  choice  due  to  superior  receiver  sensitivity, 
and  the  availability  of  Doppler  shift  information.  Further  sensitivity  improvements 
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are  possible  by  performing  post-detection  integration  over  several  transmitted  pulses, 
while  range  resolution  and  received  pulse  power  can  be  dramatically  improved  using 
pulse  compression  techniques.  A  general  block  diagram  of  the  proposed  LADAR 
system  is  shown  in  Figure  3. 

High  power  lasers  suitable  for  iong-range  IR  LADAR  were  studied,  with  the 
following  types  highlighted  as  technologically-mature,  high-power  LADAR  transmit¬ 
ters: 


Nd:YAG  S'  1.06  n m 

DF  3.8  /im 

Frequency-doubled  12C02  @  4. 5-5. 5  /im 
12C02  @  10.6  fj. m 
13C02  @  11.1  /im 

While  nC02  lasers  were  found  to  provide  similar  performance  to  12C02  , 
frequency-doubled  versions  of  12C02  are  also  available  at  conversion  efficiencies  of 
about  40%.  DF  lasers,  on  the  other  hand,  are  capable  of  single-line  operation  at 
conversion  efficiencies  of  only  20%,  while  Nd:YAG  lasers  suffer  from  poor  quality 
output  characteristics.  Therefore,  lasers  based  on  C02  appear  to  be  the  best  choices 
for  LADAR  transmitters  due  to  their  degree  of  technological  maturity,  high  power, 
narrow  linewidth,  and  small  chirp  characteristics.  Of  the  possible  C02  laser  configu¬ 
rations,  the  MOPA  configuration  was  selected  for  its  low  frequency-chirp,  high  gain, 
and  superior  control  of  pulse  characteristics. 

High  spatial  resolution  was  then  used  as  a  measure  of  effectiveness  in  reviewing 
detector  technology,  with  two-dimensional  FPAs  providing  the  best  performance. 
Array  sizes  for  current-technology  IR  FPAs  are  512  x  512  elements  for  the  3-5  /im 
band,  and  128  x  128  elements  for  the  8-12  /im  band.  Typical  square  element  sizes 
range  from  10  /im  to  50  /im. 
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The  major  design  variable  in  choosing  optics  is  the  telescope  configuration. 
The  Ritchey-Chretien  variant  of  the  Cassegrainian  arrangement  was  found  to  be 
optimum  in  terms  of  size,  alignment,  and  aberration  qualities.  Relatively  large  pri¬ 
mary  reflectors  can  now  be  accommodated  in  the  manufacture  of  telescope  optics  by 
welding  individual  mirror  segments  together.  However,  care  must  be  taken  during 
LADAR  operation  to  avoid  mirror  damage  with  high-power  laser  irradiance. 

Effective  acquisition  schemes  for  searching  geosynchronous  orbits  are  limited 
by  the  lengthy  pulse  round-trip  times  from  such  long  ranges.  However,  by  using 
two  telescopes  with  differing  primary  apertures  in  a  transmit/receive  combination, 
acquisition  time  is  greatly  reduced  as  well  as  improving  system  sensitivity. 
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IV.  The  Infrared  Environment 


Introduction 

This  chapter  examines  the  major  environmental  characteristics  which  affect 
LADAR  reception  and  IR  propagation  in  the  atmosphere.  The  chapter  covers  four 
major  areas  which  describe  (1)  computer  models  of  the  atmosphere,  (2)  background 
solar  and  thermal  radiance  across  the  IR  region.  (3)  atmospheric  transmission  of 
selected  IR  wavelengths,  and  (4)  other  effects  on  propagation  such  as  turbulence  and 
air  breakdown.  Although  the  unsuitability  of  current-technology  NdtYAG  and  DF 
lasers  for  LADAR  applications  was  highlighted  in  Chapter  III,  both  wavelengths  are 
included  in  initial  radiance  and  transmittance  calculations  for  comparison  purposes. 


Atmospheric  Models 

Several  computer  programs  exist  which  approximate  the  structure  and  char¬ 
acteristics  of  the  atmosphere  in  order  to  calculate  the  amount  of  scattering  and 
absorption  of  optical  radiation.  Perhaps  the  most  widely  used  program  in  the  sci¬ 
entific  community  is  the  LOWTRAN  atmospheric  transmittance/radiance  model 
developed  by  the  Air  Force  Geophysics  Laboratory.  The  currently  available  version 
of  the  program  is  L0WTRAN7,  and  was  developed,  according  to  Kneizys  et  al..  to 
provide  a  quick  method  of  estimating 

. . .  atmospheric  transmittance,  atmospheric  background  radiance,  sin¬ 
gle  scattered  solar  and  lunar  radiance,  direct  solar  radiance,  and  mul¬ 
tiple  scattered  solar  and  thermal  radiance.  The  spectral  resolution  of 
the  model  is  20  cm-1  (full  width  at  half-maximum)  in  steps  of  i  cm-1 
from  0  to  50,000  cm-1  (0.2  pm  to  infinity).  A  single-parameter  band 
model  is  used  for  molecular  line  absorption  and  the  effects  of  molecular 
continuum-type  absorption;  molecular  scattering,  aerosol  and  hydrome¬ 
teor  absorption  and  scattering  are  included.  Refraction  and  earth  curva¬ 
ture  are  considered  in  the  calculation  of  the  atmospheric  slant  path  and 
attenuation  amounts  along  the  path.  (46:2) 
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The  L0WTRAN7  program  contains  six  representative  atmospheric  models  for 
different  geographies  and  seasons,  and  also  includes  several  representative  aerosol, 
cloud,  and  rain  models  (40:18-46;  46:4).  However,  the  use  of  LOWTRAN  for  cal¬ 
culating  transmittance  of  laser  radiation  is  inadequate  due  to  insufficient  resolution 
for  monochromatic  calculations.  For  example,  calculation  of  transmittance  over  a 
"ertical  path  from  ground  (sea  le%-e!)  to  space  for  C02  laser  radiation  at  10.6  /zm 
(944.19402961  un'1!  is  calculated  using  LOWTRAN  to  be  about  0.87  for  clear  skies, 
mid-latitude  summer  atmosphere,  and  no  aerosol  attenuation.  Since  this  value  is  av¬ 
eraged  over  a  20  cm-1  interval,  the  full  effect  of  line  absorption  bv  atmospheric  C02 
is  not  included. 

A  detailed  calculation  of  monochromatic  transmittance  can  only  be  obtained 
by  considering  the  absorption  contribution  from  each  spectral  line  of  the  atmospheric 
constituents  (40:18-37).  Such  a  ‘line-by-line’  calculation  is  possible  using  FASCODE 
(Fast  Atmospheric  Signature  Code),  which  is  another  program  developed  by  the 
Air  Force  Geophysics  Laboratory.  Major  FASCODE  features  are  similar  to  LOW- 
TRAN,  but  the  algorithm  is  optimized  for  the  “accelerated  line  by  line  calculation 
of  spectral  transmittance  and  radiance  for  atmospheric  problems’’  (15:1).  Using  the 
FASCOD2  version  of  FASCODE,  the  calculated  transmittance  value  for  C02  laser 
radiation  yields  a  ground-to-space  transmittance  value  of  only  0.32,  with  the  same 
options  selected  as  for  the  L0WTRAN7  calculation.  Clearly,  line-by-line  calcula¬ 
tion  is  required  to  compute  accurate  transmittance  values  for  the  wavelengths  under 
study. 

In  addition  to  modeling  the  atmospheric  transmittance  of  radiation,  both 
LOWTRAN7  and  FASCOD2  include  the  capability  to  calculate  background  sky- 
radiance  for  a  selected  look  angle.  In  particular,  LOWTRAN7  includes  a  multiple 
scattering  parameterization  to  allow  calculation  of  background  solar  and  thermal 
radiance.  The  estimated  error  associated  with  this  parameterization  is  less  than 
20%,  when  all  possible  look  angles  are  considered  (46:4).  FASCOD2,  however,  does 
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not  include  multiple  scattering  effects,  and  only  uses  a  single  scattering  nodcl  when 
calculating  radiance  (15:1).  Although  the  single  scattering  model  is  a  good  approx¬ 
imation  when  look  angles  close  to  the  sun  aie  selected,  FASC0D2  radiance  calcula¬ 
tions  generally  underestimate  sky  radiance  (40:18-58).  Therefore,  the  LOWTRA.Y/ 
program  was  selected  for  all  radiance  calculations  in  this  analysis. 

Background  Radiance 

Sources  of  background  radiance  in  a  LADAR  system  include  direct  solar  ra¬ 
diance.  scattered  solar  radiance,  thermal  sky  emission,  and  thermal  emission  from 
receiver  optics.  To  simplify  calculations,  direct  viewing  of  the  sun  will  not  be  consid¬ 
ered,  and  cryogenic  cooling  of  receiver  optics  will  be  assumed  so  that  optics  emission 
is  negligible.  In  terms  of  viewing  conditions,  worst-case  background  radiance  occurs 
in  daytime  due  to  solar  radiation  and  thermal  sky  emission.  Furthermore,  sunlit 
cloud  is  detrimental  for  near  IR  wavelengths,  while  sky  emission  is  detrimental  at 
longer  wavelengths  (87:3-71).  These  factors  will  be  examined  using  LOWTRAN7. 

Several  options  are  available  for  calculating  radiance  using  LC'Vv  FRAN?,  with 
choices  between  several  reference  atmospheres,  aerosol  models,  weaker  conditions, 
and  viewing  geometries.  Models  for  the  atmosphere  and  aerosol  can  be  selected  with 
knowledge  of  the  proposed  sensor  site.  Since  co-location  with  GEODSS  is  a  logical 
choice  for  siting  of  a  LADAR  sensor  network,  characteristics  of  the  four  existing 
GEODSS  sites  are  assumed  to  be  representative  of  likely  LADAR  locations.  The 
geographic  locations  and  elevations  above  sea  level  for  the  four  GEODSS  sites  are 
shown  in  Table  4. 

The  available  LOWTRAN7  atmospheric  models  are  tropical  (15°  N),  mid¬ 
latitude  (45°  N)  summer  and  winter,  subarctic  (60°  N)  summer  and  winter,  and 
the  1976  US  Standard  atmosphere.  The  aVailable  aerosol  models  for  the  boundary 
layer  (0-2  km)  include  rural,  maritime ,  urban,  tropospheric,  advective  fog,  radiative 
fog,  and  desert.  Aerosol  models  for  the  tropospheric  (2-10  km)  and  stratospheric 
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Table  4.  GEODSS  Site  Characteristics  (Extracted  from  (47:57;  30:111-22)) 


GEODSS  Site 

Latitude 

East  Longitude 

Elevation  (km) 

Socorro,  NM 

33.3 

253.3 

1.53 

Maui,  HI 

20.7 

203.7 

3.04 

iaegu,  ooudn  Korea 

35.7 

128.6 

0.75 

Diego  Garcia 

-8.0 

73.0 

sea  level 

(10-30  km)  layers  are  usually  determined  by  the  season  of  the  model,  while  aerosols 
above  30  km  are  assigned  a  default  profile  unless  the  effect  of  volcanic  aerosol  ex¬ 
tinction  is  sought. 

The  viewing  geometry  used  foi  radiance  calculations  appears  in  Figure  5.  and 
shows  that  three  angles  must  be  specified  for  a  particular  viewing  event.  The  look 
angle  is  determined  by  the  satellite's  position  relative  to  the  sensor,  the  solar  zenith 
angle  by  the  sun's  relative  position,  and  the  relative  azimuth  angle  by  the  relative 
positions  of  satellite  and  sun  with  respect  to  the  sensor. 


Figure  5.  Geometry  for  Calculation  of  Solar  and  Thermal  Scattered  Radiance 
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Using  the  terminology  of  Figure  5.  two  cases  of  viewing  geometry  were  analyzed 
for  IR  background  radiance  using  L0WTRAN7.  Both  cases  were  calculated  using 
a  mid- latitude  summer  atmosphere,  rural  aerosol  with  23  km  visibility,  solar  zenith 
angle  of  0°,  and  a  sensor  location  ai  sea  level.  The  first  case  was  calculated  for  a 
look  angle  of  70°  against  a  clear-sky  background,  while  the  second  case  used  a  look 
angle  of  10°  against  a  cloudy  background.  Figure  6  shows  the  resultant  background 
level  of  spectral  radiance  for  both  cases  over  the  wavelength  band  of  1  /im  to  3  /im. 
The  Nd:YAG  laser  operating  wavelength  is  also  annotated. 


Figure  6.  Background  Spectral  Radiance  in  the  1-3  /im  Band,  Calculated  Using 
LOWTRAN7  for  two  Cases:  70°  Look  Angle  From  Sea-level  Against  a 
Clear  Sky,  and  10°  Look  Angle  From  Sea-level  Against  a  Cloudy  Sky 

Radiance  levels  for  the  additional  cases  of  70°  look  angle/cloudy  sky,  and  10° 
look  angle/clear  sky  were  found  to  lie  between  the  two  extreme  cases  shown  in 
Figure  6,  but  were  not  shown  for  reasons  of  plot  clarity.  Figure  6  shows  a  general 
trend  for  short  wavelengths  of  increased  levels  of  radiance  at  look  angles  close  to  the 
sun,  along  with  increased  background  radiance  from  a  cloudy  field-of-view.  At  the 


Nd:YAG  wavelength,  background  spectral  radiance  for  all  calculated  cases  is  vi-rv 
high,  ranging  from  17  \Y  m_2sr_1/mi_1  to  S3  Wnf’sf’/im'1. 

Figure  7  shows  the  spectral  radiance,  caU-  Ui-ed  under  the  same  conditions,  for 
the  3-12  p m  region,  annotated  with  DF  and  C 0 lase--  wavelengths.  In  comparison 
to  the  1-3  pm  band,  background  radiance  in  the  3-12  pm  region  is  reduced  by  about 
one  order  of  magnitude  (note  the  scale  change  between  figures).  Although  the  eases 
for  70°  look  angle/cloudy  sky,  and  10°  look  angle/clear  sky  were  again  not  shown  for 
clarity,  the  presence  of  cloud  at  each  look  angle  was  found  to  affect  radiance  levels 
only  slightly. 

.  70°.  Clear  !  j 

-  10°,  Cloud  ,  | 

|  Lasers 

i 

Solar  Zenith  =  0°  i 


3  4  5  6  7  8  9  10  11  12 

Wavelength,  pm 

Figure  7.  Background  Spectral  Radiance  in  the  3-12  pm  Band,  Calculated  Using 
LGWTRAN7  for  two  Cases:  70°  Look  Angle  From  Sea-level  Against  a 
Clear  Sky,  and  10°  Look  Angle  From  Sea-level  Against  a  Cloudy  Sky- 

Compared  to  the  small  look-angle  case,  radiance  levels  in  the  far-IR  were  found 
to  be  higher  at  large  look  angles  due  to  greater  thermal  sky  emission.  For  a  300  °K 
sky,  this  thermal  emission  peaks  at  around  10  pm.  The  crossover’  Irom  domina¬ 
tion  by  scattered  solar  radiance  to  domination  by  thermal  sky  radiance  is  at  about 
4.3  pm.  The  DF  and  frequency-doubled  CO2  lasers  radiate  near  the  crossover  point. 


at  wavelengths  which  correspond  to  the  lowest  level  of  background  radiance  for  the 
five  lasers  under  study.  Average  radiance  levels  at  DF  and  frequency-doubled  C02 
wavelengths  are  only  about  0.3  Wm'V^tm'1  and  0.7  \V  m_3sr-1/zrrr1 .  respec¬ 
tively.  while  13C02  and  I2C02  wavelengths  encounter  radiance  levels  which  are  an 
order  of  magnitude  larger  at  2.5  to  5  Wm^sf’/im'1.  depending 

on  look  angle. 

Atmospheric  Transmission 

Accurate  FASCODE  calculations  for  atmospheric  transmittance  of  monochro¬ 
matic  radiation  require  precise  specification  of  the  laser  wavelength  to  adequately 
account  for  the  strong  influence  of  line  absorption.  However,  among  the  frequencv- 
doubled  and  isotopically-substituted  C02  LADAR  systems  reported  in  the  literature, 
different  laser  transitions  are  used  from  system  to  system.  For  example,  different 
frequency-doubled  systems  reported  by  Harney  (27:99)  and  Miller  et  al.  (56:6)  use 
the  R(20)  and  R ( 1 4 )  C02  transitions,  respectively.  Among  13C02  systems,  several 
transitions  are  reportedly  used,  including  P(24)  (27:99).  P(20)  (88:82;  31:5),  and 
P ( 1 8 )  transitions  (57:103). 

Due  to  the  diversity  of  laser  transitions  in  use  for  the  C02-variant  systems, 
four  different  13C02  transitions  were  selected  for  analysis,  along  with  five  separate 
i2C02  transitions  for  the  frequency-doubled  system.  Table  5  shows  the  resultant  laser 
wavelengths  selected  for  FASCOD2  ca'culations.  All  C02  transition  frequencies  were 
taken  from  Witteman  (86:24-26). 

Although  all  I2C02  systems  found  in  the  literature  use  the  P(20)  line,  no  re¬ 
ports  of  single-line  DF  systems  were  found.  For  FASCOD2  calculations,  a  DF  wave¬ 
length  of  3.83  pm  was  selected  as  a  representative  value,  after  multiple  calculations 
showed  no  nearby  strong  absorption  lines  (implying  negligible  amounts  of  DF  in  the 
atmosphere).  Finally,  the  strongest  Nd:YAG  transition  at  1.06  /im  was  selected  to 
complete  the  list  of  candidate  LADAR  wavelengths. 


Table  5.  Laser  Transition  Wavelengths  used  in  FASC'OD'2  Calculations 


Laser 

Transition  Wavelength  (pm)  Frequency  (cm  l) 

13co2 

P  ( 24 )  11.19354451  893.37206747 

P  f  22 )  11.17130499  895.15056754 

P(  20 )  11.14939715  896.90947969 

P(  18)  11.12781727  898.64883264 

12co2 

P(20)  10.59104346  944.19402961 

12co2 

( frequency-doubled) 

R(  14)  4.65269281  2149.29298016 

R(  16)  4.64689275  2151.97564042 

R(18)  4.64122185  2154.60504026 

R(20)  4.63567900  2157.18128846 

R(22)  4.63026305  2159.70451160 

DF  average  3.83  2610.96605744 

Nd:YAG  4F3/2  -  4I;i/2  1.0642  9396.72993841 


Reduction  of  the  number  of  candidate  13C02  wavelengths  was  carried  out  by 
performing  preliminary  transmittance  calculations  for  all  selected  transitions  over  a 
vertical  path  from  ground  (sea  level)  to  space  through  various  reference  atmospheres. 
Atmospheric  profiles  for  the  four  GEODSS  locations  were  modeled  using  either  mid¬ 
latitude  or  tropical  reference  atmospheres,  with  corresponding  aerosol  profiles  mod¬ 
eled  as  rural  or  maritime,  respectively.  All  four  13C02  transitions  exhibited  similar 
transmittance  characteristics  through  all  atmospheric  models;  hence,  the  P(22)  tran¬ 
sition  was  arbitrarily  chosen  as  a  representative  13C02  operating  wavelength  for  all 
subsequent  calculations. 

Table  5  was  further  reduced  by  performing  a  second  comparison  between  tran¬ 
sitions  of  frequency-doubled  12002  for  the  same  transmission  path.  Figure  8  shows 
the  FASCOD2  results  for  three  cases;  mid-latitude  summer  model  with  rural  aerosol, 
mid-latitude  winter  model  with  rural  aerosol,  and  tropical  model  with  maritime 
aerosol.  Visibility  of  23  km  and  clear  skies  were  selected  for  all  cases. 
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Figure  8.  Atmospheric  Transmittance  of  Frequency- Doubled  CO2  Laser  Transitions 
over  a  Vertical  Path  from  Ground  to  Space  for  three  cases:  Mid- Latitude 
Summer  and  Winter  Atmospheres  with  Rural  Aerosol,  and  Tropical  At¬ 
mosphere  with  Maritime  Aerosol 

Of  the  three  cases,  poorest  transmission  occurs  through  a  tropical  atmosphere, 
with  best  transmission  through  a  mid-latitude  winter  atmosphere.  This  variation 
may  be  attributed  to  different  humidity  profiles  in  each  model,  which  cause  different 
amounts  of  water-vapor  attenuation  of  the  laser  radiation.  An  additional  feature 
of  Figure  8  is  the  strong  absorption  line  which  occurs  in  the  vicinity  of  the  R(18) 
transition,  so  that  transmittance  in  all  atmospheric  cases  is  virtually  extinguished 
around  4.64 1  /xm.  Transmittance  at  all  four  remaining  wavelengths  is  reasonably 
high;  consequently,  the  R(20)  wavelength  was  arbitrarily  chosen  as  representative 
for  subsequent  calculations. 

A  comparison  of  atmospheric  transmittance  characteristics  was  then  performed 
between  the  reduced  lis*  of  five  candidate  wavelengths  from  Table  5,  for  the  same 
transmission  path  and  models,  with  the  results  plotted  in  Figure  9.  As  expected, 
strong  self-absorption  is  evident  at  the  12C02  wavelength,  while  the  remaining  four 
laser  wavelengths  display  generally  high  transmittance.  Figure  9  also  shows  the 
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adverse  effect  that  water  vapor  has  on  transmission  in  the  8-12  \i m  band,  with 
transmittance  through  a  tropical  atmosphere  only  about  50%  of  the  value  calculated 
for  a  mid-latitude  winter  atmosphere. 


Xd-.YAG  DF  CO-i  doubled.  UC02  UC02 

Laser  Wavelength  (Not  to  Scale) 


Figure  9.  Atmospheric  Transmittance  of  Five  Selected  Laser  Wavelengths  over  a 
Vertical  Path  from  Ground  to  Space  for  three  cases:  Mid-Latitude  Winter 
and  Summer  Atmospheres  with  Rural  Aerosol,  and  Tropical  Atmosphere 
with  Maritime  Aerosol 

The  effect  of  different  atmospheric  models  on  DF  transmittance  is  seen  from 
Figure  9  to  be  only  slight.  These  excellent  transmittance  characteristics  of  DF 
radiation,  combined  with  the  high  power  possible  from  multi-line  DF  operation, 
are  strong  reasons  for  the  continued  development  of  high-efficiency,  single-line  DF 
LADAR  transmitters. 

Look  A  ngle.  The  influence  of  look  angle  on  transmittance  for  the  likely  range  of 
geographic  and  atmospheric  conditions  was  then  examined.  Viewing  conditions  were 
modeled0  in  FASCOD2  by  selecting  GEODSS  locations  at  Socorro  and  Diego  Garcia 
as  indicative  of  two  extreme  geographic/atmospheric  combinations.  Accordingly,  the 
two  cases  considered  were:  Socorro  —  slant  path  from  1.53  km  altitude  to  space, 
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through  a  mid-latitude  winter  atmosphere  and  rural  aerosol  with  23  km  visibility  — 
and  Diego  Garcia  —  slant  path  from  sea  level  to  space,  through  a  tropical  atmosphere 
and  maritime  aerosol  with  23  km  visibility.  Additionally,  to  reduce  computational 
effort,  wavelengths  of  only  the  three  C02-variant  lasers  were  examined. 

Figures  10  and  11  show  the  results  of  FASCOD2  calculations  for  the  three  laser 
wavelengths.  Comparing  transmittance  at  large  look  angles  with  vertical  transmis¬ 
sion,  both  figures  show  higher  relative  attenuation  losses  with  12C02  than  with  the 
13C02  and  frequency-doubled  lasers.  For  example,  12C02  transmittance  from  So¬ 
corro  at  a  60°  look  angle  is  40%  lower  than  vertical  transmittance.  Conversely, 
frequency-doubled  C02  transmittance  at  60°  from  Socorro  falls  off  by  5%,  while 
13C02  transmittance  is  reduced  by  only  3%.  In  the  Diego  Garcia  case,  the  transmit¬ 
tance  degradation  at  60°  look  angle  for  12C02,  frequency-doubled  C02,  and  13C02 
is  57%,  21%,  and  17%,  respectively. 

Since  two-way  atmospheric  transmittance  is  required  for  LADAR  signals,  the 
degradation  in  one-way  transmittance  at  large  look  angles  has  drastic  effects  on 
LADAR  signals.  For  example,  the  57%  reduction  in  12C02  laser  transmittance  from 
Diego  Garcia  at  a  60°  look  angle  has  the  effect  of  attenuating  LADAR  signals  by  a 
further  82%  when  compared  to  two-way  vertical  transmission.  This  sharp  reduction 
in  two-way  12C02  transmittance  at  large  look  angles  would  limit  LADAR  coverage 
to  a  small  range  of  look  angles  unless  extra  transmitter  power  was  available  to 
compensate  for  the  higher  signal  attenuation.  Therefore,  the  12C02  laser  is  rejected 
as  a  LADAR  transmitter  candidate  in  favor  of  the  13C02  and  frequency-doubled 
lasers. 

Cloud  and  Rain.  The  effect  on  transmittance  of  clouds  and  rain  can  also  be  an¬ 
alyzed  using  FASCOD2.  However,  software  anomalies  associated  with  the  available 
copy  of  both  FASCOD2  and  LOWTRAN7  prevented  execution  using  these  options. 
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Figure  10.  Atmospheric  Transmittance  at  Socorro  —  Slant  Path  from  1.53  km  Al¬ 
titude  to  Space,  Through  a  Mid-Latitude  Winter  Atmosphere  and  Rural 
Aerosol  with  23  km  Visibility 


Figure  11.  Atmospheric  Transmittance  at  Diego  Garcia  —  Slant  Path  from  Sea- 
Level  to  Space,  Through  a  Tropical  Atmosphere  and  Maritime  Aerosol 
with  23  km  Visibility 
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Nonetheless,  an  estimate  of  the  effects  of  cloud  and  rain  may  be  gained  by  considering 
the  algorithms  used  in  LOVVTRAN  and  FASCODE  to  determine  cloud  and  rain 
attenuation.  Although  details  of  the  L0WTRAN7  algorithms  were  unavailable  for 
this  analysis,  L0WTRAN6  algorithms  were  used  instead. 

The  primary  model  used  in  L0WTRAN6  for  cloud  attenuation  is  the  cirrus 
cloud  model.  Other  models  for  cumulus  and  stratus  profiles  have  been  added  to 
L0WTRAN7  and  FASCOD‘2,  but  details  of  these  were  also  unavailable  for  this 
study.  Based  on  empirical  data,  attenuation  by  cirrus  has  been  found  to  be  ‘‘pro¬ 
portional  to  the  cloud  thickness  and  independent  of  wavelength  from  ...0.317 
to  10  /.i m”  (45:58).  The  expression  used  to  calculate  cirrus  transmittance  is  given  by 
Kneizys  et  al.  as  (45:66) 

r  =  e~°‘1412  (3) 

where  L  is  the  cloud  thickness  in  km. 

Based  on  measured  cloud  statistics,  cirrus  thickness  was  found  by  Kneizys  et  al. 
to  follow  a  truncated  log-normal  distribution,  with  a  median  cirrus  thickness  of 
1.0  km  (45:60).  The  average  height  of  the  cirrus  cloud  base  was  also  defined  for  trop¬ 
ical  and  mid-latitude  winter  atmospheres  as  11  km  and  8  km,  respectively  (45:67). 
Using  the  median  thickness  value  in  Equation  (3)  yields  an  average  IR  cirrus  trans¬ 
mittance  of  0.87  for  vertical  propagation.  Assuming  a  plane-parallel  atmosphere  and 
an  ideal  cirrus  cloud  of  constant  thickness  and  uniform  density,  Equation  (3)  may 
be  modified  to  include  the  effect  of  look  angle  as 

r  =  e_0-14  secif,  L2  (4) 

where  xb,  is  the  sensor  look  angle. 

The  assumption  of,a  plane-parallel  atmosphere  allows  the  use  of  sensor  look 
angle  to  approximate  the  zenith  angle  at  which  laser  radiation  enters  the  cloud. 
Considering  the  range  of  cloud-base  altitudes  under  study,  exact  calculations  show 
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that  the  approximation  of  the  cloud-entry  angle  is  in  error  by  less  than  0.4%  for 
look  angles  as  large  as  70°.  Equation  (4)  also  neglects  refraction,  which  is  a  valid 
assumption  for  cloud-entry  zenith  angles  of  less  than  SO0  (45:13). 

The  rain  attenuation  algorithm  used  in  LOWTRAN6  relates  transmittance  to 
rain-rate,  and  is  also  independent  of  wavelength  over  the  IR  region.  The  expression 
used  to  calculate  rain  transmittance  is  given  by  Kneizys  et  al.  as  (45:69) 


—0.365  s  r0'63 


where  s  is  the  path  length,  in  km,  through  rain  falling  at  a  rate  r.  in  mm/hr. 

Calculation  of  the  path  length,  s,  can  be  written  in  terms  of  look  angle  and  alti¬ 
tude  of  the  cloud  tops  causing  the  rain,  z:,  as  s  =  :c  secrl>3.  Therefore.  Equation  (5) 
may  be  re-written  as 

r  _  e~ 0.365  2.  sec ip,  r0  63  ^ 

Rain  rate  is  categorized  in  LOWTRAN6  and  FASCOD2  as  drizzle  (2  mm/hr). 
light  rain  (5  mm/hr),  moderate  rain  (12.5  mm/hr),  heavy  rain  (25  mm/hr),  and  ex¬ 
treme  rain  (75  mm/hr)  (46:25-26).  Figure  12  shows  Equation  (6)  plotted  for  the 
Socorro  site  through  drizzle,  light  rain,  and  moderate  rain,  and  also  shows  Equa¬ 
tion  (4)  plotted  for  median  cirrus  conditions.  Note  that  the  transmittance  values 
plotted  in  Figure  12  can  be  applied  to  all  IR  wavelengths,  but  oniy  include  the  at¬ 
tenuation  caused  by  cirrus  or  rain.  A  total  atmospheric  transmittance  value  for  the 
entire  ground-to-space  path  through  cloud  or  rain  can  be  obtained  by  multiplying 
the  appropriate  Figure  12  transmittance  value  by  the  corresponding  clear-weather 
transmittance  from  Figure  10  or  11. 

Figure  12  cirrus  transmittance  was  calculated  with  a  median  cirrus  profile 
existing  between  8  km  and  9  km  altitude.  For  the  rain  conditions,  rain  rate  was 
modeled  as  constant  between  ground  level  and  2  km  altitude  (470  m  above  the 
Socorro  site).  Attenuation  from  rain  clouds  was  not  included  in  the  rain  calculation 
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Figure  12.  Transmittance  Component  due  to  Cloud  and  Rain  from  Socorro  for 
Three  Cases:  Cirrus  Cloud  from  8  km  to  9  km  Altitude.  Drizzle  Origi¬ 
nating  from  2  km,  Light  Rain  from  2  km,  and  Moderate  Rain  from  2  km 


since  a  suitable  low-altitude  cloud  model  was  unavailable,  and  cloud  attenuation 
is  smaU  <'r>mpared  to  rain  attenuation.  The  corresponding  LOWTRAN6  model  for 
these  conditions  uses  a  gradually  decreasing  rain  rate  from  ground  level  up  to  a 
final  rain  altitude  of  2  km  or  3  km  altitude,  depending  on  rain  rate.  Therefore,  the 
conditions  used  in  Figure  12  account  for  higher  rain  attenuation  than  the  equivalent 
LOWTRAN6  model. 

Although  use  of  a  constant  rain-rate  model  is  a  simplification,  some  idea  of 
the  severity  of  rain  attenuation  at  IR  wavelengths  may  be  gained  from  the  resultant 
data.  A  rain  rate  ot  only  2  mm/hr  results  in  a  68%  reduction  in  total  transmittance 
for  a  vertical  path.  This  reduction  corresponds  to  a  90%  decrease  in  two-way  total 
transmittance. 

Attenuation  effects  at  higher  rain  rates  and  larger  look  angles  are  more  severe, 
and  indicate  that  LADAR  operation  through  rain  is  impractical  due  to  the  extreme 
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amount  01  attenuation  for  two-way  signal  propagation.  For  example,  if  a  13C02 
transmitter  is  used  from  Socorro  at  a  moderate  look  angle  of  30°  through  drizzle  total 
two-way  atmospheric  transmittance  amounts  to  only  7%.  For  lower-altitude  sites, 
such  as  Diego  Garcia,  transmittance  becomes  negligible  for  even  trace  precipitations. 

On  the  positive  side,  however,  attenuation  due  to  cirrus  appears  acceptable 
for  continued  LADAR  operation.  Using  the  same  example  of  a  13C02  signal  at  30°, 
total  two-way  transmittance  is  about  70%  for  Socorro,  and  50%  for  Diego  Garcia. 
Although  the  attenuating  effects  of  lower-altitude  clouds  such  as  stratus  and  cumulus 
have  not  been  considered,  operation  of  an  IR  LADAR  through  cirrus  cloud  seems 
practical. 

Other  Effects 

Turbulence.  Apart  from  scattering  and  absorption,  laser  transmission  through 
the  atmosphere  is  likely  to  be  influenced  by  several  other  effects  produced  by  at¬ 
mospheric  turbulence.  Turbulence  induces  random  variations  in  the  atmospheric 
index  of  refraction,  via  fluctuations  in  temperature,  pressure,  and  humidity.  Theory 
developed  by  Kolmogorov  models  the  turbulent  atmosphere  as  consisting  of  mobile 
parcels  of  air.  or  eddies ,  of  differing  indices  of  refraction  (37:6-9).  While  these  eddies 
range  in  size  from  sub-centimeter  at  the  earth's  surface,  to  global  scale  in  the  open 
atmosphere,  the  size  range  of  interest  in  laser  propagation  lies  from  a  few  millimeters 
to  a  few  meters,  known  as  the  inertial  sub-range  (87:6-9). 

Turbulence  effects  experienced  by  LADAR  signals  include  beam  steering,  image 
dancing,  beam  spreading,  spatial  coherence  degradation,  and  scintillation  (40:18-71: 
69:3295).  Beam  steering  occurs  when  the  laser  beam  diameter  is  smaller  than  the 
turbulence  eddies,  and  can  result  in  deviation  of  the  beam  so  that,  in  extreme  cases, 
the  target  or  receive,  aperture  is  completely  missed  (40:18-71).  Conversely,  when  the 
beam  diameter  is  larger  than  the  eddies,  beam  spreading  occurs,  resulting  in  lower 
irradiance  values  at  the  target  or  receiver. 


Turbulence-induced  modulation  in  the  angle-of-arrival  of  the  received  wave 
can  also  occur,  resulting  in  shifting  of  the  image  focus,  or  image  dancing.  Smaller 
turbulence  eddies  cause  phase  changes  across  the  wavefront,  resulting  in  loss  of 
spatial  phase  coherence,  and  also  induce  interference  effects  within  the  beam  which 
cause  random  power  fluctuations,  or  scintillation ,  across  the  wavefront. 

Beam  steering  and  image  dancing  can  both  be  reduced  by  increasing  the  re¬ 
ceiver  aperture.  Since  prospective  LADAR  locations  such  as  the  GEODSS  sites  are 
characterized  by  low  turbulence,  the  most  serious  turbulence  effects  remaining  for 
coherent  IR  reception  are  scintillation  and  coherence  losses.  Description  of  these 
turbulence  effects  is  best  expressed  in  terms  of  the  refractive-index  structure  coeffi¬ 
cient,  Cd.  Values  of  vary  widely  according  to  prevailing  atmospheric  conditions, 
and  are  a  complex  function  of  altitude  and  wind  speed  (87:6-14;  40:18-72).  Several 
models  of  have  been  developed  from  empirical  data,  with  one  of  the  simplest 
given  by  (87:6-14) 


f  1.5  x  10~13 


Cl(z)  =  { 


below  20  km  above  sea  level 
above  20  km  above  sea  level 


(7) 


where  z  is  altitude  above  local  ground,  in  meters. 

Scintillation  effects  at  the  receiver  are  described  by  the  variance  in  the  log- 
amplitude  of  the  received  plane  wave,  given  by  (27:92) 

=  0.56  k 7/6  f  Cfc)  (3  -  =)5,s  ( j)  iz  (8) 

where  s  is  the  path  length. 

Scintillation  also  affects  the  transmitted  laser  beam,  which  can  cause  non-uniform 
irradiance  of  the  target,  depending  on  its  size.  The  combined  effect  of  scintillation 
over  both  transmit  and  receive  paths  (which  may  be  correlated  (14:514))  makes 
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estimation  of  cr  difficult.  Typically,  a  signal-to-noise  ratio  ’scintillation  margin'  is 
added  to  design  calculations  to  account  for  power  fluctuations  (87:23-16). 

Signal  coherence  losses  are  usually  described  in  terms  of  r0,  the  transverse 
phase  coherence  length ,  which  measures  the  lateral  distance  over  which  the  radiation 
remains  constant  (-10:18-71;  87:6-29).  In  the  transmitter  case,  when  the  primary 
aperture  exceeds  r01  "the  effective  transmitter  beamwidth  will  be  turbulence-limited 
rather  than  diffraction-limited”  (69:3295);  for  the  receiver  case,  when  the  primary 
aperture  exceeds  r0,  heterodyne  efficiency  is  reduced,  thereby  reducing  signal-to-noise 
ratio  (87:6-30).  Harney  expresses  the  transverse  coherence  length  for  L lie  Iransinitter- 
to-target  path,  rot,  as  (27:92) 

r  fzt  /  -\5/3  1-3/5 

rot  =  2.91  k2  C2n(z)  [l-f -J  dz  (9) 

where  k  —  2 rr / A  ,  and  zt  is  the  target  altitude  above  local  ground. 

Similarly,  Harney  gives  the  transverse  coherence  length  for  the  target-to-receiver 
path,  ror,  as  (27:92) 

r  -  -1  -3/5 

r„=  2.9 \k3  f"cl(z)  (i-V  dz  (10) 

Jo  \Zt/ 

Using  the  simplified  model  of  C 2  from  Equation  (7)  in  Equations  (9)  and 
(10)  allows  a  simplified  calculation  of  r„  for  selected  LADAR  site  elevations  and 
wavelengths.  Figure  13  shows  the  results  for  rot  and  ror,  calculated  over  the  1-13  /im 
band  for  the  Diego  Garcia  site  elevation.  The  figure  shows  that  coherence  lengths  for 
the  receive  path  are  generally  about  five  times  longer  than  for  the  transmit  path.  At 
the  ,3C02  wavelength,  the  ra(  value  —  and  hence,  maximum  transmitter-aperture 
diameter  —  is  about  0.9  m,  while  the  maximum  receiver-aperture  diameter  is  about 
4.5  m.  For  the  frequency-doubled  C02  wavelength,  maximum  aperture  diameters 
are  one-third  of  the  size,  at  0.3  m  and  1.5  m,  respectively. 
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Figure  13.  Transverse  Phase  Coherence  Length  for  Transmit  and  Receive  Paths 
versus  Wavelength 

Site  elevations  of  up  to  2  km  yield  virtually  the  same  curves  since  the  simplified 
expression  used  for  corresponds  to  a  turbulence  profile  commencing  from  local 
ground  rather  than  sea  level.  Thus,  an  increase  in  site  elevation  causes  the  20  km 
altitude  turbulence-cutoff  to  reduce  turbulence  from  the  ‘weak’  higher  end  of  the 
model  profile.  In  reality,  reduced  thermal  effects  and  increased  wind  at  higher  site 
elevations  leads  to  lower  turbulence  than  that  predicted  by  this  simplified  model. 
Therefore,  higher  site  elevations  wiii  generally  ailow  larger  values  of  r0.  Also,  since 
turbulence  is  random  in  nature,  there  will  always  be  a  finite  probability  chat  coher¬ 
ence  lengths  will  exceed  the  calculated  values,  no  matter  how  accurate  the  model. 
Wolfe  et  al.  state  that,  typically,  atmospheric  turbulence  limits  the  transmitter  beam 
divergence  to  about  10  /rrad  at  a  wavelength  of  10  /jm  (87:23-12). 

High-Power  Beam  Effects.  Two  additional  propagation  factors  need  to  be  con¬ 
sidered  when  the  transmitted  laser  beam  is  of  high  power.  The  first  factor  is  thermal 
blooming ,  which  occurs  when  the  beam  contains  enough  power  to  significantly  heat 


the  atmosphere  (87:6-35).  Although  serious  for  high-power  CW  signals,  thermal 
blooming  may  become  significant  in  pulsed  systems  which  use  long  pulse  widths  or 
high  repetition  rates.  The  second  high-power  effect,  air  breakdown ,  occurs  when  the 
laser  pulse  interacts  with  aerosol  particles  which  lower  the  air-breakdown  threshold. 
This  phenomenon  depends  on  particle  size,  beam  diameter,  and  pulse  width:  for 
example,  at  10.6  p m,  breakdown  has  been  showm  to  occur  at  an  irradiance  level 
of  100-200  MW/cm2,  which  is  attributed  to  aerosol  particles  of  4-10  pm  diameter 
(5:1421). 

Summary 

Analysis  of  the  environment  in  which  a  deep  space  IR  LADAR  must  oper¬ 
ate  has  highlighted  several  factors  which  influence  wavelength  selection,  including 
atmospheric  transmission,  background  radiance,  and  turbulence  effects.  TV  LOW- 
TRAN7  and  FA8COD2  computer  programs  were  used  to  compare  transmittance 
and  radiance  conditions  at  Socorro  and  Diego  Garcia.  Nd:YAG  and  DF  laser  wave¬ 
lengths  were  included  in  the  comparison  to  demonstrate  their  excellent  transmission 
performance.  However,  due  to  the  poor  output  characteristics  of  Nd:YAG  and  DF. 
the  final  transmitter  comparison  w'as  confined  to  C02-variant  lasers. 

The  final  comparison  showed  that  poorest  transmission  occurs  through  the 
tropical  atmosphere  at  Diego  Garcia,  with  transmittance  decreasing  with  larger  look 
angles  in  all  cases.  In  particular,  the  P(22)  wavelength  of  13C02  at  11.17  pm  exhib¬ 
ited  the  best  transmittance  under  all  conditions,  closely  followed  by  the  frequency- 
doubled  R(20)  wavelength  of  12C02  at  4.64  pm.  These  two  laser  candidates  represent 
the  best  examples  from  the  3-5  pm  and  8-12  pm  bands  in  terms  of  transmittance.  All 
other  candidate  C02  laser  wavelengths  were  rejected  as  potential  LADAR  transmit¬ 
ters  due  to  poor  atmospheric  transmittance  over  the  range  of  conditions  examined. 

The  effects  of  cloud  and  idin  were  found  to  be  independent  of  wavelength, 
based  on  the  cloud  and  rain  models  in  LOWTRAN’6.  Transmittance  was  found 
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to  be  severely  degraded  by  even  a  small  amount,  of  rain:  hence,  LADAR  operation 
through  rain  appears  impractical.  However,  LADAR  operation  through  average 
cirrus  conditions  was  shown  to  be  theoretically  feasible,  with  only  minor  attenuation 
effects. 

Background  radiance  levels  at  medium  to  far-IR  wavelengths  were  determined 
to  be  generally  low  compared  to  levels  at  shorter  wavelengths.  Of  the  two  laser 
candidates,  radiance  at  4.64  //m  is  much  lower  than  at  11.17  //m. 

Several  effects  of  turbulence  on  signal  propagation  were  also  examined,  with 
scintillation  and  coherence  degradation  highlighted  as  the  most  serious  effects  for 
long-range  LADAR  applications.  Compensation  of  scintillation  can  be  adequately 
performed  by  the  addition  of  a  scintillation  margin  to  system  signal-to-noise  calcu¬ 
lations.  Loss  of  phase  coherence,  on  the  other  hand,  may  be  minimized  by  limiting 
the  size  of  the  transmitter  and  receiver  apertures  to  the  value  of  r0,  the  transverse 
phase  coherence  length.  r0  is  wavelength-dependent  and  generally  allows  larger  re¬ 
ceiver  apertures  than  transmitter  apertures.  Figure  13  provides  rough  estimates  of 
aperture  size  limits  using  a  simplified  atmospheric  model. 
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V.  Parametric  Analysis 


Introduction 

The  design  of  a  ground-based  LADAR  which  is  designed  to  detect  small,  mov¬ 
ing,  deep-space  targets  while  ’looking'  through  a  turbulent,  hazy,  sunlit  atmosphere 
requires  high  performance  specifications.  Careful  design  of  all  sensor  elements  is 
essential,  from  telescope  mount  to  processing  electronics,  so  that  maximum  laser 
energy  is  deposited  on  target,  and  maximum  reflected  signal  is  collected  at  the  de¬ 
tector  for  processing.  Once  the  target  is  detected,  estimates  of  target  range,  angular 
position,  and  radial  velocity  are  obtained  from  the  detected  pulse.  Detection  of 
successive  pulse  returns  also  allows  estimation  of  target  angular  rates  and  lateral 
velocity,  which  may  be  used  to  derive  an  orbital  element  set  for  the  target. 

With  these  overall  design  goals  in  mind,  a  comprehensive  set  of  equations  is 
required  which  relate  laser,  telescope,  and  detector  characteristics  to  LADAR  per¬ 
formance  parameters.  This  equation  set  can  then  be  used  for  the  specification  and 
evaluation  of  a  candidate  ground-based  deep-space  LADAR  design;  the  LADAR 
configuration  used  as  the  basis  for  this  parametric  analysis  is  the  candidate  system 
outlined  in  Chapter  III,  with  the  transmitted  signal  consisting  of  a  train  of  XI  rect¬ 
angular  pulses.  Detection  analysis  concentrates  initially  on  a  single  square  detector 
element,  but  is  then  optimized  for  two-dimensional  FPA  detection  with  arbitrary 
pixel  and  gap  dimensions.  Analysis  of  the  heterodyne  process  examines  beam  mis¬ 
alignment  effects  and  also  considers  different  LO  characteristics  and  their  effect  on 
receiver  efficiency. 

Detection  Probability 

The  process  of  detecting  a  LADAR  signal  in  the  presence  of  noise  is  statistical 
in  nature  due  to  the  randomly  fluctuating  noise  component  (73:2-7).  This  detec¬ 
tion  process  is  conventionally  described  using  two  parameters  —  the  probability  of 
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detection ,  pj,  which  represents  the  probability  of  detecting  a  real  target,  and  the 
probability  of  false- alarm,  p;a,  which  represents  the  probability  of  mistakenly  iden¬ 
tifying  a  noise  fluctuation  as  a  real  target  (68:253;  73:2-7).  Both  parameters  are 
specified  prior  to  system  design  calculations,  and  are  examined  first  since  they  form 
the  basis  from  which  minimum  system  performance  specifications  are  derived. 

To  enable  further  analysis  of  the  detection  process.  Figure  14  depicts  a  block 
diagram  of  typical  LADAR  detection  circuitry  (6:38).  The  received  signal  is  het¬ 
erodyne  detected  before  passing  through  a  matched  filter.  The  output  of  the  filter 
is  envelope  detected  and  compared  to  a  threshold  reference,  providing  a  true/false 
output.  The  threshold  reference  is  set  by  combining  a  measurement  of  the  noise 
photocurrent  with  a  specified  p;a  (6:39). 


Figure  14.  LADAR  Detection  Circuitry  (Adapted  from  Bachman  (6:39)) 


The  sources  of  noise  in  a  heterodyne  LADAR  include  dark  current,  thermal 
noise,  background  photocurrent,  and  local  oscillator-induced  shot  noise  (41:141). 
However,  if  the  local  oscillator  (LO)  power  is  adjusted  so  that  the  weak  received 
signal  is  mixed  with  a  strong  LO  signal,  the  LO-induced  shot  noise  can  be  made  to 
dominate  all  other  noise  sources  (6:39).  Operation  of  the  LADAR  in  this  manner 
is  said  to  be  shot-noise  limited.  This  mode  of  operation  is  preferable  to  background 
limited  operation  since  background  noise  generally  fluctuates  in  an  unpredictable 


fashion  and  is  difficult  to  model.  Due  to  the  random  nature  of  shot  noise,  however. 
Gaussian  statistics  may  be  used  to  model  the  detection  process.  In  this  case  the 
noise  current  can  be  written  as 


Jo{t) 


n(t)  =  r(t)e 


where  n(t)  is  the  time-dependent  noise  current,  and  r(t)  and  o(t)  are  the  amplitude 
and  phase  of  the  current,  respectively,  at  time  t. 

The  Gaussian  joint  probability  density  function  of  r  and  d>  can  then  be  written 
as  (6:39) 

r 

p(r,4>)  =  7 — i  e  2<7"  (12) 

'  t  rr  ^ 


where  er2  is  the  average  noise  power. 

The  corresponding  noise  envelope  at  the  output  of  the  envelope  detector  with 
no  signal  present  can  then  be  found  by  integrating  Equation  (12)  to  remove  the 
phase  term  (8:59): 

r  _jd_ 

p{r)  =  —  e  (13) 

This  Rayleigh-distributed  probability  density  function  can  then  be  used  to 
provide  an  expression  for  the  probability  of  false  alarm,  p/a  (6:40): 


[OO 

fa  =  /  p(r)  dr  = 

Jr, 


e 


where  rt  is  the  decision  threshold. 

A  tynical  value  of  p>a  for  radar  detection  is  10-6.  However,  a  more  practical 
representation  of  false-alarm  occurrence  is  the  mean  time  between  false-alarms ,  tja. 
If  a  train  of  ,1/  received  pulses  is  detected  with  range  gates,  and  then  integrated  by 
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post-detection  circuitry,  the  mean  time  between  false-alarms  is  given  by  (73:2-17) 


Mtj 

P/a  (  1  —  ) 


15) 


where  tg  is  the  range  gate  on-time,  and  60r  is  the  fraction  of  the  time  for  which  no 
gates  are  open  (dead  time). 

For  example,  for  M  =  10  ,  t3  =  10  ps,  pja  =  10~6  ,  and  60r  =  0.99  .  the  theoreti¬ 
cal  mean  time  between  false-alarms  is  2.3  hours. 

From  Equation  (14),  /?/a  is  seen  to  depend  on  the  decision  threshold  and  the 
noise  power.  The  probability  of  detection,  on  the  other  hand,  is  influenced  by  the 
threshold,  noise  power,  and  signal  power.  Some  knowledge  of  the  reflected  signal 
characteristics  from  the  target  is  therefore  necessary  to  adequately  model  the  detec¬ 
tion  process.  Since  the  targets  under  consideration  have  generally  diffuse  surfaces, 
the  reflected  signal  is  described  by  Bachman  as 


.  . .  the  sum  of  the  field  scattered  from  each  of  the  individual  elemen¬ 
tary  scatterers,  each  of  which  has  a  uniformly  distributed  random  phase 
as  a  consequence  of  the  random  spatial  distribution  of  the  surface. 

Hence,  by  the  central  limit  theorem,  in-phase  and  phase-quadrature 
components  of  the  field,  as  well  as  the  photocurrents  they  induce,  have 
independent  Gaussian  distributions  with  a  common  variance,  a2d  [cr\  is 
the  average  signal  power  received  from  a  diffuse  target],  and  zero  mean. 
(6:40) 


Since  the  Gaussian  signal  distribution  is  added  to  the  Gaussian  noise  distri¬ 
bution,  the  envelope  detector  output  with  both  signal  and  noise  present  follows  a 
Rayleigh  probability  distribution  given  by 


p(r')  = 


(16) 


where  rr2  =  cr2n  +  a\. 
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In  a  similar  manner  to  Equation  (14),  the  probability  of  detection  can  be 
described  as 


Figure  15  provides  a  pictorial  representation  of  pja  and  p d.  The  value  of  a\  is 
easily  obtained  via  the  noise  measurement  process  shown  in  Figure  14,  but  the  value 
of  a\  relies  on  information  concerning  both  noise  and  signal  and  can  be  written  as 
(6:40) 

<y]  =  +  =  a2n{\+SNR)  (IS) 

where  SNR  is  the  signal-to-noise  power  ratio  from  a  diffuse  target. 


Figure  15.  Graphical  Representation  of  pd  and  pja 


SNR  is  now  the  only  unknown  variable  in  relating  pd,  pja,  and  rt.  Equa¬ 
tions  (14),  (17),  and  (18)  can  be  combined  to  obtain  an  explicit  expression  for  min¬ 
imum  SNR  as  follows: 

SNRmtn  =  ^±-  1  (19) 

In 

From  Equation  (19),  pd  and  pja  are  seen  to  be  related  via  SNR  to  measurable 
system  parameters  —  received  signal  and  noise  powers.  For  multi-pulse  detection, 
significant  improvement  in  SNR  can  be  achieved  through  pulse  integration,  where 
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the  magnitude  of  improvement  is  dictated  by  target  characteristics  and  pulse  inte¬ 
gration  techniques.  Calculation  of  p 4  and  p/a  in  this  event  is  complicated  since  the 
distribution  of  noise  and  signal  powers  at  the  integrator  output  is  required  (73:2-18). 
For  a  non-fluctuating  received  signal  and  post-detection  integration.  Figure  16  shows 
the  minimum  SNR  requirement  for  pj  =  0.95  and  three  values  of  pja  (73:2-21). 


Figure  16.  Minimum  SNR  Requirement  using  Post-Detection  Multi-Pulse  Integra¬ 
tion,  and  pi  =  0.95  for  various  pja  (Extracted  from  Skolnik  (73:2-21)) 


A  relation  between  SNR  and  system  parameters  is  required  so  that  design 
variables  can  be  adjusted  to  achieve  a  desired  p^  and  p/0.  Assuming  that  the  noise 
contributions  are  uncorrelated,  SNR  at  the  detector  can  be  defined  in  terms  of 
generated  photocurrent  (19:2399): 


SNR 


(^)2 

BO2 


(20) 


where  (i,)  is  the  mean  photocurrent  generated  by  the  detected  IF  signal,  and  (in)  is 
the  mean  photocurren*  generated  by  the  nlh  noise  process. 


In  a  shot-noise  limited  system,  SNR  reduces  to 


SNR 


(21) 


where  (iq)  is  the  mean  photocurrent  generated  by  LO-induced  shot  noise. 


The  shot  noise  current  term  of  Equation  (21)  can  be  determined  by  considering 
the  amount  of  LO  power  converted  to  noise  photocurrent  over  the  detector  and  IF 
passband.  Therefore,  (i7)2  may  be  written  as  (8:111,  78:87-37) 


(Plo)Tlq  '\Bq2 
he 


(22) 


where 

(Pi0)  is  the  mean  LO  power, 

r/q  is  the  detector  quantum  efficiency, 

A  is  the  photon  wavelength, 

B  is  the  detector  and  IF  bandwidth, 
q  is  the  charge  per  electron, 
h  is  Planck’s  constant,  and 
c  is  the  speed  of  light. 


The  signal  photocurrent  term  of  Equation  (21)  can  be  determined  by  using 
the  IF  component  of  the  photo-mixed  received  and  LO  signals.  Thus,  (z,)2  may  be 
written  as  (19:2399;  73:37-38) 

(h)2  =  (P,;)!  )'  (23) 
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or 


(i,)2  =  (pto)(pr)  r-^ 


2 


(24) 


where  (P,/)  and  (Pr)  are  the  mean  IF  power  and  mean  received  signal  powers,  re¬ 
spectively. 

Therefore,  using  Equations  (21),  (22),  and  (24),  SNR  may  be  expressed  by 


SNR 


(Pr)  Xrjq 
hcB 


l-’Jj 


Equation  (25)  implies  that,  for  selected  detector  hardware  and  wavelength  in  a 
shot-noise  limited  system,  the  SNR  value  is  d  u-  rmined  solely  by  received  power.  In 
this  situation,  p j  and  ptz  now  depend  directly  on  (Pr).  The  value  of  ( PT )  is  commonly 
expressed  in  terms  of  transmitter  and  receiver  parameters,  target  characteristics,  and 
transmission  properties  using  the  LADAR  Equation  (73:1-4;  37-4). 


The  LADAR  Equation 

Several  forms  of  the  LADAR  Equation  are  used  in  the  literature  (41;  88;  8: 
6),  their  primary  differences  being  related  to  detection  of  either  a  ‘resolved’  target 
(transmitted  beamwidth  smaller  than  the  target  dimension),  or  ‘unresolved’  target 
(beamwidth  larger  than  the  target).  Since  the  extremely  long  ranges  considered  in 
this  analysis  mean  that  targets  will  be  unresolved,  the  following  form  of  the  LADAR 
Equation  from  Lange  and  Evans  will  be  used  (48:156): 


<^>  = 


Pi  Gt  Tt o  tA 

4t r/?2 


G l  C 


rA  X2  Gr  Tro 
(47 tRY 


(26) 


where 


Pt  is  the  transmitter  power, 

Gt  is  the  transmitter  power  gain, 
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Tto  is  the  transmitter  optics  transmittance, 
ta  is  the  one-way  atmospheric  transmittance. 

R  is  the  range  to  the  target, 

<7(CJ  is  the  target  LADAR  cross-section, 

Gr  is  the  receiver  power  gain,  and 
rro  is  the  receiver  optics  transmittance. 

Equation  (26)  has  been  expressed  for  co-located  transmitter  and  receiver,  but 
with  possibly  different  transmitter  and  receiver  optics.  The  first  term  in  Equa¬ 
tion  (26)  represents  the  amount  of  power  in  the  transmitted  antenna  pattern  at  the 
target’s  range.  The  last  term  represents  the  amount  of  power  reflected  from  the 
target  which  is  detected  by  the  receiver. 

Variations  in  (Pr)  will  occur  from  detection  to  detection  due  to  differences  in 
target  and  atmospheric  characteristics.  These  variations  will  appear  in  the  terms 
ta  ,  <7 ic3 ,  and  R.  Additional  variations  arise  in  the  antenna  gain  terms.  Gt  and  Gr , 
depending  on  the  location  of  the  target  within  the  antenna  beamwidth.  Selection 
of  a  suitable  transmitter  laser  and  optics  will  serve  to  define  the  ‘fixed’  terms  in  the 
LADAR  Equation  —  P£,  A,  r£o,  and  tto.  However,  expressions  for  transmitter  and 
receiver  antenna  power  gain  require  further  analysis. 

Transmitter  Antenna  Gain 

The  characteristics  of  the  antenna  which  explicitly  affect  LADAR  performance 
are  represented  by  the  gain  term  of  the  LADAR  Equation,  which  is  analogous  to 
the  gain  of  conventional  microwave  radars.  Gain  describes  the  ability  of  an  antenna 
to  concentrate  transmitter  energy  in  a  given  direction  and,  in  particular,  pov;er  gain 
also  accounts  for  all  losses  associated  with  the  antenna  (73:9-5).  Antenna  power 
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gain  may  then  be  written  as  (73:9-12) 


~  4t,1( 

G,  =  — 


where  ,4(  is  the  transmitter  antenna  primary  aperture  area,  and  Tjt  is  the  transmitter 
antenna  efficiency. 

The  term  4x.4(/A2  is  the  ideal  gain  term  formulated  by  Silver  for  a  uniformly 
irradiated,  unobscured  circular  aperture  with  no  losses  (71:177).  Since  ideal  gain  rep¬ 
resents  the  maximum  possible  gain,  qt  represents  all  losses  associated  with  truncation 
of  the  transmitted  laser  beam  by  the  antenna  aperture,  obscuration  by  secondary 
telescope  elements,  and  diffraction  effects  in  the  radiated  pattern  (43:2134-2135). 

Determination  of  r/t  requires  analysis  of  the  radiated  antenna  power  pattern. 
This  radiation  pattern  can  be  characterized  by  the  Poynting  vector.  5,  which  rep¬ 
resents  the  irradiance  on  any  surface  perpendicular  to  the  direction  of  electric  and 
magnetic  fields  (10:9).  This  irradiance  has  a  magnitude  which  is  proportional  to 
the  square  of  the  magnitude  of  the  scalar  field  derived  using  the  Fresnel-Kirchoff 
diffraction  formula  (73:9-7),  and  has  an  average  value  of  (50:465) 

(S)=[-ct0E2  (28) 


where  c0  is  the  permittivity  of  free-space,  and  E  is  the  peak  amplitude  of  the  electric 
field. 

For  this  analysis,  th<-  Poynting  vector  definition  may  be  used  in  the  develop¬ 
ment  of  an  expression  for  the  far-field  radiation  pattern  resulting  from  laser  beam 
transmission  through  an  optical  telescope.  To  assist  with  this  development,  consider 
an  aperture  irradiated  by  a  wave,  E{Eq),  generated  from  point  P0  on  the  telescope 
axis,  as  shown  in  Figure  17.  The  resultant  scalar  diffraction  field  at  point  P ,  L  (P), 
is  found  by  integration  over  the  aperture  area  of  all  area  elements,  d£  dq,  which  com- 
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prise  the  aperture  field.  i'(P)  is  given  by  the  Fresnel-KirchofF  diffraction  integral  as 
(10:382:  73:9-7) 

f:  o 

L'{P)  =  -  jE[c.rj)t-~[cos(nJ)-cos(A.snd^dr]  (29) 

AREA 

where 

k  =  2-/  A, 

s  is  the  distance  from  the  aperture  to  the  observation  point. 

(n.r)  is  the  angle  between  the  aperture  normal  and  the  r  direction,  and 
(n,  i)  is  the  angle  between  the  aperture  normal  and  the  j  direction. 


Figure  17.  LADAR  Transmitter  Aperture  Geometry 


Several  simplifications  to  Equation  (29)  can  be  made  when  considering  LADAR 
applications.  First,  the  far-field  regime  of  an  antenna  pattern  is  conventionally  de¬ 
fined  as  starting  from  a  range  of  2D2/A  to  infinity,  where  D  is  the  primary  aper- 
ture  diameter  (73:9-9).  For  the  range  of  values  of  D  and  A  under  consideration, 
the  far  field  region  extends  outward  from  a  range  of  between  20  and  250  kilometers. 
These  values  allow  confinement  of  the  analysis  to  far-field  radiated  power,  so  that 
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(n,s)  ~  rr.  Additionally,  since  LADAR  telescopes  concentrate  energy  along  the  c 
axis,  (rc.r)  %  0.  Therefore.  Equation  (29)  reduces  to 

U{P)  =  ~\S  (30) 

area 

Furthermore,  the  distances  5  and  s'  can  be  written  in  terms  of  Cartesian  coor¬ 
dinates  as 

s  =  [(*  -  O1  +  (y  -  n)‘  + 
s'  =  (*2  +  r  +  =-)m 

Hence 

s  =  (V2  +  *12  +  C  ~  2x£  -  2 yrj}  '  (32) 

In  the  far-held  regime,  th<-  bincm'.C  *.*.*. pension  can  be  used  io  obtain  the  fol¬ 
lowing  approximation  for  s  and  i/s  (13:20,10:383): 

a  %  s'  —  (x£  +  yrj)/s' 

1  /s  ~  i/s' 

where  the  approximation  for  s  is  used  in  the  exponential  term  and  l/s  is  used  for 
the  amplitude  term.  Note  that  s  cannot  be  approximated  by  s'  in  the  exponential 
transverse-phase  term  since  the  effect  of  neglecting  the  smaller  terms  is  amplified 
through  multiplication  and  the  exponential  power  effect  (13:20).  Using  these  ap¬ 
proximations,  Equation  (30)  may  be  written  as 

U[P)  =  —  J  J £(£.»/)  PXP  (xs  +  Vl)  dV  (34) 

AREA  L  J 

As  outlined  in  Chapter  III.  the  candidate  system  under  consideration  includes 
a  secondary  telescope  element  which  partly  omscures  the  transmitted  and  received 
signal.  An  example  of  such  a  design  is  illustrated  in  figure  4  in  Chapter  III,  showing 
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a  typical  Cassegrainian  arrangement.  Since  this  analysis  assumes  diffraction  at  an 
obscured  circular  aperture,  polar  coordinates  are  appropriate  in  representing  points 
in  the  radiation  pattern.  The  following  substitutions  will  be  used  for  the  points 
(p.  0)  in  the  aperture,  and  (p' ,  0')  in  the  observation  plane: 


£  =  p  cos  0  77  =  p  sin  v 

x  =  p'  cos  0'  y  —  p'  sin  0' 


(35) 


Using  Equation  (34),  and  replacing  s'  by  p'/smO  in  the  argument  of  the 
integral  s  exponential  function,  the  scalar  field  becomes 


U(P)  =  /2T  f  E(p,  0)  sin  $  cosfp  -  v')  pdpdv  (36) 

As'  do  di 


where  a  is  the  primary  aperture  radius,  b  is  the  secondary  aperture  radius,  and  the 
difference-of-cosines  identity  has  been  used  to  simplify  the  exponential  argument. 

The  form  of  Equation  (36)  can  be  further  simplified  by  using  the  Bessel  function 
integral  relation  (10:395): 

]—  f2*  ejx  cos  B  ejn.J  dg  _  jn(x)  (37) 

2  7T  do 


where  the  Bessel  function  of  the  first  kind  of  order  zero,  J0(x),  is  required. 

Additionally,  if  the  incident  field  amplitude  is  circular  symmetric  and.  there¬ 
fore,  independent  of  0,  Equation  (36)  becomes 


U(P) 


jkeJk 


/  rn 

- - J I  E{p)  JQ(kp  sin  0)  p  dp 


This  expression  represents  the  scalar  field  at  a  point  in  the  far-field  of  an 
obscured  circular  aperture.  The  incident  field.  E(p),  can  be  described  by  a  Gaus¬ 
sian  amplitude  function  corresponding  to  the  fundamental  TEMoo  laser  mode.  The 


following  representation  of  E(p)  is  taken  from  Klein  and  Degnan  and  is  corrected 
to  include  the  free-space  impedance  term.  l/(ce0).  in  the  normalization  constant 


-j  is  the  beam  radius,  measured  from  the  axis  to  the  l/e2  power  point,  and 
R.y  is  the  phase-front  curvature  at  the  aperture. 

The  normalization  constant  is  chosen  so  that  the  incident  power.  Pinc.  over 
the  entire  aperture  plane  is  unity.  This  may  be  expressed  mathematically  using 
the  following  expression  —  also  taken  from  Klein  and  Degnan  and  corrected  for 
free-space  impedance  (43:2135): 

P,nc  =  —  | E(p)\2  pdpdv  =  1  (40) 

t.’sirig  Equation  (39),  the  propagating  field  distribution  can  be  determined  at 
a  distance  s'  from  the  aperture  to  be 


E(p)  =  C0  exp 


(11) 


Combining  this  result  with  Equation  (38)  provides  the  field  distribution  at  a  distant 
,4>servat  ion  point,  P,  resulting  from  irradiation  of  an  obscured  circular  aperture  by 
a  ( iaussian  field.  This  distribution  is  given  by 


Based  on  the  Poynting  vector  definition  of  Equation  (28),  the  irradiance  distribution 
at  P  can  now  be  expressed  bv 


up)  =  ~  mp)\' 


The  antenna  power  gain  can  then  be  expressed  in  terms  of  irradiance  using 
Klein  and  Degnan’s  definition  as  follows  (43:2135): 


G,[s'.9)  = 


Us'.  9) 


where  I0  is  defined  as  the  irradiance  at  a  distance  s'  from  a  unit  power  isotropic 
radiator  (since  the  incident  field  has  been  normalized  fo  unit  power),  and  is  given 
bv  (13:2135) 


Combining  Equations  (42),  (43),  (44),  and  (45),  the  transmitter  antenna  gain 
can  be  written  as 


G.ls',  9) 


jkp2  /II 


exp  V  ^  )  exp  I - ~  +  Mkp  sin  9)  p  dp  |  (46) 


This  result  agrees  with  the  gain  expression  given  by  Klein  and  Degnan,  since 
the  free-space  impedance  terms  cancel.  Equation  (46)  can  be  further  simplified  for 
the  far-field  by  setting  {I/s'  +  1 //£*,)  to  zero,  since  phase-front  curvature  h's  infinite 
radius  in  the  far-field.  This  simplifies  Equation  (46)  to 


Jrt-V)  -  npr  |£ JM-v.'") 


wliere 


.V  =  ka  sin#  =  off-axis  parameter. 


.-If  =  -a2  =  area  of  primary  aperture, 

a  =  a/u,\ 

7  =  b/a  =  obscuration  ratio,  and 

U  =  p  jo  . 

Comparison  between  Equations  (27)  and  (47)  allows  transmitter  efficiency  to 
be  written  as 

rjt  =  2a2  I  f  e~a2u  J0(Xul/ 2)  du  (48) 

I  I 

Depending  on  the  laser  beam  radius  used,  different  values  of  gain  can  be 
achieved.  This  variable  can  be  eliminated  from  Equation  (48)  by  finding  the  a 
value  corresponding  to  maximum  on-axis  gain  as  a  function  of  7.  Setting  .Y  equal 
to  zero  and  performing  the  Equation  (47)  integration  to  determine  an  expression  for 
on-axis  gain  yields 

C,(0)  =  tv4i'2Qi(e-“I-e-a!T!)2  (49) 

Equation  (49)  can  be  differentiated  and  set  to  zero  to  find  an  optimum  expres¬ 
sion  for  a  in  terms  of  7.  Using  second  order  perturbation  theory,  Klein  and  Degnan 
obtained  the  following  solution,  accurate  to  within  ±1%  for  7  <  0.4  (43:2136): 

a  «  1.12  -  1.307“  +  2.1274  (50) 

Using  this  expression  for  a  in  Equation  (48)  allows  transmitter  efficiency  to 
be  determined  as  a  function  of  the  off-axis  parameter,  .Y,  for  a  specific  obscuration 
ratio.  Figure  18  shows  Equation  (48)  plotted  for  three  values  of  7.  Efficiency  is  seen 
to  be  degraded  significantly  for  values  of  ,Y  above  0.5.  For  example,  a  transmitter 
aperture  of  60  cm  diameter  with  20%  obscuration,  and  a  wavelength  of  5  /mi,  will 
have  a  radiated  power  pattern  which  is  theoretically  degraded  irom  boresight  gain 
by  about  6  dB  at  3  microradians  off-axis. 
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Figure  IS.  Transmitter  Efficiency  Variation  against  Off-Axis  Parameter 
Receiver  Antenna  Gain 

Using  the  reciprocity  theorem  of  Helmholtz  (10:381),  gain  for  a  given  antenna 
is  the  same  for  receiving  or  transmitting  operations,  provided  the  same  source  is 
used.  In  this  analysis,  however,  the  transmitter  source  function  follows  a  Gaussian 
amplitude  distribution,  whereas  the  received  signal  from  a  far-held  point  source 
consists  of  plane  waves  at  the  aperture.  Therefore,  assuming  aberration-free  optics 
with  perfect  transmission  (since  rtptics  transmittance  is  already  included  in  the  rt0 
and  rro  terms  of  the  LADAR  Equation),  the  only  losses  in  a  receiver  antenna  are 
those  due  to  obscuration  and  photon  detection  efficiency  at  the  focal  plane.  This 
may  be  expressed  ma  .tically  as 


G 


r 


4?r  Ar 


4  ir  Ar 

A2 


/  i 

l 1 


72)  TlH 


(51) 


where  Ar  is  the  receiver  antenna  primary  aperture  area,  T]r  is  the  receiver  antenna 
efficiency  factor,  and  rjn  is  the  heterodyne  detection  efficiency. 
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In  a  similar  manner  ro  transmitter  efficiency,  receiver  efficiency  can  be  analyzed 
by  considering  the  radiation  pattern  at  the  focal  plane.  Using  the  geometry  of 
Figure  19,  the  incident  signal  at  the  receiver  aperture  consists  of  infinite  plane  waves 
of  the  form 

E{z)  -  Er0e]kz  (52) 


where  z  is  the  distance  along  the  telescope  axis,  and  Era  is  the  peak  amplitude. 


Figure  19.  Detector  Geometry 


Equation  (38)  can  then  be  used  to  obtain  an  expression  for  the  field  incident 
on  the  focal  plane  as  follows: 

,i./r  P]k{s'  +  c)  n 

Er(p) - - 2 - ; - /  Jo(kp  sin  9)  pdp  (53) 

s  Jb 

where  p  is  the  radial  coordinate  in  the  aperture  plane. 

Integration  can  be  completed  using  the  following  recurrence  relation  for  Bessel 
functions  ( 10:395): 

{xa+1  ,/n+1(.r)}  =  xn+1  Jn(x)  (54) 
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Thus,  Equation  (53)  produces 


Er(p)  =  ~ 


jkEroej^(s  +  ~)  [a  J,  (ka  sin  9) 


k  sin  9 


b  Jx{kb  sin  9) 
k  sin  9 


where  J\  is  the  Bessel  function  of  the  first  kind  of  order  one. 

Using  the  nomenclature  of  Figure  19,  where  pd  is  the  radial  coordinate  in  the 
detector  plane  located  at  the  effective  focal  length,  from  the  aperture,  the  scalar 
field  magnitude  at  the  detector  can  be  written  as 


or,  in  Cartesian  coordinates: 


\P(  m  £r0  f  ,  \ka{x\  +  yj)l/ 2 


-  bj , 


f  kb(xj  + 


where  and  yd  are  the  Cartesian  coordinates  in  the  detector  plane. 

The  form  of  Equation  (57)  is  seen  to  be  that  of  an  Airy  pattern  modified  by  the 
diffraction  pattern  due  to  the  central  obscuration.  The  resultant  Airy  disk  radius, 
/?4,  is  defined  as  the  distance  from  the  center  of  the  Airy  pattern  to  the  first  null. 
That  is,  substituting  RA  for  {x\  +  yd)1^2  in  Equation  (57)  produces  the  condition 
of  zero  field  magnitude,  which  can  be  written  as 

Ji(<t)  -  7  J,(7cr)  =  0  (58) 


.  k  Ra  .  ,  c,  /«ff  F  , 

where  a  —  ^  ,  with  t,  =  —  =  antenna  r -number. 


Equation  (58)  relates  Airy  disk  radius,  antenna  F-number,  and  obscuration 
ratio,  which  is  expressed  graphically  in  Figure  20.  For  a  constant  F-number,  the  size 


of  the  central  Airy  disk  decreases  as  the  size  of  the  secondary  telescope  element  is 
enlarged. 


Figure  20.  Airy  Size  Parameter  Variation  with  Obscuration  Ratio 


The  incident  field  described  by  Equation  (57)  mixes  with  the  LO  field  at  the 
detector  face.  The  irradiance  on  the  detector  surface  is  given  by  the  superposition  of 
received  and  LO  fields,  and  can  be  expressed  using  the  time-varying  representation 
of  the  incident  fields  (19:2398)  as  follows: 

liPi.t)  =  y  ^{EApJe^1}  +  xi[Edp<)eJ(u‘t  ~  (59) 

where  denotes  the  real  component,  u;t  and  a?;  are  tL„-  approximately-equal 

transmitter  and  LO  angular  frequencies,  and  <t>  is  an  arbitrary  phase  difference  be¬ 
tween  the  two  signals. 

Note  that  Equation  (59)  is  a  simplistic  model  since  perfect  phase-front  align¬ 
ment  between  the  two  signals  is  assumed.  In  practice,  this  assumption  may  be 
difficult  to  realize,  particularly  at  short  wavelengths,  since  the  maximum  allowable 
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'lag  angle’  between  the  two  beams  for  efficient  coherent  detection  is  directly  propor¬ 
tional  to  wavelength  (82:3427).  Therefore,  the  assumption  of  perfect  alignment  will 
be  relaxed  in  a  later  section  so  that  misalignment  effects  can  be  examined  separately. 
With  this  assumption.  Equation  (59)  may  be  expressed,  after  some  algebra,  as 


I(PdJ)  = 


CC°  ^  ~ ^  +  Er(Pi)  Ei(pd)  cos  [(u;e  -  t  +  <?] 


,  EliPd)  ,n  „  ,  EfiPd)  rn,  t 
+  - — - —  cos(2u.’(  t)  H - - —  cos  [2  (a;/ 1  -  6)\ 


+  Er(p,i)  Et(pd)  COS  [(c 


i)t-0 1  }  (60) 


The  first  term  in  Equation  (60),  which  represents  the  amplitude  of  the  re¬ 
ceived  signal,  may  be  neglected  in  comparison  to  the  amplitude  of  the  LO  signal. 
The  last  three  terms  in  Equation  (60),  which  represent  signal  components  of  about 
2u/f,  can  also  be  neglected  since  they  are  outside  the  detector  bandpass.  Therefore, 
Equation  (60)  can  be  simplified  to 

Ilpd,  t)  =  ~  |  y-  +  Er(pd)  Ei(pd)  COS  [(to,  -  U>i)  t  +  6] |  (61) 


The  instantaneous  signal  power  corresponding  to  this  irradiance  can  be  ex¬ 
pressed  as  four  times  the  power  from  one  detector  quadrant.  That  is, 


[L  rl 

P(t)  -  4/  /  I(xd,yd,t)dxddyd 

Jo  Jo 


(62) 


where  the  detector  dimensions  are  2 L  x  2 L,  or 


P{t)  =  r.f0  I  f  \E,(xd,yd)\2  dxddyd 
Jo  JO 

nL 

\Er{xd,yd)  Ei{ xd,yd)\  cos  [(u;,  +u y)t  +  6\  dxddyd  (63) 

l 
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The  first  term  in  Equation  (63)  can  be  interpreted  as  the  mean  LO  power, 
while  the  second  term  is  the  instantaneous  IF  power.  That  is. 


P(t)  =  (Pt0)  +  Ptf(t) 


(64) 


The  power  equations  presented  thus  far  are  sufficient  to  determine  an  expres¬ 
sion  for  receiver  heterodyne  efficiency,  tjh.  Heterodyne  efficiency  can  be  defined  as 
the  actual  signal  detected,  divided  by  the  maximum  possible  signal.  The  maximum 
signal  possible  is  achieved  with  perfect  mixing  of  the  LO  beam  and  the  power  received 
at  the  aperture.  This  received  aperture  power  may  be  described  by  (19:2398) 


<P,)  =  ~#0(-*2)(  1-72) 


(66) 


If  the  actual  detected  signal  power  is  characterized  by  the  IF  power,  then  the 
heterodyne  efficiency  may  be  expressed  in  a  dimensionless  form  as  follows: 


UH  = 


(P,)2 

(Pr)(Pto) 


(66) 


Using  Equations  (63)  and  (64),  the  mean  square  IF  power  is  extracted  as 


(Pf)2  =  j 


2 cea  f  [  \Er(xd,yd)  Et(xd,ydll  dxddyd 
Jo  Jo 


(67) 


The  explicit  expression  for  Er(xd,yd)  may  now  be  substituted  from  Equa¬ 
tion  157),  ano  terms  re-arranged  to  obtain 


( Pi/i 2  =  (P) 


4  ct0 


tt(1  -  q2) 


L  fL  Et(xd,yd) 


!  j 


(*$  +  y})1/2  [ 
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k(xUylYr‘ 


2  Fs 


7jJyk{x]  +  lA) 


2U/2' 


2  F, 


dxddijd 


(68) 
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Heterodyne  efficiency  may  then  be  found  from  Equation  (66)  using  the  (Plo) 
term  from  Equations  (63)  and  (64),  and  the  term  ( Plf)2/(Pr )  from  Equation  (68). 
Using  a  generalized  form  for  the  LO  distribution,  the  expression  for  is  now 


1  fL  Ei{xd,yd) 


/  7 


l4  +  yl)'/2 


( 

k(xd  +  y2d)l!2 

-  7  Jx 

l 

2  Ft 

k{x-  +  yj) 
2  F, 


d/2 


tt(1  -  72) 


[L  rL 

I  I  F((xd,  V d)  dxd  d }Jd 
0  Jo 


dxd  dyd  | 

2  (69) 


Equation  (69)  is  similar  in  form  to  that  of  Degnan  and  Klein  (19:2399).  except 
that  a  circular  detector  and  polar  coordinates  led  to  a  simpler  expression  in  their 
analysis.  A  choice  of  normalized  variables  in  similar  vein  to  Degnan  and  Klein 
provides  a  simpler  form  of  Equation  (69)  for  further  analysis  and  plotting,  as  follows: 


VH 


r  r  E‘(x'y) 

Jo  Jo  (i2  +  y2)K2 


Jx 


{x2  +  P) 


-  i-AMi2  +  s2)172]) 


2 


dx  dy 


where 

x  =  kxd/{2F,), 
y  =  kyd/(2F,),  and 
r  =  L/Ra- 

Equation  (70)  shows  that  careful  selection  of  an  LO  distribution  will  optimize 
heterodyne  efficiency.  Maximum  efficiency  is  obtained  when  the  spatial  distribution 
of  the  LO  signal  exactly  matches  the  received  signal  (82:3422).  This  situation  could 
only  be  realized  for  the  case  under  study  by  enlarging  the  LO  beam  and  focusing  it 
through  a  centrally  obscuied  aperture  onto  the  detector  (19:2401). 
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Two  other  LO  distributions,  which  are  simpler  to  implement,  are  the  Gaussian 
and  uniform  profiles  (82,  19,  41).  The  Gaussian  LO  distribution  would  be  used  in 
LADARs  where  a  single  detector,  or  quadrant  tracking  detector,  is  used  in  the  focal 
plane.  The  LO  beam  irradiance  on  the  detector  would  then  be  adjusted  so  that  the 
LO  beam  diameter  is  of  the  same  order  of  magnitude  as  the  detector  width. 


The  uniform  LO  distribution  represents  the  case  where  either  the  LO  beam 
diameter  is  much  greater  than  the  single  detector  width,  or  where  a  mosaic  focal  plane 
array  is  used.  In  order  to  adequately  irradiate  the  entire  array,  the  LO  beam  diameter 
is  increased  to  cover  all  array  pixels,  thereby  causing  the  LO  distribution  over  each 
pixei  to  be  approximately  uniform.  For  the  FPA  case,  a  Gaussian  LO  distribution 
for  each  pixel  would  more  closely  match  the  received  Airy  pattern  than  a  uniform  LO 
beam,  but  would  be  difficult  to  implement.  The  difference  in  performance  between 
the  two  cases  will  be  examined  more  closely,  using  Equation  (70)  as  the  basis. 

Gaussian  LO  Distribution.  For  the  Gaussian  case,  the  LO  beam  can  be  de¬ 
scribed  in  a  similar  fashion  to  Equation  (39).  r  sing  an  arbitrary  peak  amplitude  of 
C'i .  the  LO  amplitude  is 

\Ei{x,y)\  =  C,e 


where  ajl  is  the  LO  beam  radius,  measured  from  the  axis  to  the  1/e2  power  point. 
Assuming  that  the  Gaussian  LO  beam  and  the  Airy  disk  share  the  same  center, 
may  then  be  written  as 


where  t  —  Ra/^'l- 
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Note  thdt  Tj //  only  depends  on  distribution  of  the  LO  signal,  and  is  independent  of 
LO  amplitude. 


The  solution  to  Equation  (72)  may  be  found  numerically  for  a  given  ratio 
of  detector  dimension  to  Airy  disk  radius,  r,  and  a  specific  obscuration  ratio,  7. 
Knowing  7  enables  the  determination  of  a  via  Equation  (58).  However,  in  calculating 
rjH  using  Equation  (72),  the  ratio  of  Airy  disk  radius  to  LO  beam  radius,  t,  must 
be  optimized  for  maximum  heterodyne  efficiency.  Therefore,  an  iterative  technique 
is  required  to  determine  the  optimum  value  of  t  for  each  combination  of  r  and  7.  so 
that  maximum  heterodyne  efficiency  can  then  be  calculated. 


Variation  in  the  optimum  t  value  is  shown  against  r  in  Figure  21  for  different 
obscuration  ratios.  As  r  increases,  corresponding  to  increasing  detector  size  for 
fixed  Airy  radius,  the  LO  beam  radius,  ljl,  must  decrease  for  optimum  heterodyne 
efficiency.  This  is  because  more  of  the  Gaussian  tail  is  being  captured  by  the  detector, 
which  increases  shot  noise,  but  not  the  signal  level. 
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Figure  21.  Optimum  Airy-to-LO  Ratio  versus  Detector-to- Airy  Ratio  for  Gaussian 
LO  Distribution 
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The  corresponding  curves  for  rjH  are  shown  in  Figure  22.  As  the  detector  size 
increases,  more  of  the  Airy  pattern  is  received,  which  increases  efficiency.  As  L/ RA 
increases  above  0.S35,  however,  the  additional  shot  noise  generated  from  capturing 
more  LO  irradiance  tends  to  dominate  any  gain  in  received  Airy  signal.  This  is 
due  to  the  more  rapid  fall-off  in  Airy  pattern  power  as  opposed  to  the  Gaussian  LO 
distribution. 
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Figure  22.  Heterodyne  Efficiency  versus  Detector- o  -  Airy  Ratio  for  Gaussian  LO 
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Uniform  LO  Distribution.  For  the  case  of  a  uniform  LO  distribution,  Equa¬ 
tion  (70)  becomes 


Oh  = 


J, 

) 


(x?+y? 


F7r 


dx i  dyi 


x(  1  —  ~fi)a2r2 


'3) 


figure  23  shows  t]h  plotted  against  r  for  specific  values  of  7.  The  optimum 
ratio  of  detector-to- Airy  size  is  seen  to  be  approximately  0.635.  This  value  can  be 
interpreted  similarly  to  the  Gaussian  case  as  a  trade-off  between  signal  Dower  and 
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LO-ir.duced  shot  noise.  As  the  detector  size  is  increased  from  r  =  0.3  ,  more 
Airy  pattern  irradiance  is  available  to  mix  with  LO  irradiance,  thereby  providing 
higher  efficiency.  However,  as  the  detector  size  is  increased  beyond  r  =  0.635  to 
capture  the  much  lower  irradiance  in  the  'wings'  of  the  central  Airy  disk,  additional 
shot  noise  generated  from  the  uniform  LO  irradiance  starts  to  degrade  performance. 
This  effect  is  more  pronounced  in  a  square  detector  when  compared  to  a  circular 
detector,  such  as  that  analyzed  by  Degnan  and  Klein  (which  vielded  an  optimum  r 
value  of  0.75).  due  to  the  growing  'bare'  corner  areas  on  the  square  detector  which 
receive  predominantly  LO  irradiance. 
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Figure  23.  Heterodyne  Efficiency  versus  Detector- to- Airy  Ratio,  for  Uniform  LO 

Comparison  of  Figures  22  and  23  shows  that,  for  optimum  detector  size  and 
obscuration  ratio,  use  of  a  uniform  LO  distribution  will  decrease  heterodyne  efficiency 
to  a  level  of  about  Sb%  of  the  efficiency  gained  from  using  a  Gaussian  LO.  This 
represents  a  loss  of  about  0.7  dB  in  receiver  gain. 


Phase.  Front  MisaPqnmr.nt .  The  preceding  efficiency  equations  developed  for 
both  Gaussian  and  uniform  LO  distributions  have  assumed  peifect  alignment  be- 


tween  received  beam  and  LO  beam  phase  fronts.  As  discussed  previously,  dose 


physical  alignment  of  optics  and  detector  is  difficult,  but  possible.  Ho.vCver.  mis¬ 
alignment  may  also  occur  in  a  normal  detection  event,  depending  on  target  location. 
Energy  reflected  from  targets  located  off  the  boresight  axis  may  still  be  successfullv 
detected,  due  to  the  width  of  the  receiver  antenna  beam  in  the  far-field.  This  energy 
will  arrive  back  at  the  receiver  aperture  as  a  plane  wave  inclined  at  an  angle  to  the 
aperture  plane.  Figure  24  shows  the  resultant  relationship  between  the  LO  beam 
incident  on  the  detector  surface,  and  the  incoming  signal  (82:3422).  The  wave  vec¬ 
tor.  k,  represents  the  received  beam  with  a  different  phase  to  the  LO  beam,  where 
the  LO  beam  is  parallel  to  the  c  axis.  When  the  returned  signal  is  received  from  a 
diffuse  target,  the  detection  efficiency  is  independent  of  the  angle  A  (82:3422). 


Figure  24.  LO- Received  Beam  Misalignment  Geometry  (Reproduced  from  Wang 
(82:3422)) 

The  misalignment  factor  may  be  included  in  Equation  (53)  for  the  field  incident 
on  the  focal  plane  as  follows: 


Er{p) 


jkEra  eJk(s'  +  z) 


J  e-ik  '  P  J0(kps\n9)  p  dp 


(74) 


where  eJk  '  P  represents  the  misalignment  between  the  two  beams,  and  p  is  the  radial 
coordinate  in  the  receiver  aperture.  The  mixing  process,  and  hence,  misalignment. 


can  be  considered  to  occur  in  the  receiver  aperture  due  to  the  antenna  theorem  for 
het^odyne  reception  (69:3293). 

I’sing  Equation  (67),  the  mean  square  IF  power  is  now 

[  f\-iakps\nS  r  \kP(xl  + ylY/2 

XU  J°l  2 T, 

For  a  generJized  distribution,  the  equation  for  heterodyne  efficiency  can  then 
be  written  as 


where  3  =  ka  sin<5. 

Using  a  uniform  LO  distribution,  and  expanding  the  exponential  misalignment 
term  with  Euler’s  Formula,  Equation  (76)  yields 

m  =  ^'(T-V)  (l  J  J  dxJ dyd l  cos^p)Jo[p(xl  +  ylYl2\  pdp 

+  Jo  Jo  dXd  dyd  Jy  Sin(^)'/o[/5(Id  +  ^)1/2]  prfpj  |  ("") 

Equation  (77)  can  be  evaluated  for  various  obscurations  using  numerical  tech¬ 
niques,  provided  values  of  r  and  a  are  specified.  The  optimum  value  of  r  =  0.635 
can  be  taken  from  Equation  (73)  (plotted  in  Figure  23)  for  a  uniform  LO  distribu¬ 
tion  and  no  misalignment.  The  a  value  is  set  from  the  specified  obscuration  using 
Equation  (58)  (plotted  in  Figure  20).  Using  these  values  of  r  and  cr,  Equation  (77)  is 
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plotted  in  Figure  25  against  the  misalignment  parameter.  3.  for  various  obscurations. 
A  general  decrease  in  heterodyne  efficiency  is  observed  as  misalignment  increases. 
Interestingly,  higher  obscuration  ratios  are  not  as  sensitive  to  misalignment,  possibly 
due  to  additional  diffraction  effects  from  the  secondary  element.  The  transmitter  ef¬ 
ficiency  factor,  rjt,  can  now  be  easily  obtained  from  Equations  (73)  and  (77)  by  using 
the  following  relation: 

m  =  (1  -  Y)vh  ("$) 
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Figure  25.  Heterodyne  Efficiency  versus  Misalignment  Parameter,  3  —  ka  sine 


Having  finally  developed  the  expressions  for  transmitter  and  receiver  gain  in 
terms  of  efficiency  factors,  all  LADAR  Equation  variables  are  now  defined  in  terms 
of  system  design  parameters.  However,  the  calculated  value  of  received  power,  (Pr), 
cannot  then  be  used  with  Equation  (25)  to  find  SNR ,  since  not  all  received  power 
is  actually  detected.  This  loss  of  power  is  due  to  finite  detector  size,  and  the  FPA 
detection  process. 
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FPA  Detection  Analysis 

As  already  described,  the  detected  signal  is  compared  to  a  decision  threshold 
which  is  set,  in  turn,  by  combining  a  noise  measurement  with  a  desired  p<a.  For  FPA 
detection,  this  process  is  repeated  for  every  pixel  in  the  array.  If  the  received  image 
falls  over  many  FPA  elements,  signal  power  from  several  pixels  could  conceivably 
cross  the  threshold.  Of  the  pixels  which  have  crossed  the  threshold,  the  element 
with  the  largest  signal  is  discriminated  so  that  the  target  location  is  confined  to  one 
pixel  field-of-view. 

This  method  contributes  to  the  high  angular  resolution  of  FPAs,  but  is  only 
practical  if  the  maximum  signal  fraction  falling  on  any  single  pixel  is  sufficient  to  re¬ 
liably  cross  the  threshold.  Calculation  of  the  maximum  signal  fraction  for  each  likely 
detection  event  is  required  to  determine  an  expected  value  of  ( Pr )  for  FPA  detection. 
Therefore,  determination  of  the  expected  image  location,  and  resultant  maximum 
signal  fractions  is  required  before  the  expected  mean  received  power,  E{(Pr)},  can 
be  derived. 

To  simplify  the  analysis  of  image  detection  on  an  FPA,  the  received  Airy  pat¬ 
tern  can  be  approximated  by  a  disk  of  uniform  irradiance.  The  radius  of  the  ap¬ 
proximating  disk  can  be  set  by  considering  the  distribution  of  Airy  image  power. 
From  Equation  (56),  the  expression  for  Airy  power  normalized  to  peak  power  may 
be  written  as 

|£(0)|J 

where  pA  is  the  Airy  image  radial  coordinate  at  the  detector. 

Figure  26  shows  normalized  Airy  disk  power  plotted  against  normalized  image 
radius,  p,\/ Ra  ,  for  various  obscurations.  The  point  at  which  the  power  falls  to  1  /e  of 
the  peak  value  is  a  common  definition  for  the  ‘edge’  of  a  beam  or  image  (59:14).  This 
l/e  rau;us,  denoted  by  u>A,  can  be  seen  from  Figure  26  to  correspond  to  an  image 
radius  of  approximately  50%  to  53%  ot  the  A i i >  radius,  depending  on  obscuration. 


Ji{apAl Ra)  -  c?Pa/Ra ) 

<7{pa!Ra  a) 
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Therefore,  the  disk  approximation  represents  a  smaller  amount  of  received  power 
than  Ra  since  the  Airy  disk  tail  is  trimmed  off. 


Figure  26.  Normalized  Airy  Disk  Power  versus  Normalized  Image  Radius 


The  actual  fraction  of  signal  power  contained  within  a  radius  of  u>a  can  be 
determined  through  consideration  of  the  integrated  Airy  pattern  power  as  a  function 
of  image  radius.  Figure  27  shows  the  normalized  Airy  power,  integrated  from  the 
central  peak  (lower  integration  limit)  outwards,  and  plotted  as  a  function  of  outer 
radius,  pA  (upper  integration  limit).  The  Airy  disk  radii  are  marked  for  each  obscu¬ 
ration,  and  represent  approximately  84%,  73%,  and  49%  of  the  total  image  power 
for  obscurations  of  0,  0.2,  and  C.4,  respectively.  Conversely,  Figure  27  also  shows 
that  the  u>a  radii  represent  only  58.4%,  51.5%,  and  35.2%  of  the  power,  respectively. 

An  additional  FPA  characteristic  which  can  cause  a  high  false-alarm  rate  or 
poor  detectior  capability  is  detector  responsivity,  as  described  in  Chapter  III.  Assum¬ 
ing  that  each  pixel  has  been  individually  calibrated  to  achieve  uniform  responsivity, 
and  that  no  ‘dead’  pixels  exist,  the  fraction  of  power  encircled  by  uja  will  now  be  used 
to  determine  the  expected  maximum  signal  power  for  FPA  signal  detection.  Recall- 
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Figure  27.  Normalized  Integrated  Airy  Power  as  a  Function  of  Outer  Image  Radius 


ing  the  optimum  detector-to-Airy  size  ratio  for  a  uniform  LO  of  r  =  Lj  RA  =  0.635  . 
an  optimum  ratio  of  uA  i ,o  L  can  be  determined  as  follows: 


UA  _  UA  RA 
L  Ra  L 


i jja  I 

fll' 0.635 
0.787  for  7  —  0 

1  0.804  for  7  =  0.2 

0.834  for  7  =  0.4 


(80) 


The  assumption  that  the  Airy  disk  can  be  characterized  by  a  circular  disk  of 
uniform  irradiance  with  radius  u>\  allows  the  ratio  uiAl  L  to  be  fixed  for  each  obscu¬ 
ration.  FPA  detection  analysis  may  now  proceed  by  considering  a  two-dimensional 
FPA  with  pixel  dimensions  2 L  x  2 L  and  a  gap  between  pixels  of  width  g.  Figure  28 
depicts  possible  image  locations  on  a  section  of  FPA,  where  the  ratio  uA/ L  has 
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been  fixed  for  an  obscuration  ratio  of  0.2,  and  the  gap  width  has  been  set  at  0.61. 
C’leariy,  the  image  can  overlap  a  maximum  of  four  pixels.  Therefore,  assuming  that 
the  image  falls  randomly  on  the  FPA,  the  expected  maximum  received  signal  is  given 
by  (53:5) 

E{Sn}  =  (SI) 

n—  1 


where  5n  is  the  maximum  power  possible  from  a  single  pixel  within  a  group  of  n  pixels 
which  contain  all  of  the  image  disk,  pn  is  the  probability  that  all  power  is  contained 
within  n  pixels,  and  £{5/?)  is  the  expected  value  of  the  maximum  received  signal. 


Figure  28.  Possible  Image  Locations  on  an  FPA  for  7  =  0.2  and  Ljg  =  1.67 


The  total  power  contained  in  the  disk  and  Airy  image,  which  both  vary  with 
obscuration  ratio,  will  be  denoted  by  Pd„k{l)  and  P,m{ 7),  respectively.  The  value  of 
Sn  when  all  of  Pdi,k(l)  is  received  by  n  pixels  can  then  be  calculated  geometrically, 
by  considering  the  maximum  fraction  of  an  image  disk  overlapping  a  pixel.  This 
maximum  fraction  depends  on  the  disk  location,  which  can  be  uniquely  specified  in 
terms  of  the  location  of  the  disk  center,  since  ^4/  L  is  fixed.  Figure  29  gives  the 
distribution  and  geometry  of  image  centers  in  terms  of  array  and  image  dimensions, 
and  is  valid  for  g  <  w.4. 

The  expected  location  of  the  image  center  is  the  centroid  of  each  area  of  image 
centers.  Therefore,  for  n  =  1,  the  expected  image  center  is  the  center  of  the  pixel, 
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Figure  29.  Distribution  of  Image  Center  Locations  for  n  overlapped  pixels 


and  Si  =  Pdi,k{l)  for  all  values  of  g  <  uA.  The  expressions  for  S2,  S3,  and  SA  are 
considerably  more  complicated  and  are  developed  in  Appendix  A.  Figure  30  displays 
the  expressions  graphically  as  Sn/ Pdt,k(~()  plotted  against  the  detector  dimension 
parameter,  l  /  g,  for  an  obscuration  ratio  of  0.2.  Graphs  for  other  obscurations  appear 
in  Appendix  A,  and  show  similar  general  trends.  Since  the  expressions  are  only  valid 
for  g  <  ^4,  the  graphs  are  only  valid  for  the  range  L/r  >  1.2  to  1.3,  depending  on 
7 .  Figure  30  shows  that  the  maximum  signal  fraction  increases  as  L/g  increases,  as 
expected,  and  that  less  signal  is  detected  as  the  image  overlays  more  pixels. 

The  corresponding  pn  values  can  also  be  obtained  geometrically  by  dividing 
the  area  of  all  possible  image  center  locations  shown  in  Figure  29  by  the  total  area 
of  a  detector  element.  In  this  case,  the  area  of  a  detector  element  consists  of  the 
pixel  area  plus  half  of  the  adjoining  gap  area,  because  ea^h  pixel  can  be  considered 
to  ‘own'  half  of  the  surrounding  gap.  From’ Figure  29,  px  is  therefore 


=  [2(L-ua  +  g)}2 
(2  L  +  g)2 


(82) 
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Figure  30.  Maximum  Signal  Fraction  when  Image  Overlays  n  Pixels  and  7  =  0.2 

The  remaining  pn  expressions  are  also  developed  in  Appendix  A,  accompanied 
by  graphs  corresponding  to  7  =  0  and  7  =  0.4.  Figure  31  is  plotted  for  an 
obscuration  ratio  of  0.2.  The  major  variation  with  L/g  involves  px  and  p4.  As  L/g 
increases,  corresponding  to  a  smaller  gap,  the  probability  of  the  disk  falling  on  only- 
one  pixel  decreases  sharply,  since  the  disk  is  more  likely  to  overlap  adjacent  elements. 
Since  a  shrinking  gap  also  sharply  reduces  the  dead  area  situated  at  the  ‘crossroads' 
of  four  pixels  (proportionally  with  g 2  ),  p4  increases  much  faster  than  p3,  and  at  the 
expense  of  p2. 

From  Equation  (81),  both  Sn  and  pn  are  needed  to  determine  the  expected 
maximum  received  signal,  E{Sr).  Expressing  E{Sr}  as  a  fraction  of  Ptm( 7), 
Figure  32  shows  PIS’/*}  as  a  function  of  L/g,  for  various  obscurations. 

The  resultant  overall  trend  for  signal  fraction  to  increase  as  L/g  increases  is  as 
expected.  The  ‘dip’  in  the  curves  at  values  of  Lj g  <  2  can  be  explained  by  the  effect 
of  a  sharply  decreasing  px  combined  with  a  signal  fraction,  „  .  =  1  .  Since  the 

*1 liufcv  W 

higher  order  signal  fractions  are  much  lower  in  this  region,  the  overall  signal  ‘dips' 
before  rising  again  as  low-er  order  signal  fractions  become  significant. 


Detector  Dimension  Parameter.  Lj  g 


Figure  31.  Probability  of  Image  Overlaying  n  Pixels  with  7  =  0.2 

As  obscuration  increases,  the  lower  expected  signal  can  be  attributed  to  a 
smaller  Airy  disk  containing  less  power.  As  an  example,  choosing  an  arbitrary  L/g 
value  of  2  and  an  obscuration  ratio  of  0.2,  the  expected  maximum  fraction  of  Pim(*  ) 
received  is  about  29.4%.  This  fraction  of  ( Pr )  can  then  be  used  to  calculate  SNR, 
and  hence,  p.i  and  pja. 

The  preceding  analysis  has  developed  the  theory  necessary  to  determine  the 
major  LADAR  design  parameters,  including  transmitter  gain,  receiver  gain,  expected 
received  power,  SNR ,  pj,  and  p/a.  These  equations  will  provide  an  estimate  of 
whether  detection  is  possible,  and  with  what  likelihood.  Other  degrading  effects  such 
as  mild  turbulence,  atmospheric  scintillation,  target  depolarization,  target  scintilla¬ 
tion,  weather,  variations  in  manufacturers’  specifications,  component  deterioration, 
optics  misalignment,  and  signal  processing  losses  are  difficult  to  quantify,  with  some 
factors  exhibiting  wide  variation  between  detection  events.  Common  design  practice 
leads  to  the  addition  of  a  SNR  system  margin  to  account  for  such  factors. 

Several  other  limitations  exist  which  can  cause  complete  loss  of  signal  under 
certain  conditions.  They  include  strong  atmospheric  turbulence,  telescope  pointing 
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Figure  32.  Expected  Maximum  Signal  Fraction  for  FPA  detection 

accuracy,  telescope  mount  stability,  and  high  power  beam  effects  such  as  theimal 
blooming.  These  factors  have  been  discussed  in  Chapters  III  and  IV.  and  must 
be  included  in  any  comprehensive  LADAR  design  analysis;  however,  many  of  the 
deleterious  effects  can  be  minimized  considerably  by  careful  selection  of  LADAR 
site,  wavelength,  and  .specification  of  high-performance  hardware. 

I.  A  DAP,  Accuracy 

Although  reliable  signal  detection  is  the  primary  LADAR  design  goal,  simply 
knowing  that  a  target  is  'out  there’  is  insufficient,  information  for  surveillance  pur¬ 
poses.  The  target's  position  must  be  identified  accurately  by  determining  its  spatial 
coordinates  and  velocity  so  that  an  orbital  element  set  can  be  calculated.  The  accu¬ 
racy  with  which  this  can  be  performed  is  characterized  by  the  detector  element  size 
and  the  transmitted  waveform. 

Since  a  target  can  be  positioned  only  to  within  one  detector  element  on  the 
FPA,  the  angular  field-of- view  of  the  pixel  plus  half  of  the  surrounding  yap  deter¬ 
mines  the  system’s  angular  accuracy.  Conversely,  the  transmitter  pulse  width  can 

DO 


be  thought  of  as  a  travelling  wave  of  length  cr.  where  c  is  the  speed  of  light,  and  r  is 
the  pulse  width.  Neglecting  the  effects  of  pulse  dispersion,  a  target  of  large  spatial 
extent  can  generate  signal  reflections  from  anywhere  within  the  distance  cr  and  can 
therefore  only  be  positioned  with  a  range  accuracy  of  cr/ 2. 

Additionally,  the  measurement  accuracy  of  the  target’s  radial  velocity  is  de¬ 
termined  by  the  magnitude  of  the  Doppler  shift  and  the  performance  of  the  post¬ 
detection  circuitry.  For  a  target  with  radial  velocity  component  vr.  the  Doppier 
shift  is  given  by  2c  /A.  The  measurement  of  lateral  velocity,  or  angular  rate,  is 
performed  by  successive  detections  —  assuming  no  target  maneuvering  —  with  a 
measurement  accuracy  determined  by  the  time  between  independent  detections,  and 
the  post-detection  circuit  performance. 

LADAR  Resolution 

An  additional  detection  problem  arises  when  two  or  more  targets  are  located 
closely  together  in  angle,  range  and  velocity.  The  LADAR.' s  ability  to  discriminate 
between  targets  is  described  by  the  LADAR’s  angular  resolution,  range  resolution, 
and  velocity  resolution. 

LADAR  resolution  may  be  specified  in  terms  of  design  parameters  and  is  nor¬ 
mally  fixed  as  a  result  of  trade-off  between  waveform  and  circuit  complexity,  and 
the  parameters  which  affect  SNR,  such  as  optics  diameter,  obscuration  ratio,  wave¬ 
length.  and  pulse  width.  As  for  system  accuracy  calculations,  the  parametr’c  equa¬ 
tions  for  resolution  will  be  presented  for  the  case  of  a  simple  pulse- modulated  sinusoid 
of  constant  frequency. 

Using  the  same  argument  as  for  range  accuracy,  range  resolution  is  limited  to 
cr/ 2.  For  a  transmitter  pulse  width  of  1  p s,  two  targets  can  be  resolved  in  range  as 
close  as  150  m  apart.  The  velocity  resolution  for  a  heterodyne  system  is  given  by 
Bachman  as  approximately  A/27’  ,  where  T  is  the  duration  of  the  entire  transmitted 
waveform  (6:31). 


Angular  Resolution  Angular  resolution  is  commonly  specified  using  the  Ray¬ 
leigh  definition,  which  states  that  two  point  images  are  resolved  if  the  center  of 
the  diffraction  pattern  generated  by  the  first  image  falls  at  the  first  null  of  the  sec¬ 
ond  image's  pattern  (6S:1S8).  Without  considering  other  limitations,  this  definition 
sets  the  lower  resolution  lim't,  and  systems  operating  at  this  limit  are  said  to  be 
diffraction-limited.  Adopting  the  Rayleigh  definition,  the  angular  resolution  for  a 
diffraction-limited  system,  Qdiff.  may  be  written  as 


Odifr 


Ra 

ft ft 


a 
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For  example,  an  obscuration  ratio  of  0.2,  wavelength  of  5  pm.  and  receiver 
aperture  primary  radius  of  30  cm  yields  an  angular  resolution  of  approximately 
9.7  microradians.  This  system  would  theoretically  be  able  to  discriminate  between 
two  targets  as  close  as  350  m  apart  at  geosynchronous  orbit. 

Two  other  limitations  to  angular  resolution  which  may  dominate  for  systems 
not  diffraction-limited  are  atmospheric  turbulence  and  detector  size.  The  scintil¬ 
lation  effects  of  turbulence  can  be  estimated  from  the  transverse  phase  coherence 
length  described  in  Chapter  IV.  For  example,  a  wavelength  of  10  pm  and  coherence 
length.  r0,  of  3.6  m  yields  a  resolution  of  2.4  microradians,  while  for  A  =  0.5  pm  , 
and  r0  =  .1  m,  the  resolution  limit  is  4.S  microradians  (87:6-29). 

The  effect  of  detector  size  is  more  severe,  since  each  Airy  pattern  must  not  only 
be  imaged  on  separate  pixels  for  successful  detection,  but  must  also  be  resolved  in 
a  separate  range  and/or  Doppler  cell  to  differentiate  between  a  single  Airy  pattern 
overlapping  two  pixels,  and  two  separate  patterns  on  adjacent  pixels.  Assuming 
both  conditions  are  successfully  discriminated,  then  the  system  resolution  can  be 
defined  in  terms  of  the  detector  instantaneous  field-of-view  (IFOV).  If  the  IFOV  of  a 
square  detector  element  is  defined  as  the  IFOV  of  a  pixel  plus  half  of  the  surrounding 
gap ,  then  the  resolution  limit  for  a  detector-limited  system,  Odeti  is  identical  to  the 
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system's  angular  accuracy,  and  is  given  by 


■2L+g 

/rff 


For  example,  a  pixel  size  of  50pm  x  50  pm  and  a  gap  of  12.5  pm,  combined 
with  a  3  m  focal  length,  yields  a  resolution  limit  and  angular  accuracy  of  only 
20.8  microradians. 

Clearly,  minor  variations  in  detector  or  telescope  parameters  may  cause  the 
angular  resolution  to  be  limited  by  either  diffraction  or  detector  factors,  with  turbu¬ 
lence  becoming  dominant  only  under  extremely  strong  turbulence  conditions. 


Summary 

Parametric  analysis  of  the  LADAR  detection  problem  has  produced  a  set  of 
equations  which  can  be  used  as  the  basis  for  a  first-level  theoretical  design  of  a 
candidate  ground-based  deep-space  LADAR.  The  effects  of  beam  truncation  and 
diffraction,  heterodyne  misalignment,  off-boresight  detection,  optics  transmittance, 
and  atmospheric  transmittance  are  included  in  the  analysis.  However,  the  effects  of 
atmospheric  turbulence  and  scintillation,  target  scintillation,  telescope  mount  sta¬ 
bility,  and  high  power  beam  effects  are  not  included. 

A  summary  of  the  parametric  relationships  is  now  outlined  for  the  case  of  \t 
transmitted  pulses,  obscured  circular  transmitter  and  receiver  apertures,  heterodyne 
receiver  configuration,  FPA  detection,  and  uniform  LO  distribution  over  the  FPA 
surface.  The  FPA  consists  of  a  two-dimensional  array  of  square  pixels,  with  gaps 
between  each  pixel. 

For  a  specified  mean  time  between  false  alarms,  tja,  the  probability  of  false 


larm,  p/a,  is 


Pfa  = 


L(o  (  1  -  <5off) 
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where  A/,  t3 ,  and  £0fr  are  chosen  by  the  designer  and  are  constrained  by  transmitter 
performance  and  receiver  complexity.  With  A/.  pja,  and  p,{  set.  Figure  16  is  thei  used 
to  find  the  minimum  required  SXR  to  meet  detection  specifications.  An  arbitrary 
system  margin  is  then  added  to  S X R  to  cover  miscellaneous  system  losses. 


Minimum  average  received  power  is  then  calculated  to  meet  S X R  requirements 


using 


(Pr)  = 


SXR  hcB 
A  nq 


I S6) 


where  the  wavelength  and  bandwidth  are  fixed  from  the  selection  of  laser  and  detec- 
tor/IF  hardware. 

With  minimum  required  (Pr)  determined,  the  LADAR  Equation  is  then  used 
in  an  iterative  process  to  determine  the  optimum  combination  of  transmitter,  tele¬ 
scope.  and  detector  characteristics  which  will  meet  performance  specifications,  with 
consideration  given  to  complexity,  cost,  and  technology  availability.  The  final  form 
of  the  LADAR  Equation  is 

/  n  \  E{ Pn}  Pt  At  Ar  T;o  Tro  7}  a,cs  ( l  -  72)  Vt  ’Ih 

17xr<  (8,) 


where  the  value  of  is  determined  from  Figure  32,  based  on  detector  dimen¬ 

sions  and  telescope  obscuration.  The  value  of  rjt  used  in  the  LADAR  Equation  is 
determined  using 
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which  is  plotted  in  Figure  IS. 


The  value  of  pn  used  in  the  LADAR  Equation  is 
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calculated  for  a  uniform  LO  distribution  with  no  misalignment  from 
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which  is  plotted  in  Figure  23.  Misalignment  effects  can  be  included  in  the  heterodyne 
efficiency  term,  and  are  determined  using 


On 


Jq  dxi  dyd  J  cos(3p)  J0[p(x‘d  +  yd)1/2  pdp 
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which  is  plotted  in  Figure  25  using  the  optimum  LjR,\  ratio  of  0.635. 

Based  on  a  specified  detection  probability,  the  preceding  equation  set  allows 
theoretical  analysis  of  detailed  LAD  AR  design  feasibility,  subject  to  the  assumptions 
already  outlined.  LADAR  performance  in  terms  of  accuracy  and  resolution  is  then 
calculated  using  the  following  equations: 


Range  Accuracy  = 

cr 

2 

Doppler  Shift  = 

2vr 

A 

Range  Resolution  = 

cr 

2 

Velocity  Resolution  = 

A 

2 f 

Angular  Resolution  = 

Ha  __  * 

Ait  ka 

;  diffraction-limited 

= 

2  L+g 

Jrff 

;  detector-limited 

The  parametric  analysis  summarized  by  this  set  of  design  equations  has  ex¬ 
amined  the  losses  associated  with  the  transmitting  and  receiving  telescopes,  and 
detection  efficiency  for  a  wide  range  of  detector  arrangements.  LO  configurations, 
and  beam  misalignment  conditions.  These  equations  allow  calculation  of  the  primary 
LADAR  performance  parameters  in  terms  of  the  LADAR  hardware  specifications, 
and  include  optimization  factors  for  detector  dimensions.  LO  beam  size,  and  tele¬ 
scope  dimensions. 
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VI.  Candidate  Design 


Introduction 

This  chapter  applies  the  LADAR  theory  presented  in  Chapter  V  to  the  prob¬ 
lem  of  deep-space  satellite  detection.  LADAR  component  specifications  are  based  on 
current  technology  as  outlined  in  Chapter  III,  using  the  two  candidate  laser  wave¬ 
lengths  drawn  from  Chapter  IV.  Calculation  results  are  presented  for  two  LADAR 
sdes,  using  both  wavelengths,  with  detection  under  ideal  and  worst-case  conditions. 
Calculations  are  performed  for  two  transmitter  waveforms:  an  unmodulated  pulse 
train,  and  a  linear-FM  chirp  pulse  train.  Two  alternative  design  concepts  are  then 
discussed  which  are  aimed  at  greatly  increasing  the  amount  of  received  power. 

Design  Calculations 

The  fundamental  question  in  this  analysis  is  whether  deep  space  objects  can 
be  detected  by  a  ground-based  LADAR  and,  if  so,  under  what  conditions.  Design 
feasibility  depends  not  only  on  successful  detection  of  deep-space  objects,  but  on 
reliable  detection  of  most  likely  targets  over  a  wide  range  of  viewing  geometries 
and  ‘seeing’  conditions.  If  overall  LADAR  performance  cannot  be  demonstrated  to 
be  significantly  better  than  existing  systems  such  as  GEODSS,  then  further  devel¬ 
opment  is  not  justified.  ‘Significantly  better’,  in  this  case,  means  that  substantial 
improvements  in  weather-tolerance  and  daylight  operation  must  be  achieved  while 
maintaining  a  detection  capability  comparable  to  GEODSS. 

The  diversity  in  conditions  at  the  various  GEODSS  sites  is  an  indication  of 
the  range  of  weather  and  atmospheric  conditions  likely  to  affect  LADAR  operation. 
For  each  GEODSS  location,  two  extreme  detection  cases  are  considered  for  design 
feasibility  calculations.  The  ‘best’  detection  case  is  defined  as  LADAR  detection  of 
a  geosynchronous  target  at  a  36,000  km  range  which  is  positioned  on  the  telescope 
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boresight  and  located  at  the  sensor  zenith.  Other  best-case  conditions  include  clear 
skies  and  perfect  heterodyne  alignment. 

Similarly,  the  ‘worst’  case  is  defined  as  detection  of  a  geosynchronous  target 
located  off-boresight  at  a  look  angle  of  60°  and  range  of  40.000  km.  Since  a  typical 
beamwidth  for  a  long-range  LADAR  is  about  10  /mad  (77:153),  the  worst-case  off- 
boresight  angle  is  set  at  5  p  rad  to  coincide  with  the  beam  ‘edge'.  For  the  receive 
path,  the  corresponding  worst-case  heterodyne  misalignment  is  also  set  at  5  p rad. 
Operation  through  rain  is  not  considered  due  to  the  associated  severe  attenuation; 
however,  cirrus-cloud  attenuation  is  included  in  the  worst-case  scenario. 

Since  atmospheric  attenuation  varies  between  GEODSS  sites,  Socorro  and 
Diego  Garcia  are  used  as  two  examples  of  typical  seeing  conditions.  Therefore,  de¬ 
sign  calculations  are  required  for  four  combinations  of  detection  conditions  at  each 
of  the  two  candidate  wavelengths  selected  in  Chapter  IV. 

The  general  methodology  used  in  design  calculations  is  shown  in  Table  6  as  a 
sequence  of  design  steps  aimed  at  determining  two  values  for  comparison:  (1)  the 
required  minimum  received  power  for  successful  target  detection,  and  (2)  the  calcu¬ 
lated  received  power  from  the  candidate  design.  The  overall  methodology  is  based 
on  the  procedure  summarized,  in  Chapter  V,  which  is  constrained  by  target  char¬ 
acteristics,  atmospheric  effects,  and  available  technology  as  summarized  in  previous 
chapters.  In  particular,  the  maximum  size  of  telescope  apertures  is  limited  by  both 
turbulence  effects  and  available  technology.  Also,  transmitter  power  is  limited  by 
both  available  technology  and  atmospheric  effects,  and  is  restricted  to  operation  at 
11.17  pm  or  4.64  pm. 

In  addition  to  hardware  constraints,  input  specifications  are  required  for  pd 
and  p/a.  Equation  (85)  in  Chapter  V  may  be  used  to  calculate  a  value  of  pja 

o 

from  a  desired  tja  after  selecting  a  transmitter  waveform  and  receiver  range-gate 
characteristics. 


Table  6.  Candidate  LADAR  Design  Sequence 


Design  Step 

Relationship 

Parameters 

Performance  Requirements 

1.  Set  pd  and  tja 

Pdi  tfa. 

2.  Choose  Transmitter 

A,  Pt,  rep.  rate 

3.  Choose  Waveform 

A/  ,  T  ,  tg  , 

4.  Calculate  pja 

Eq  (85) 

Pf* 

5.  Determine  SNR 

Figure  16  or  (73:2-22) 

6.  Set  System  Margin 

7.  Calculate  Minimum  SNR 

SNR+margin 

8.  Define  Target  Characteristics 

k*r ,  &lc$i  Ry 

V>sy  <5 

9.  Calculate  Bandwidth 

Doppler +  chirp+linewidth 

B 

10.  Choose  Detector 

Vqi  L,  L  / Q 

11.  Calculate  Minimum  Power 

Eq  (86) 

Required  Pr 

LADAR  Equation  Calculations 

12.  Determine  Maximum  Signal  Fraction 

Figure  32 

E{SR] 

f’.m('l') 

13.  Estimate  Coherence  Lengths 

Figure  13 

^*0  ti  ^*o  r 

14.  Choose  Optics 

-Ar,  7"f0> 

Tro  i  7  >  /eff 

15.  Optimize  L/ Ra  for  max  t]h 

(repeat  from  Steps  10/14  if  needed) 

Figure  23 

L/Ra 

16.  Define  Weather  along  Signal  Path 

Figures  10,  11,  12 

ta 

17.  Calculate  Transmitter  Efficiency 

Figure  18 

Pt 

18.  Calculate  Heterodyne  Efficiency 

Figure  25 

Ph 

19.  Calculate  Design  Power 

Eq  (87) 

Design  Pr 

20.  Check  for  Power  Ratio  1 
(repeat  from  Step  1  if  needed) 

Design/Required  Pr 

21.  Calculate  Accuracy  and  Resolution 

Eq  (91) 

- - 1 
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Unmodulated  Pulse  Train 


For  the  first  candidate  design,  a  transmitter  waveform  consisting  of  SI  rect¬ 
angular  pulses  of  constant  frequency  is  considered.  Choosing  a  transmitter  pulse 
repetition  rate  of  300  pulses  per  second,  pulse  width  of  1  /zs,  and  range-gate  on-time. 
t3 ,  of  20  fis  implies  a  receiver  dead-time  fraction,  <5off,  of  0.994.  For  a  transmit¬ 
ted  pulse  train  of  M  =  110  pulses,  and  a  time  between  false-alarms  of  one  hour. 
Equation  (85)  yields  a  pja  value  of  10-4.  These  specifications  are  combined  with 
an  arbitrary  p 4  value  of  0.95  to  obtain  an  SNR  value  from  Figure  16  in  Chapter  V 
of  -3  dB. 

An  arbitrary  system  margin  of  3  dB  is  then  added  to  compensate  for  scin¬ 
tillation  and  miscellaneous  losses.  The  resultant  minimum  SNR  value  of  0  dB  is 
used  in  Equation  (85)  to  calculate  the  minimum  required  received  power.  Other 
input  parameters  required  are  IF  bandwidth  and  detector  quantum  efficiency.  Since 
geosynchronous  targets  will  have  only  small  radial  components  of  velocity,  bandwidth 
will  be  limited  to  100  MHz  to  cover  expected  Doppler  shifts.  Based  on  available  t^ch 
oology,  rjq  is  set  to  0.85.  Table  7  summarizes  the  design  specifications  assumed  for 
all  detection  cases.  Depending  on  which  wavelength  is  used,  the  minimum  required 
received  power  can  then  be  calculated.  This  value  must  be  met  or  exceeded  by  the 
value  of  expected  received  power  calculated  by  the  LADAR  Equation  in  Chapter  V. 

The  first  input  parameter  for  the  LADAR  Equation  is  the  expected  maximum 
signal  fraction  received  using  FPA  detection,  .  For  an  obscuration  ratio  of  0.2 

and  detector  dimension  ratio  of  4,  a  signal  fraction  of  0.32  is  obtained  from  Figure  32 
in  Chapter  V.  Secondly,  the  transmitter  power  must  be  defined.  Considering  the 
pulse  repetition  rate  and  pulse  width  selected,  a  laser  transmitter  capable  of  5  MW 
output  is  probably  the  limit  of  current  or  near-future  technology. 

Telescope  characteristics  must  now  be  specified  by  trading  off  beam  divergence 
and  aperture  area  with  loss  of  phase-coherence  and  available  technology.  Figure  13 
in  Chapter  IV  provides  rough  estimates  for  maximum  aperture  sizes  based  on  turbu- 
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Table  7.  Candidate  Design  Specifications  for  all  Detection  Cases  Using  an  Unmod¬ 
ulated  Transmitter  Pulse  Train 


Parameter  Value 

Parameter  Value 

Target  Cross-Section,  aiC3  4  m2 

Probability  of  Detection,  p^  0.95 

Probability  of  False-Alarm,  p/a  10-4 

Obscuration  Ratio,  7  0.2 

IF  Bandwidth,  B  100  MHz 

Detector  Dimension  Ratio,  Lj g  4 

Received  Signal  Fraction.  p^S'^  0-32 

Detector  Quantum  Efficiency,  0.85 

Pulse  Repetition  Rate  300  pps 

Number  of  Pulses,  M  110 

Range-Gate  On-Time,  tg  20  ps 

Receiver  Dead-Time,  0.994 

Pulse  Width,  r  1  ps 

Optics  Transmittance  0.85 

System  Margin  3  dB 

Minimum  SNR  Required  0  dB 

lence  effects,  but  current  technology  and  high  manufacturing  costs  suggest  a  practical 
upper  limit  for  aperture  size  of  about  1.5  -  2  meters  in  diameter.  Therefore,  max¬ 
imum  transmitter  aperture  size  is  limited  by  turbulence,  while  maximum  receiver 
aperture  size  is  limited  by  technology.  Apart  from  telescope  elements,  several  other 
elements  exist  in  the  optical  train  which  can  attenuate  the  LADAR  signal.  A  com¬ 
bined  one-way  transmittance  factor  of  0.85  for  the  optical  system  will  account  for 
these  losses. 

The  atmospheric  transmittance  value  used  in  the  LADAR  Equation  is  taken 
from  Figures  10,  11,  and  12  in  Chapter  IV,  depending  on  LADAR  site  and  cloud 
conditions.  The  remaining  input  parameters  —  transmitter  efficiency  and  heterodyne 
efficiency  —  are  derived  from  Figures  18  and  25  in  Chapter  V,  which  include  off- 
boresight  and  misalignment  effects.  Tables  8  and  9  summarize  the  viewing  conditions 
and  remaining  input  parameters  for  calculation  of  the  expected  received  power  at 
wavelengths  of  11.17  ^m  and  4.64  pm,  respectively. 

As  can  be  seen  from  the  tables,  heterodyne  and  transmitter  efficiencies  are 
more  sensitive  at  shorter  wavelengths  to  off-boresight  detection  and  misalignment. 
The  transmitter  power  at  4.64  pm  is  set  to  40%  of  the  13C02  output,  based  on 
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TabieS.  Received  Power  Calculations  at  11.17  //m  Using  an  Unmodulated  Pulse 
Train 


Diego  Garcia 

Socorro 

Best 

Worst 

Best 

Worst 

Sky  Conditions 

Clear 

Cirrus 

Clear 

Cirrus 

Look  Angle,  v. 

0° 

O 

O 

CO 

0° 

O 

O 

CO 

Target  Range,  R 

36,000  km 

40.000  km 

36,000  km 

40,000  km 

Boresight  Angle,  8 

0  /mad 

5  /mad 

0  /mad 

5  /mad 

Atmospheric  Transmittance.  rA 

0.58 

0.25 

0.98 

0.71 

Misalignment  Angle,  6 

0  /mad 

5  /mad 

0  /mad 

5  /mad 

Transmitter  Aperture  Diameter 

1  m 

1  m 

1.4  m 

1.4  m 

Receiver  Aperture  Diameter 

1.75  m 

1.75  m 

1.75  m 

1.75  m 

Transmitter  Power,  pt 

5  MW 

5  MW 

5  MW 

5  MW 

Transmitter  Efficiency,  r/( 

0.71 

0.45 

0.71 

0.38 

Heterodyne  Efficiency,  r/^ 

0.61 

0.5 

0.61 

0.47 

Design  Received  Power,  Pr  (W) 

5  x  10"16 

3  x  10~17 

3  X  10~15 

4  x  10~16 

Required  Received  Power  (W) 

2  x  10-12 

2  x  10~12 

2  x  10~12 

2  X  10"12 

Power  Ratio:  Design/Required 

2.2  x  lO"4 

1.4  x  10-5 

1.2  x  10~3 

9.8  x  10~* 

reported  conversion  efficiencies  for  frequency-doubling.  The  LADAR  Equation  was 
then  used  to  calculate  the  expected  received  power  from  each  design.  The  ratio  of 
'design’  received  power  to  ‘required’  received  power  must  equal  or  exceed  unity  for 
the  candidate  design  to  meet  the  specified  probability-of-detection  and  probability- 
of-false-alarm  criteria. 

Both  tables  show  that  reliable  LADAR  detection  of  geosynchronous  targets  us¬ 
ing  an  unmodulated  pulse  train  is  theoretically  infeasible,  even  for  best-case  detection 
conditions.  Calculated  values  of  received  power  at  11.17  /;m  are  about  four  orders  of 
magnitude  smaller  than  the  required  minimum,  and  about  five  orders  of  magnitude 
below  minimum  at  4.64  /im.  Since  efficiency  and  transmittance  factors  for  best-case 
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Table  9.  Received  Power  Calculations  at  4.64  //m  Using  an  Unmodulated  Pulse 
Train 


j 

Diego  Garcia 

Socorro 

Best 

Worst 

Best 

Worst 

Sky  Conditions 

Clear 

i-  1  ■  '■  - 

Cirrus 

Clear 

Cirrus 

Look  Angle,  v. 

0° 

O 

O 

tO 

0° 

60° 

Target  Range,  R 

36.000  km 

40.000  km 

36.000  km 

40.000  km 

Boresight  Angle.  9 

0  /mad 

5  /rrad 

0  /mad 

5  //rad 

Atmospheric  Transmittance, 

O.S 

0.43 

0.96 

0.68 

Misalignment  Angle,  S 

0  /mad 

5  /mad 

0  //rad 

5  ura.  d 

Transmitter  Aperture  Diameter 

0.3  m 

0.3  m 

0.4  m 

0.4  m 

Receiver  Aperture  Diameter 

1  m 

1  m 

i 

1.2  m 

1.2  m 

Transmitter  Power.  pt 

2  MW 

2  MW  ! 

i 

2  MW 

2  MW 

Transmitter  Efficiency.  rjt  j 

0.71 

0.53  ! 

0.71 

0.47 

Heterodyne  Efficiency,  rjfj 

0.61 

0.35 

0.61 

0.25 

Design  Received  Power.  Pr  (W) 

2  x  lO-16 

3  x  10-17 

1  X  10~15 

1  x  10"1* 

Required  Received  Power  (W) 

ox  10“ 12 

5  x  10~12 

5  x  10~12 

5  x  10~12 

Power  Ratio:  Design/Required 

4.S  x  10-5 

" 

5.3  x  10-* 

2.1  x  lO"4 

1.9  x  10-’ 

conditions  are  close  to  the  maximum  possible,  large  variations  in  received  power  can 
only  be  achieved  by  altering  transmitter  power  or  aperture  sizes.  However,  scaling 
transmitter  power  up  by  four  orders  of  magnitude  for  LADAR  applications  cannot 
be  done  with  current  technology,  and  is  impractical  due  to  possible  optics  damage 
and  air  breakdown.  Also,  the  aperture  sizes  used  in  design  calculations  are  close  to 
the  maximum  practical  for  acceptable  coherence  losses  and  manufacturing  cost. 

Some  improvement  may  be  gained  by  relaxing  the  probability-of-detection  cri¬ 
terion;  however,  the  huge  increase  in  received  power  required  would  result  in  a  value 
of  pi  which  is  too  low  for  practical  sensor  operation. 
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Alternative  Designs 

Examination  of  the  LADAR  Equation  shows  that  three  other  major  parameters 
affect  received  power:  wavelength,  range,  and  target  cross-section.  Extensive  anal¬ 
ysis  has  already  shown  that  the  two  selected  wavelengths  are  the  optimum  choices 
available  for  LADAR  applications.  Therefore,  range  and  target  cross-section  are  the 
only  remaining  design  variables.  By  altering  these  variables,  two  methods  of  im¬ 
proving  detection  performance  are  possible:  range  reduction  by  using  a  space-based 
sensor,  and  enhancement,  of  target  cross-section  by  using  a  corner-cube  reflector. 

Space-Based  LADAR.  free  space  loss  over  the  two-way  signal  path  is  ex¬ 
pressed  by  the  RA  term  in  the  LADAR  Equation,  and  represents  an  extremely  large 
power  loss  for  the  range  to  geosynchronous  orbit.  Reduction  of  range  can  only  be 
accomplished  through  the  use  of  a  space-based  LADAR  sensor,  which  has  added 
advantages  such  as  removal  of  turbulence,  aerosol,  and  weather  effects  on  signal 
propagation. 

For  space-based  detection  under  the  same  best-case  detection  conditions  as  in 
the  ground-based  case,  a  minimum-altitude  sensor  orbit  can  be  determined  which  will 
provide  the  desired  minimum  received  power.  Setting  both  transmitter  and  receiver 
aperture  diameters  to  a  more  conservative  value  of  1.4  m.  a  power  ratio  exceeding 
unity  is  achieved  at  both  wavelengths  with  a  range  of  6.000  km.  This  corresponds 
to  a  sensor  orbit  of  30,000  km,  which  means  that  one  sensor  could  be  used  for 
surveillance  of  the  entire  geosynchronous  ring  due  to  a  shorter  orbital  period.  Since 
both  apertures  are  the  same  size,  a  single  telescope  could  he  used,  with  automatic 
detection  carried  out  according  to,  sav,  a  search/t.rack  schedule  which  is  updated 
daily  from  a  single  ground-control  station.  Also,  space-based  LADAR  transmitter 
candidates  could  include  the  'TO>  laser  at  10.6  /nn  since  transmittance  is  no  longer 
a  selection  criterion.  Further  reductions  in  transmitter  power,  and  relaxation  of 
alignment  and  aperture-size  requirements  would  he  possible  at  higher  sensor  orbits. 
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Several  major  hurdles  must  be  o'  ercome  before  such  a  scheme  can  be  imple¬ 
mented.  Sophisticated  search  and  track  techniques  would  be  required  for  autonomous 
operation  of  the  sensor  platform.  A  large-capacity  power  supply,  efficient  thermal 
control  system,  and  stable  platform  are  also  required  for  sensor  support.  Probable 
the  two  most  critical  requirements  for  a  space-ba,sed  sensor  platform  of  this  type  are 
pointing  accuracy  and  platform  stability.  Some  method  of  relating  telescope  look 
angle  to  an  inertial  or  earth-centered  reference  frame  is  required,  along  with  isola¬ 
tion  of  the  telescope  mount  from  equipment  vibration  and  plattorm  mo'-ernent.  so 
that  pointing  accuracies  on  the  order  of  a  few  microradians  are  maintained.  Further 
analysis  of  space-based  sensor  performance  is  beyond  the  scope  of  this  study. 


Target  Enhancement.  The  majority  of  geosynchronous  objects  are  unenhance  ’. 
with  relatively  small  LADAR  cross-sections.  However,  if  a  corner-cube  reflector  de¬ 
signed  to  reflect  IR  radiation  is  attached  to  the  object  and  faces  earthward,  target 
cross-section  can  be  increased  by  several  orders  of  magnitude.  Although  military 
satellites  would  not  normally  be  fitted  with  such  devices  for  obvious  reasons,  the 
majority  of  geosynchronous  objects  consist  of  communications  satellites  and  rocket 
stages  of  commercial  origin.  Therefore,  it  is  conceivable  that  future  commercial 
satellites  could  be  fitted  with  re  Pec  tors  t,o  assist  in  their  detection  and  tracking. 


The  LADAR  cross-section  of  a  corner-cube  reflector  is  expressed  bv  Skolnik  as 


( 74:47-55) 


cries 


3.T 


(92) 


where  d  is  the  length  of  an  edge  of  the  corner-cube  reflector. 

For  a  reflector  of  dimension  a  =  5  cm,  the  corresponding  cross-sections  are 
2  x  10’  rid  and  1  x  10'’  m2  .  for  wavelengths  of  11.17  /mi  and  4.64  pm,  respectively. 
Table  10  shows  the  resultant  received  power  for  all  detection  cases  and  both  candi¬ 
date  wavelengths.  The  power  ratio  values  show  that  target  enhancement  can  easily 
provide  the  required  received  power  in  most,  detection  cases  by  using  a  relatively 
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small  reflector.  Since  reflector  cross-section  scales  proportionally  with  dA .  a  small 
increase  in  reflector  size  allows  successful  detection  for  all  cases  and  wavelengths. 
With  larger  reflectors,  smaller  transmitters  and  apertures  can  be  used  to  improve 
cost  efficiency. 


Table  10.  Received  Power  for  Enhanced  Target  and  an  Unmodulated  Pulse  Train 


Diego  Garcia 

Socorro 

Best 

Worst 

Best 

Worst 

A  =  1 1.17/mi 

y 

Target  Cross-Section.  cqCJ  (m*) 

2.10  x  105 

2.10  x  10J 

2.10  x  105 

2.10  x  105 

Design  Received  Power.  Pr  (W) 

2  x  10-u 

2  x  10-12 

1  x  10-’° 

2  x  10-11 

Required  Received  Power  (W) 

2  x  10"12 

2  x  10"12 

2  x  10"12 

2  x  10~12 

Power  Ratio:  Design/Required 

11.62 

0.74 

65.05 

9.50 

A  =  4.64/im 

Target  Cross-Section,  <7ic,  (m2) 

Design  Received  Power,  Pr  (W) 

.  1 

1 

1 

Required  Received  Power  (W) 

a 

1 

Power  Ratio:  Design/Required 

14.51 

1.61 

63.98 

5.71 

FM-Chirp  Pulse  Train 

As  discussed  in  Chapter  III,  pulse  compression  using  an  FM-chirp  pulse  train 
provides  high  range-resolution  due  to  an  effective  reduction  in  the  pulse  width  by 
t  the  compression  ratio.  Another  useful  characteristic  of  pulse  compression  is 

that  peak  power  of  the  compressed  pulse  is  also  increased  by  r  Bp\j.  The  level  of 
power  amplification/pulse  compression  is  limited  by  the  practical  bandwidth  of  the 
FM  chirp  and  the  transmitter  puisne  width. 

Due  to  the  narrow  linewidth  of  CO2  laser  transitions,  FM  bandwidths  of  up  to 
600  MHz  appear  to  be  practical  (6:92).  Pulse  widths,  on  the  other  hand,  depend  on 
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laser  excitation  and  cavity  design  (66:63).  Moody  et  al.  report  a  13CCU  amplifier 
with  t  —  30  fis.  pulse  power  of  7  MW.  and  pulse  repetition  rate  of  10  Hz  (57:102). 
Based  on  this  level  of  technology,  I3C02  amplifier  characteristics  with  r  =  50  /is, 
pulse  power  of  6  MW.  and  pulse  repetition  rate  of  100  Hz  appear  feasible  now  or  in 
the  near-future. 

Using  these  specifications,  design  calculations  were  repeated  for  the  two  sites 
under  ideal  and  worst-case  conditions.  Table  11  shows  the  results  of  these  calcula¬ 
tions  for  Bf  \f  =  600  MHz  and  r  =  50  /is.  IF  bandwidth  was  set  at  650  MHz  to 
allow  for  expected  Doppler  shifts  in  the  reflected  signal,  and  transmitter  power  at 
4.64  /im  was  again  reduced  to  2  MW  due  to  conversion  efficiency.  Aperture  sizes 
and  all  other  specifications  are  identical  to  those  listed  in  Tables  S  and  9. 

Table  11.  Received  Power  Calculations  Using  a  Linear-FNl  Chirp  Pulse  Train 


Although  detection  using  pulse  compression  at  4.64  /im  is  still  unsuccessful, 
received  power  is  within  two  orders  of  magnitude  of  the  required  minimum.  At 
11.17  /im,  Table  11  shows  that  detection  is  successful  under  most  conditions,  with 
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insufficient  received  power  only  under  worst-case  conditions  at  Diego  Garcia.  Re¬ 
calling  that  worst-case  conditions  correspond  to  off-boresight  detection  of  a  target 
through  cirrus  at  60°  look-angle  with  significant  heterodyne  misalignment,  the  cal¬ 
culations  show  that  the  pulse  compression  method  leads  to  reliable  target  detection 
under  most  realistic  operational  conditions.  The  high  GEODSS  site  elevations  at 
Taegu  and  Maui  could  also  be  expected  to  exhibit  similar  atmospheric  transmittance 
and  turbulence  characteristics  to  Socorro,  implying  that  detection  performance  at 
these  sites  is  likely  to  be  similar  to  the  calculations  shown  for  Socorro. 

Resolution  and  Accuracy.  Using  Equation  (91)  from  Chapter  V,  the  resolu¬ 
tion  and  accuracy  of  the  successful  pulse  compression  LADAR  design  can  be  calcu¬ 
lated.  For  an  obscuration  ratio  of  0.2,  primary  aperture  diameter  of  1.75  m.  and 
A  =  11.17  pm,  the  diffraction-limited  angular  accuracy  and  resolution  is  7.4  /zrad. 
Conversely,  for  a  detector  dimension  ratio  L/ g  =  4  ,  L  =  25  p m,  and  an  effective 
receiver-telescope  focal-length  of  4  meters,  the  detector-limited  angular  accuracy  and 
resolution  is  14.1  ^rad. 

Assuming  average  turbulence  conditions,  the  angular  performance  of  this  de¬ 
sign  is  detector-limited.  That  is,  the  design  will  theoretically  be  able  to  differentiate 
between  two  targets  which  are  only  500  m  apart  at  geosynchronous  orbit,  and  be 
able  to  determine  their  angular  position  to  within  500  m.  If  the  FPA  consists  of 
128  x  128  elements,  the  linear  FOV  at  36,000  km  range  is  about  65  km.  This  small 
■footprint’  means  that  search  rates  will  be  slow,  necessitating  a  reasonably  accurate 
a  priori  estimate  of  target  orbital  elements. 

Since  pulse  compression  is  used,  the  effective  pulse  width  is  given  by  1  /  Bf\f, 
which  corresponds  to  1.7  ns.  This  extremely  short  pulse  width  provides  a  range 
resolution  and  accuracy  of  25  cm  —  more  than  enough  for  satellite  detection  re¬ 
quirements.  Also,  since  many  long  pulses  are  transmitted  at  a  low  repetition  rate, 
the  velocity  resolution  given  by  A/2 T  is  extremely  small  at  about  5  pm/s. 


However,  the  long  duration  of  the  transmitted  waveform  and  the  resolution  of 
Doppler-filter  circuits  would  probably  be  limiting  factors  to  velocity  and  range  resolu¬ 
tion.  During  the  1.1-second  period  in  which  the  transmitted  waveform  is  interacting 
with  the  target,  the  target  may  move  significantly,  causing  inaccurate  measurement 
of  instantaneous  Doppler-shift  and  range. 

Sensitivity  and  Coverage.  The  sensitivity  of  LADAR  detection  depends  on 
target  characteristics,  viewing  geometry,  atmospheric  conditions,  turbulence,  and 
site  elevation.  With  good  conditions,  target  cross-sections  considerably  smaller  than 
4  m2  could  theoretically  be  detected,  especially  from  high  elevations. 

The  calculation  of  SNR  and,  hence,  received  power,  is  based  on  the  underlying 
assumption  of  shot-noise  limited  detection.  This  condition  occurs  if  enough  LO 
power  is  provided  so  tha*:  the  background  signal  is  negligible  in  comparison  to  the 
LO-induced  shot  noise.  Background  power  may  be  calculated  using  Figure  7  in 
Chapter  V.  The  expected  worst-case  background  spectral  radiance  at  11.17  pm, 
assumed  constant  across  the  IF  bandwidth,  is  about  5  W  m_2sr-1pm_1-  For  the 
FPA  dimensions  used,  the  receiver  FOV  is  about  2.54  micro-steradians,  while  the 
receiver  aperture  area  is  2.4  m2  and  IF  bandwidth  is  650  MHz.  These  values  yield  an 
expected  background  signal  power  of  about  8  nW.  Therefore,  providing  LO  power 
above  about  1  pW  will  ensure  shot-noise  limited  operation. 

The  maximum  look  angle  used  in  worst-case  design  calculations  was  60°,  cor¬ 
responding  to  GEODSS  maximum  angular  coverage.  Apart  from  Diego  Garcia,  all 
other  sites  should  be  capable  of  geosynchronous  coverage  equal  to  or  better  than 
GEODSS.  The  Diego  Garcia  site  may  be  restricted  to  smaller  look  angles  under 
worst-case  conditions. 
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Summary 


Design  calculations  show  that  reliable  detection  of  unenhanced  geosynchronous 
objects  using  an  unmodulated  pulse  train  from  a  ground-based  IR  LADAR  is  cur¬ 
rently  not  possible  at  either  of  the  candidate  wavelengths.  State-of-the-art  LADAR 
performance  parameters  were  used  to  calculate  received  signal  power  under  both 
ideal  and  worst-case  detection  conditions  for  the  GEODSS  sites  at  Socorro  and  Diego 


,VoiSt-».asc  couuitiouJ  included  maximum  look -angle  'dewing 


cloud,  with  a  target  located  5  /mad  off-boresight  and  heterodyne  misalignment  of 
5  /mad.  These  calculations  show  that  an  increase  in  received  power  by  at  least  four 


orders  of  magnitude  is  necessary  for  reliable  detection  of  a  typical  geosynchronous 
object.  The  low  level  of  received  power  is  primarily  due  to  the  extremely  large 
amount  of  free-space  loss  experienced  by  a  LADAR  signal  reflected  from  a  range  of 


36.000  km. 


Two  alternative  design  concepts  using  an  unmodulated  pulse  train  were  pre¬ 
sented  which  are  aimed  at  increasing  received  power.  The  first  entails  location  of  the 
LADAR  sensor  in  a  deep-space  orbit  to  reduce  the  target  range.  Calculations  show 
that  successful  detection  is  possible  for  a  sensor  located  at  an  altitude  of  at  least 
30,000  km.  Although  pointing  accuracy  and  platform  stability  are  critical  to  the 
success  of  the  space-based  sensor,  a  comprehensive  analysis  of  space-based  satellite 
detection  is  beyond  the  scope  of  this  research. 

The  second  alternative  design  concept  assumes  that  the  target  is  fitted  with  a 
small  corner-cube  reflector  which  is  reflective  at  IR  wavelengths.  Calculations  show 
that  the  resultant  dramatic  increase  in  target  cross-section  allows  successful  target 
detection  under  nearly  all  conditions  for  both  candidate  wavelengths. 

Finally,  design  calculations  were  performed  for  an  unenhanced  target  using 
the  original  ground-based  IR  LADAR  with  pulse  compression.  The  results  show 
that  successful  detection  is  possible  under  most  operational  conditions  at  11.17  pm, 
while  received  power  at  4.64  pm  is  two  orders  of  magnitude  below  the  required 
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minimum.  Resolution  and  accuracy  of  the  11.17  /jm  system  is  very  good,  and  deep- 
space  LADAR  coverage  is  at  least  as  good  as  GEODSS. 


Ill 


VII.  Conclusions 


The  purpose  of  this  research  was  to  determine  the  minimum  performance  re¬ 
quirements  for  a  ground-based,  IR  LADAR  for  deep-space  satellite  detection,  and 
to  select  a  candidate  sensor  design  using  available  technology.  The  first  stage  of  re¬ 
search  included  a  review  of  typical  deep-space  target  characteristics  and  long-range 
IR  LADAR  techniques.  A  typical  target  cross-section  cf  4  m2  ar.H  range  of  36.000  - 
40,000  km  was  selected  for  design  calculations,  along  with  LADAR  design  architec¬ 
ture  based  on  heterodyne  FPA  detection,  a  MOPA  transmitter  configuration,  and  a 
Ritchey-Chretien  telescope  design. 

Two  transmitter  waveforms  were  considered  for  the  LADAR  application:  an 
unmodulated  rectangular  pulse  train,  and  a  linear-FM  ‘chirp’  pulse  train.  The  FM- 
chirp  pulse  train  was  emphasized  since  it  allows  use  of  pulse  compression  techniques 
in  the  receiver  to  compress  the  effective  pulse  width,  and  increase  the  effective  pulse 
amplitude,  by  several  orders  of  magnitude.  A  target  acquisition  scheme  using  sepa¬ 
rate  transmit  and  receive  telescopes  was  also  presented  which  greatly  reduces  acqui¬ 
sition  time  and  improves  system  sensitivity. 

Characteristics  of  several  high-power  pulsed  lasers  were  then  reviewed,  includ¬ 
ing  Nd:YAG,  DF,  12C02,  13C02,  and  frequency-doubled  I2C02.  Of  these,  the  C02- 
variant  lasers  were  selected  for  further  study,  based  on  their  technological  maturity, 
high  power,  narrow  linewidth,  and  small  chirp  characteristics. 

The  second  stage  of  research  involved  calculation  of  background  radiance,  at¬ 
mospheric  transmission,  and  turbulence  effects  for  typical  operational  conditions 
across  the  IR  band.  Of  the  candidate  C02  lasers,  the  13C02  laser  using  the  P ( 2 2 ) 
transition  at  11.17  fim  and  the  frequency-doubled  12C02  laser  using  the  R(‘20)  tran¬ 
sition  at  4.64  /im  were  selected  as  optimum  in  terms  of  transmission  under  a  wide 
range  of  atmospheric  conditions. 


Other  results  derived  were  that  average  cirrus  cloud  conditions  have  relatively 
minor  effects  on  transmission,  while  rain  attenuation  is  severe  at  all  IR  wavelengths. 
Two  turbulence  effects  —  scintillation  and  coherence  loss  —  were  highlighted  as  se¬ 
rious  for  LADAR  applications.  However,  addition  of  a  system  signal-to-noise  margin 
serves  to  overcome  scintillation  effects,  while  limiting  the  size  of  telescope  apertures 
avoids  coherence  loss  under  most  turbulence  conditions. 

The  third  stage  of  research  consisted  of  parametric  analysis  of  the  LADAR 
detection  problem  using  diffraction  theory  for  obscured  circular  telescope  apertures, 
a  Gaussian-profile  laser  beam,  heterodyne  detection,  and  a  two-dimensional  array  of 
square  detector  elements  separated  by  finite  gaps.  The  effects  of  Gaussian  and  uni¬ 
form  LO-irradiance  over  the  detector  were  included  in  the  analysis,  along  with  beam 
truncation  and  obscuration,  heterodyne  misalignment,  off-boresight  detection,  optics 
transmittance,  atmospheric  transmittance,  and  a  diffuse  target  cross-section.  Addi¬ 
tional  equations  were  presented  to  allow  calculation  of  LADAR  accuracy  and  resolu¬ 
tion  performance.  Although  the  parametric  analysis  was  carried  out  for  ground-based 
deep-space  IR  LADAR  applications,  the  resultant  equation  set  is  sufficiently  general 
to  be  used  for  most  LADAR  designs  using  heterodyne  detection. 

The  set  of  derived  parametric  equations  was  then  used  in  the  final  stage  of 
research  to  assess  the  performance  of  several  candidate  system  designs.  LADAR 
component  specifications  for  the  designs  were  based  on  current  technology  using  the 
two  final  candidate  laser  wavelengths,  with  calculations  performed  for  both  ideal  and 
worst-case  detection  conditions  from  two  typical  LADAR  sites.  Using  the  unmod¬ 
ulated  pulse  train  waveform,  the  analysis  shows  that  an  increase  in  received  power 
by  at  least  four  orders  of  magnitude  is  necessary  for  reliable  detection  of  a  typical 
geosynchronous  object.  The  major  reason  for  low  received  power  is  the  amount  of 
free-space  loss  incurred  over  the  long  range  to  deep-space.  Additional  limitations 
exist  for  transmitter  power  and  aperture  sizes  due  primarily  to  atmospheric  effects 
«u<~h  as  turbulence  and  thermal  blooming. 


113 


Two  alternative  design  concepts  using  an  unmodulated  pulse  train  were  pre¬ 
sented  which  are  aimed  at  increasing  received  power.  The  first  concept  entails  loca¬ 
tion  of  the  LADAR  sensor  in  space  at  an  altitude  of  at  least  30.000  km  to  achieve 
successful  detection;  however,  pointing  accuracy,  platform  stability,  sensor  control, 
and  power  generation  issues  were  not  addressed  in  the  analysis.  The  second  alter¬ 
native  concept  assumes  that  a  small  corner-cube  reflector  is  fitted  to  the  target, 
providing  enhanced  signal  returns  which  allow  successful  detection  under  all  condi¬ 
tions.  Although  both  alternative  concepts  appear  theoretically  feasible,  considerable 
problems  exist  with  their  practical  implementation. 

Finally,  design  calculations  were  performed  for  an  unenhanced  target  using  a 
ground-based  IR  LADAR  with  pulse  compression.  The  results  show  that  successful 
detection  is  possible  under  most  operational  conditions  at  11.17  pm,  while  received 
power  at  4.64  pm  is  two  orders  of  magnitude  below  the  required  minimum.  Reso¬ 
lution  and  accuracy  of  the  11.17  pm  system  is  very  good,  and  deep-space  LADAR 
coverage  is  at  least  as  good  as  GEODSS. 

Selection  of  heterodyne  detection  and  pulse  compression  techniques,  combined 
with  recent  improvements  in  C02  laser  and  FPA  technology  has  produced  a  ‘first- 
level’  LADAR  system  design  which  is  theoretically  feasible.  Compared  to  GEODSS 
performance,  a  network  of  LADAR  sites  positioned  at  the  same  locations  as  GEODSS 
promises  much-improved  availability  due  to  the  capability  of  operation  in  daylight 
and  cloudy-sky  conditions.  However,  the  narrow  beamwidth  and  relatively  long  sig¬ 
nal  round-trip  time  to  deep-space  means  that  a  reasonably  accurate  a  priori  estimate 
of  the  target  orbital  elements  is  required  so  that  search-time  is  reduced. 

Recommendations 

Additional  preliminary  analyses  and  design  calculations  are  required  before 
implementation  of  the  successful  candidate  design.  More  specific  data  concern¬ 
ing  atmospheric  seeing  conditions  at  prospective  LADAR  sites  is  needed  in  order 
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to  determine  maximum  aperture  sizes,  and  to  accurately  assess  the  likely  effect  of 
haze  and  turbulence  on  system  performance.  Assuming  co-location  with  GEODSS. 
much  of  this  data  may  be  available  from  previous  GEODSS  environmental  studies. 
Additionally,  a  more  comprehensive  analysis  of  acquisition,  tracking,  and  pointing 
performance  is  required  in  light  of  the  large  telescope  size  and  narrow  beamwidths 
required.  Finally,  further  development  of  I3C02  laser  technology  is  desirable  so 
that  pulse  repetition  rate  is  increased  from  the  currently-reported  levels,  while  high 
pulse  power  is  maintained. 
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Appendix  A.  FPA  Detection  Analysis 


The  first  objective  of  FPA  detection  analysis  is  to  determine  the  maximum 
possible  power,  Sn,  from  a  single  pixel  within  a  group  of  n  pixels  which  contain  all 
of  the  image  disk.  The  image  disk,  in  this  case,  is  a  disk  of  constant  irradiancc  with 
radius  u\i,  and  is  used  to  approximate  the  received  Airy  pattern.  Exact  values  of  u 
for  different  obscuration  ratios  can  be  found  in  Chapter  V'.  The  location  of  the  image 
disk  center  is  confined  to  a  finite  region  when  the  image  overlaps  n  pixels,  with  the 
expected  location  of  the  image  center  given  by  the  centroid  of  the  area  of  possible 
image  center  locations.  Therefore,  Sn  is  found  by  first  calculating  the  centroid  of  the 
area  of  image  centers,  and  the  evaluating  the  expected  fraction  of  image  disk  which 
overlaps  an  active  detector  pixel. 

The  second  objective  is  to  find  the  corresponding  probability,  pn,  that  all  image 
power  is  contained  within  n  pixels,  so  that  Sn  and  pn  can  then  be  combined  to  cal¬ 
culate  the  expected  value  of  the  maximum  received  signal,  E{Sr}  —  Tins i  Snpn  . 
pn  is  obtained  by  calculating  the  area  of  all  regions  which  could  contain  the  image 
center  for  n  overlapped  pixels,  and  then  dividing  by  the  total  area  of  a  detector 
element.  Using  Figure  29  in  Chapter  V,  values  of  Sn  and  pn  for  n  —  2,  3,  and  4 
will  now  be  developed.  Values  for  n  =  1  are  given  in  Chapter  V,  while  values  for 
n  >  4  are  not  required,  since  this  analysis  constrains  the  image  disk  size  to  values 
less  than  the  pixel  size. 

Two  Overlapped  Pixels 

Figure  33  provides  the  detector  geometry  for  the  case  of  two  pixels  overlapped 
by  the  received  image.  The  pixels  have  dimensions  of  2 L  x  2 L  ,  and  are  separated 
by  a  gap  of  width  g.  Only  one  of  the  four  areas  of  possible  image  centers  requires 
analysis  since  the  other  three  have  identical  shape  and  size,  and  are  symmetrically 
placed. 


Figure  33.  Image  Center  Geometry  for  the  Case  of  Two  Overlapped  Pixels 


The  area  of  image  centers,  ,42,  shown  in  Figure  33  can  be  written  as  the  area 
of  a  rectangle  plus  the  area  of  two  diverging  wedge-shaped  sections  as  follows: 


The  centroid  of  the  area  of  image  centers  is  now  required  so  that  the  expected 
location  of  the  image  can  be  fixed.  Clearly,  the  y  centroid  coordinate  is  y  =  0  .  The 
x  value  can  be  obtained  by  calculating  x  =  My/ A 2  ,  where  My  is  the  area  moment 
with  respect  to  the  y  axis.  My  can  be  written  as: 


XL 


=  2  J  x(L  —  lja  +  g)  dx  +  2  J  x  L  +  g  -  {u\  —  x2) 


dx 


-  — (£  +  9)  ~  {9  ~  ua)2{L  ~  +  9) - 3^  + 


2u>3  ,  2  (u2A-g2/i) 


(94) 
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Therefore,  the  centroid  coordinates  are  (x,y)  =  (~^S0)  .  The  corresponding 
expected  area  of  overlapped  pixel,  02,  when  the  image  is  centered  at  (x.y).  is: 


Since  the  image  disk  area  is  .  the  expected  fraction  of  total  disk  power 
can  be  expressed  as: 


d>2  02  1 

r.T 

-  -  arcsin 

(  -V/y  kl 

A  U 

2 

P-i<A  7J  ~^A  ~ 

ho,'^  ,4  2 ) . 

T  .420,'*) 

V  .42  )  _ 

1/2 


(96) 


In  a  similar  fashion  to  the  calculation  of  px  in  Chapter  V,  the  value  of  p2  can 
be  found  by  dividing  the  area  of  all  possible  image  center  locations  by  the  area  of  a 
detector  element.  In  this  case,  the  area  of  a  detector  element  consists  of  the  pixel 
area  plus  half  of  the  adjoining  gap  area.  That  is. 


P2 


E42 

(2  L+gy 


(97) 


Equations  (96)  and  (97)  contribute  to  the  calculation  of  E{Sr)  for  the  case  of 
two  overlapped  pixels.  The  corresponding  equations  for  the  cases  of  three  and  four 
overlapped  pixels  are  now  developed  in  a  similar  manner. 
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Three  Overlapped  Pixels 


Figure  34  depicts  the  image  center  geometry  for  the  case  of  three  overlapped 
pixels.  The  area  of  image  centers.  .43,  can  be  written  as  the  sum  of  areas  formed 
from  three  '-ertical  slices  as  follows: 
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After  some  algebra,  the  expression  for  A3  reduces  to: 


The  centroid  for  the  n  =  3  case  is  located  at  (.r ,  y)  =  ■  However. 

My  =  Mr  =  M3  due  to  the  region's  symmetrical  nature.  Using  three  vertical  slices 
again,  the  expression  for  .\/3  is  given  by: 


After  some  algebra,  the  expression  for  .\/3  reduces  to: 


u>“ 
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(99) 


The  corresponding  area  of  overlapped  pixel,  0 3,  for  the  image  centered  at 
(jjU  can  be  found  from  two  different  cases  of  detection  geometry  corresponding 
to  two  different  ranges  of  L/g  ratio.  Figure  35  shows  both  cases,  where  the  centroid 
is  located  either  on  the  pixel  itself  [g  —  ^  <  o),  or  in  the  gap  region  (<7  —  >  oj. 
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Figure  35.  The  Two  Possible  Image-Pixel  Geometries  for  the  Cases  of  Three  and 
Four  Overlapped  Pixels 


For  the  case  of  g  —  ^  <  0 , 03  is  expressed  as  the  sum  of  the  four  pixel  areas 
separated  by  the  x  and  y  axes  as  follows: 


which  yields: 
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while  for  (<7  —  ^  >  Cl)  ,  the  expression  for  0 3  is: 


The  expected  received  fraction  of  total  disk  power  is  then  found  in  a  similar 
fashion  to  S2  as: 

=  (102) 

where  O3  is  determined  from  Equations  (100)  or  (101),  depending  on  FPA  dimen¬ 
sions.  The  expression  for  p3  is  found  in  the  same  manner  as  for  p2,  with  four  areas 
of  possible  image  center  locations  given  by  Equation  (98).  The  p3  expression  is  then 
written  as: 

4/1 3 


P3  =  ~. 


(2L  +  gV 
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Four  Overlapped  Pixels 

Figure  36  shows  the  image  center  geometry  for  the  case  of  four  overlapped 
pixels.  Area  A 4  can  be  written  as  follows: 


(104) 
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which  reduces  to 
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Figure  36.  Image  Center  Geometry  for  the  Case  of  Four  Overlapped  Pixels 

As  for  the  n  =  3  case,  A fy  —  Mx  =  A/4  due  to  the  symmetrical  region, 

\fA  is  expressed  as: 


12  4  3 


The  corresponding  area  of  overlapped  pixel,  0A,  is  found  to  be  identical  in  form 
to  the  case  of  three  overlapped  pixels.  The  image-pixel  geometry  used  to  calculate 
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O4  is  given  by  Figure  35  and  the  relevant  expression  for  0 4  when  (jg  -  <  oj  is 

written  as 


04  = 


21  1/2 


—w2a  arcsin 


(107) 


and  for  {jj  —  >  o)  ,  the  expression  is 
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From  Equations  (102)  and  (103),  S4  and  p 4  are  then  determined  using: 


54  _  O4 

Pdukh)  'wil 


(109) 


and 


4  A, 

r'~  {2L+g)l 


(110) 


Summary 

Equations  have  been  developed  which  specify  the  maximum  power  possible,  Sn, 
from  a  single  pixel  within  a  group  of  n  pixels,  which  contain  all  of  the  image  disk. 
These  equations  are  shown  graphically  in  Figure  30  (Chapter  V)  for  an  obscuration 
ratio  of  0.2,  and  in  Figures  37  and  38  for  obscuration  ratios  of  0  and  0.4,  respectively. 
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Figure  37.  Maximum  Signal  Fraction  when  Image  Overlays  n  Pixels  and 


0 


An  additional  set  of  equations  has  been  developed  for  the  corresponding  prob¬ 
ability,  pn  ,  that  all  image  power  is  contained  within  n  pixels.  These  equations  are 


shown  graphically  in  Figure  31  (Chapter  V)  for  an  obscuration  ratio  of  0.2.  and  in 
Figures  39  and  40  for  obscuration  ratios  of  0  and  0.4.  respectively. 


Figure  39.  Probability  of  Image  Overlaying  n  Pixels  with  7  =  0 
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The  equations  for  5n  and  pn  can  then  be  combined  to  determine  the  maximum 
received  signal,  E{Sr }.  Figure  32  in  Chapter  V  displays  the  result  for  various 
obscuration  ratios. 
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